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ABSTRACT 


The  aerothermoelastlc  characteristics  of  a tactical  missile  operating  at 
near  zero  angle  of  attack  In  the  low  hypersonic  speed  range  have  been  studied. 
An  Integrated  digital  computer  program  has  been  developed  for  the  aeroelastlc 
response  of  sharp-edge  aerodynamical ly  heated,  cantilevered  lifting  surfaces 
representative  of  tactical  missile  fins.  The  nonuniform  flow  generated  by 
a spherically  blunted  missile  nose  Is  accounted  for.  Based  upon  the  applica- 
tion of  the  computer  program  to  selected  cases,  recommendations  are  made  re- 
garding the  analysis  procedures.  It  Is  found  to  be  best  from  the  standpoint 
of  computer  economy  to  run  a full  nonlinear  quasl-statlc  analysis  for  each 
trajectory,  and  to  perform  flutter  analyses  only  at  discrete  times  along  the 
trajectory,  selected  by  the  analyst.  Appropriate  prestress  conditions  are  ob- 
tained from  the  quasl-statlc  solution.  In  a nonlinear  flutter  analysis  the 
solution  Is  obtained  by  direct  numerical  Integration  of  the  coupled  nonlinear 
equations.  However,  In  Its  present  form  this  solution  Is  markedly  more  time 
consuming  (on  the  computer)  than  standard  type  linear  flutter  analyses. 

An  analytic  theory  Is  developed  which  can  predict  the  nonunlfoxm  flow  effects 
of  nose  bluntness  In  the  complete  hypersonic  speed  range  N ^ 3.  The  predicted 
adverse  effects  of  missile  nose  bluntness  on  vehicle  djmamlcs  are  In  excellent 
agreement  with  available  experimental  results.  On  a tactical  missile  the  non- 
uniform  flow  field  generates  large  spanwlse  dynamic  pressure  gradients  over 
the  fins.  At  nonzero  angle  of  attack  this  Invlscld  shear  flow  makes  the  fin 
forces  displacement  sensitive.  Simple  viscous  perturbation  methods  are  de- 
veloped which  can  predict  the  viscous  Induced  effects  on  the  unsteady  fin 
aerodynamics  with  the  desired  accuracy. 

Thermos  true  tural  numerical  results  are  compared  with  existing  closed  form 
solutions  and  show  how  nonlinear  coupling  effects  cause  a deviation  from  the 
Idealized  cases  considered  In  closed  form  solutions.  Typical  results  Illus- 
trate the  thermally  Induced  stiffness  changes  In  wings  and  demonstrate  the 
Importance  of  boundary  conditions , temperature  sensitive  elastic  moduli.  Initial 
Imperfections,  planform  and  thickness  variations. 
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Section  1 
INTRODUCTION 


When  considering  the  aerothermoelastlc  characteristics  of  a finned  tactical 
missile  operating  in  the  low  hypersonic  speed  range,  which  is  the  subject  of 
the  present  study,  one  encounters  several  problem  areas.  The  aerodynamic 
heating  causes  a degradation  of  structural  properties  and,  in  combination  with 
thermally  induced  stresses,  accentuates  nonlinear  structural  behavior.  Aero- 
dynamically  the  usual  nonlinear  airfoil  characteristics  of  the  fins  are  com- 
plicated by  the  nonuniform  flow  field  generated  by  the  slight  nose  bluntness 
necessary  for  the  survival  of  the  missile  nose  cone.  The  goal  of  the  present 
study  is  to  combine  these  nonlinear  features  in  one  computer  program  that  will 
be  able  to  predict  divergence  and  flutter  speeds  for  the  fins  on  a tactical 
missile  flying  at  near  zero  angle  of  attack  at  Mach  numbers  between  <=  3 

and  M * 6. 

00 

One  of  the  features  that  distinguishes  the  aeroelastlc  behavior  of  lifting 
surfaces  at  subsonic  and  low  supersonic  speeds  is  that  isothermal,  linear 
theories  generally  suffice  to  define  both  structural  and  aerodynamic  behavior. 
There  are  exceptions  such  as  the  nonlinear  flow  behavior  at  transonic  speeds 
and  nonlinearities  at  all  speeds  caused  by  separated  flow.  In  most  cases, 
however,  the  equations  can  be  linearized  making  it  possible  to  conduct  the 
different  phases  of  aeroelastlc  analyses  independently.  Because  the  effects 
of  prestress  and  small  geometric  imperfections  are  ignored  in  the  lineariza- 
tion process,  both  strength  and  stiffness  are  Invariant  with  time.  This  leads 
to  a rather  simple  definition  of  critical  loading  conditions  by  means  of  the 
dynamic  pressure  and  the  velocity  normal  load  factor,  and  the  response  to 
gusts  and  turbulence  can  be  considered  separately.  Above  all,  it  is  the  re- 
moval of  time  as  a variable  in  the  definition  of  the  aerothermoelastlc  environ- 
ment which  results  in  a greatly  simplified  analysis. 
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Once  the  flight  regime  is  extended  to  high  supersonic  speeds,  aerodynamic 
heating  must  be  included  and  the  analysis  becomes  more  complicated.  The 
temperatures  vary  with  time,  and  since  thermal  stresses  change  the  structural 
stiffness,  the  aeroelastic  problem  cannot  be  formulated  in  terms  of  the  current 
flight  condition  even  if  one  assumes  linear  structural  behavior.  The  classic 
aeroelastic  triangle  of  interacting  Inertial,  Elastic,  and  Aerodynamic  forces 
used  by  Collar  (Ref.  1)  becomes  the  aeroelastic  rectangle  shown  by  Rogers 
(Ref.  2)  and  reproduced  here  as  Figure  1.1  with  only  the  strong  couplings 
indicated. 

In  the  mid-fifties  numerous  studies  were  initiated  on  the  influence  of  aero- 
dynamic heating  on  aeroelastic  behavior  in  the  low  supersonic  speed  range. 

Usually  these  studies  were  based  on  linearized  theories.  The  thermal,  aero- 
dynamic, or  structural  effects  of  the  problem  were  usually  considered  separately. 
Only  a few  investigations  were  devoted  to  their  interactions.  These  Included  a 
wide  ranging  study  on  static  phenomena  at  Bell  Aircraft  which  took  full  advan- 
tage of  the  digital  computer  capability  existing  at  that  time  (Refs.  3-7). 

With  the  aid  of  an  early  model  digital  computer,  significant  progress  was  also 
made  at  M.I.T.  toward  the  understanding  of  the  nonlinear  effects  of  thermally 
Induced  stresses  (Refs.  8-10).  In  addition  to  these  and  other  similar  studies 
with  very  limited  distribution,  numerous  closed  form  solutions  have  been  pub- 
lished. They  are  all  based  on  simplifying  assumptions  such  as  infinite  aspect 
ratio,  zero  chordwise  flexibility,  and  idealized  boundary  conditions.  Such 
solutions  are  satisfactory  for  demonstration  of  certain  phenomena  but  their 
usefulness  for  analysis  of  practical  flight  structures  is  limited.  Harder 
et  al  (Ref.  11)  used  a combination  of  digital  and  analog  methods  to  investigate 
the  flutter  of  heated  monocoque  wing  structures. 

With  the  extension  of  the  flight  regime  to  hypersonic  speeds  nonlinear  effects 
become  dominant.  Over  the  past  decade  many  studies  have  been  made  of  the 
aeroelastic  behavior  of  a wide  range  of  vehicle  configurations  operating  at 
hypersonic  speeds.  These  include  studies  of  static  aerodynamics  and  vehicle 
dynamics  of  reentry  vehicles  and  manned  spacecraft,  flutter  analyses  of 
lifting  surfaces  accounting  for  body  flexibility,  and  panel  flutter  analyses. 

By  the  early  nineteen  fifties,  the  nonlinear  relationship  between  the  local 
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angles  of  attack  of  a body  and  the  pressure  distribution  was  well  understood 
for  inviscid  flow.  The  Importance  of  thermal  stresses  was  appreciated.  It 
was  also  common  knowledge  that  for  thin  platelike  lifting  surfaces,  trans- 
verse deflections  of  the  same  order  as  the  thickness  result  in  nonlinear  force- 
deflection  relationships.  It  was  possible  to  evaluate  the  effects  of  each  of 
these  complications  separately,  but  the  state  of  technology  did  not  allow  the 
inclusion  of  all  facets  of  the  aeroelastic  problem  in  a realistic  analysis  of 
the  flight  vehicle. 

An  important  paper  by  Ashley  and  Zartarian  (Ref.  12)  demonstrates  the  use  of 
piston  theory  in  various  applications  to  lifting  surfaces  with  infinitely 
sharp  leading  edges.  Although  the  thickness  of  the  airfoil  has  been  shown  to 
have  significant  nonlinear  effects  on  the  pressure  distribution,  the  equations 
of  motion  can  still  be  linearized  and  the  classical  method  of  solving  linear 
equations  can  be  used.  This  is  done  for  simple  chordwise  rigid  sections  in  a 
series  of  MIT  reports  which  are  based  on  the  unsteady  shiick  expansion  method 
and  the  unsteady  Newtonian  theory  (Refs.  13-17).  It  was  found  that  for  a 
"typical  section"  with  sharp  leading  edge  the  effect  of  the  aerodynamic  non- 
linearity on  the  flutter  speed  was  small.  During  the  late  fifties  and  early 
sixties,  the  Bell  Aerospace  Corporation,  under  Air  Force  Systems  Command  sponsor- 
ship, produced  a series  of  reports  on  static  aero-thermoelastic  phenomena  at 
hypersonic  speeds  (Refs.  18-20).  In  Ref.  20,  geometrically  nonlinear  struc- 
tural behavior  and  thermal  stresses  are  accounted  for.  Various  inviscid  aero- 
dynamic theories  are  presented  in  simple  form,  such  as  tangent  wedge,  modified 
Newtonian,  shock  expansion,  and  blast  wave  theories.  Simple  means  to  correct 
for  leading  edge  sweep,  nose  bluntness,  and  boundary  layer  effects  are  provided. 

This  sampling  of  the  literature  has  been  made  to  show  what  types  of  analyses 
are  available.  They  are  all  for  highly  idealized  conditions  and  are,  therefore, 
insufficient  for  determination  of  the  aeroelastic  characteristics  of  a tactical 
missile  in  the  low  hypersonic  flight  regime  ((3  <:  5 6),  the  objective  of 

the  present  study.  Even  if  one  considers  only  the  fins,  the  aerodynamics  can 
become  rather  complicated  for  a typical  geometry  (see  Fig.  1.2).  A nonuniform 
inviscid  flow  field,  the  "entropy  wake",  is  generated  by  the  curved  bow  shock. 
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causing  the  vehicle  aerodynamics  to  become  highly  nonlinear.  The  entropy 
gradients  in  the  "entropy  wake"  have  large  adverse  effects  on  the  vehicle 
dynamics  due  to  the  convective  time  lag  associated  with  the  "entropy  wake" 
effects  (Refs.  21-23).  When  only  the  local  strip-aerodynamics  on  the  fin  are 
considered,  the  "entropy  wake"  effect  is  to  introduce  a spanwlse  dynamic  pres- 
sure variation  and,  at  non-zero  angle  of  attack,  a locally  nonuniform  flow  in 
the  "2-D  strip  planes"  (see  Fig.  1.2).  In  the  present  study  an  already  de- 
veloped hypersonic  theory  (Refs.  21-23)  is  extended  down  to  = 3. 

Viscous  effects  are  usually  important  in  hypersonic  flow.  The  large  effects 
of  boundary  layer  transition  on  reentry  vehicles  are  well  known  (Refs.  24  and 
25).  The  transition  effects  are  of  special  concern  when  local  temperature  is 
important,  such  as  is  the  case  on  an  ablating  vehicle  (Refs.  21  and  26)  or  on 
thermally  distorted  surfaces  (Ref.  25).  The  fin  shock  will  cause  the  fuselage 
boundary  layer  to  separate,  generating  complex  aero-thermodynamic  interactions 
(Ref.  27),  which  can  have  profound  effects  on  the  unsteady  aerodynamics  of  the 
fins.  Within  the  scope  of  the  present  study  only  an  exploratory  analysis  of 
the  viscous  flow  effects  is  possible. 

It  is  desirable  to  develop  an  aeroelastic  analysis  capability  for  the  hyper- 
sonic range  that  includes  all  important  nonlinear  effects.  Experience  indi- 
cates that  below  hypersonic  speeds  mechanically  induced  stress  can  be  neglected, 
as  the  geometric  changes  are  small  for  flight  velocities  significantly  below 
those  that  cause  aeroelastic  instability.  This  makes  the  structural  and  aero- 
elastic  analyses  relatively  Independent  of  each  other,  contrary  to  what  is  the 
cane  at  hypersonic  speeds. 

The  nonlinear  coupling  between  structural  deformation  and  aero-thermodynamic 
behavior  is  probably  the  most  important  interaction  which  can  be  neglected 
only  at  subsonic  and  low  supersonic  speeds.  For  example,  in  a symmetric  air- 
foil section  the  inplane  stresses  do  not  produce  any  out-of-plane  distortions, 
and  hence  no  change  in  distribution  of  aerodynamic  pressures,  until  thermal 
buckling  or  aeroelastic  instability  is  encountered.  In  practical  applications 
a number  of  complications  arise.  Syu^netric  profiles  with  imperfections,  as 
well  as  nonsymmetric , cambered,  or  deliberately  warped  sections,  experience 
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thermally  Induced  lateral  deformations,  which  cause  the  lift  distribution  to 
become  dependent  upon  the  vehicle  flight  history.  The  fact  that  any  plate~ 
like  wing  subjected  to  bending  will  exhibit  chordwlse  curvature  due  to  the 
Poisson  effect  also  provides  thermoelastic  coupling  not  considered  In  linear 
theory. 

Since  In  the  hypersonic  domain  larger  deformations  must  be  allowed,  the  non- 
linear Interdependence  between  deformation,  aerodynamic  forces,  and  temperatures 
becomes  continually  stronger  as  the  Mach  number  Is  Increased.  This  makes  It 
necessary  to  Include  thermal,  aerodynamic,  and  structural  analyses  In  the  same 
program  for  evaluation  of  aeroelastlc  stability.  Here  such  a program  Is  de- 
veloped for  the  analysis  of  a lifting  surface  on  a tactical  missile.  The 
computer  program  Is  built  around  the  STAGS  code  (Ref.  28),  which  Includes 
geometrically  nonlinear  structural  effects.  A modified  STAGS  code  constitutes 
the  heart  of  the  program  package,  with  the  other  programs  attached  to  ^t.  In 
what  follows,  the  analyses  leading  to  the  different  computer  codes  used  in  the 
Integrated  aerothermoelastlc  computer  program  are  described  first.  After  that 
the  numerical  results  obtained  with  the  Integrated  program  are  discussed  and 
compared  with  the  limited  amount  of  pertinent  data  that  Is  available.  Using 
these  results,  an  optimum  procedure  Is  formulated  for  aeroelastlc  analysis  of 
a hypersonic  tactical  missile,  based  upon  current  state  of  the  art.  Recom- 
mendations are  made  for  further  Improvement  of  this  analytic  capability.  In- 
structions for  the  use  of  the  aerothermoelastlc  computer  program  are  given 
In  Appendix  B. 
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Fig.  1.1  A'^rcthermoelastic  Rectangle  (from  Ref.  2) 

Arrows  denote  strong  coupling. 
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Figure  1.  2 Hypersonic  missile  flow  Reid 
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Section  2 

STRUCTURAL  ANALYSIS 


The  thermoelastic  analysis  In  the  computer  program  Is  carried  out  by  use  of  a 
combination  of  the  STAGS  program  for  elastic  analysis  and  a program  for  de- 
termination of  the  temperature  history  In  the  structure.  These  two  codes  are 
briefly  described  in  the  following.  Specific  numerical  results  are  presented  to 
verify  that  these  codes  are  suitable  for  aeroelastlc  analysis  of  tactical  missiles . 

2.1  The  STAGS  Code 

STAGS  is  Intended  for  analysis  of  shell  or  plate  type  structures.  The  struc- 
ture to  be  considered  may  consist  of  up  to  30  different  shell  branches,  each 
treated  by  use  of  finite  difference  discretization.  In  addition,  the  struc- 
ture may  Include  some  finite  elements:  elastic  bars,  shear  panels,  beams, 

and  nonlinear  plate  elements.  The  computer  code  is  discussed  in  detail  in 
Ref.  28.  Only  a brief  description  is  given  here. 

The  analysis  In  STAGS  Is  based  on  a finite  difference  energy  method.  That  Is, 
the  total  potential  energy  Is  expressed  in  terms  of  the  applied  mechanical  or 
thermal  loads  and  the  displacement  components  at  a number  of  node  points.  This 
is  possible  after  derivatives  occurring  In  the  expression  for  the  strain  energy 
have  been  substituted  by  their  finite  difference  equivalents  and  the  energy 
has  been  numerically  Integrated  over  the  entire  structure.  The  total  potential 
energy  Is  minimized  with  respect  to  the  degrees  of  freedom  of  the  system 

dv/ax^  = a/aXj^(u-Q)  (2.1) 

Here  U is  the  strain  energy,  0 the  work  done  by  external  (or  thermal) 
forces,  and  the  represent  the  displacement  unknowns.  The  minimization 

leads  to  a set  of  equations  of  the  form 


K(X)  = fF] 


(2.2) 


where  K is  a nonlinear  algebraic  operator  (of  the  third  order).  The  force 
vector  [ F}  consists  of  thermal  and  mechanical  loads,  and  also  contains  the 
strains  caused  by  inelastic  deformation.  The  equilibrium  equations  for  a 
system  in  motion  are  obtained  through  the  addition  of  the  effects  of  inertia 
and  damping. 


[M]fX}  + [D]fX}  + K(X)  = fF}  (2.3) 

Options  are  given  in  the  program  that  allow  for  omission  or  inclusion  of  the 
dynamic  terms.  The  perturbation  technique  is  used  to  formulate  the  eigen- 
value problems  for  bifurcation  buckling  and  small  amplitude  vibrations. 

Consequently,  the  code  can  be  used  for: 
o Linear  stress  analysis 

o Geometrically  nonlinear  elastic  stress  analysis 

o Inelastic  stress  analysis,  geometrically  linear  or  nonlinear 

o Bifurcation  buckling  analysis  with  linear  or  geometrically 
nonlinear  prestress  (elastic) 

o Small  amplitude  vibration  analysis  with  prestress  based  on 
linear  or  geometrically  nonlinear  analysis  (elastic) 

o Transient  response  analysis,  linear  or  geometrically  non- 
linear, elastic  or  inelastic 

Any  combination  of  point  forces,  line  loads,  and  distributed  surface  tractions 
can  be  applied.  Loading  by  specification  of  displacements  or  thermal  gradients 
(through  the  shell  wall  and  over  the  shell  surface)  is  also  permitted.  Any  set 
of  boundary  conditions  or  other  displacement  constraints  can  be  included. 

Restart  capability  exists  both  in  the  nonlinear  static  analysis  and  in  the 
transient  analysis.  Intermediate  data  can  be  saved  on  tape  or  file,  permitting 
the  analysis  to  be  continued  a later  time.  The  input  data  is  automatically 
checked  for  certain  types  of  errors  and  for  inconsistencies.  If  so  desired, 
the  user  can  suppress  execution  of  the  program  in  order  to  obtain  a check  on 
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the  input.  The  output  from  such  a preliminary  run  may  include  a graphic  pre- 
sentation of  the  shell  geometry  with  a display  of  the  grid  lines.  The  results 
of  a STAGS  run  can  be  transferred  to  tape  or  disk  file  to  be  used  as  input  for 
a postprocessor  to  obtain  contour-  or  cross-plots  of  displacements,  stress  re- 
sultants, and  stresses.  Displacement  histories  can  also  be  oloi.ted  in  the 
case  of  transient  or  nonlinear  static  analysis. 

The  STAGS  code  accepts  the  input  in  two  different  forms,  regular  data  cards 
and  user  written  subroutines.  The  latter  type  of  input  may  sometliaes  be  more 
demanding  on  the  user.  However,  it  makes  possible  the  definition  of  input 
parameters  such  as  loads,  temperatures,  etc.  by  functional  relationships. 

This  tends  to  increase  the  generality  of  the  computer  code  and  also,  in  many 
cases,  results  in  a more  compact  input  deck.  User-written  subroutines  can 
be  used  for 

o Initial  imperfections 
o Loads 

o Load  factor  history  (for  transient  analysis) 

0 Shell  wall  properties  and  temperatures 
o Reference  surface  geometry 
o Grid  generation 

o Boundary  conditions  and  other  constraints 
o Stiffener  geometry  and  temperature 
o External  damping 

The  shell  configuration  is  described  by  use  of  a reference  surface.  The  shell 
itself  can  be  offset  from  this  reference  surface  by  a small  amount  that  is 
allowed  to  vary  ever  the  surface.  The  location  of  any  point  on  the  reference 
surface  can  be  uniquely  defined  by  use  of  a set  of  two  independent  parameters, 
X,  Y.  The  geometry  of  the  surface  is  given  after  a functional  relationship 
has  been  established  that  defines  the  coordinates  in  a Cartesian  system  x,y,z 
as  functions  of  the  values  of  these  two  parameters. 
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X = x(X,  Y),  y * y(X,  Y)  , z = x(X,  Y) 


(2.4) 


The  parameters  X and  Y defining  the  position  of  a point  on  the  surface  are 
referred  to  as  surface  coordinates.  For  shell  geometries  not  included  as 
standards  in  STAGS  the  user  must  define  the  relation  given  in  Eq.  (2.4)  in 
a user-written  subroutine. 

If  one  of  the  surface  coordinates  is  held  constant  and  the  other  is  varied 
over  its  range,  a trace  is  obtained  on  the  surface.  Such  a trace  is  referred 
to  as  a coordinate  line.  Coordinate  lines  have  a two-fold  purpose  in  the 
analysis.  The  directions  of  the  in-plane  displacements  are  defined  by  the 
coordinate  lines.  In  addition,  a set  of  coordinate  lines  are  used  to  define 
the  surface  grid  that  is  needed  in  a finite  difference  analysis.  In  most 
cases  the  two  sets  of  coordinate  systems  are  identical,  but  sometimes  the 
natural  coordinates  for  the  shell  are  not  suitable  for  grid  generation.  In 
that  case  a separate  set  of  surface  coordinates,  denoted  E,  T|,  can  be  used 
to  define  the  grid  through  a user-written  subroutine.  It  is  also  possible  to 
define  directly  the  X,  Y coordinates  for  the  grid  points,  without  use  of  an 
additional  coordinate  system,  either  by  use  of  regular  input  data  cards  or  by 
a user-written  subroutine. 

For  definition  of  shell  geometry  a Cartesian  system  x,y,z  (see  Eq.  (2.4)) 

is  defined  for  each  shell  branch.  One  Cartesian  system  x , y , z is  chosen 

SEE 

as  a global  coordinate  system.  Orientation  parameter  cards  define  the  orien- 
tation of  the  individual  coordinate  systems  of  each  of  the  branches  with  re- 
spect to  the  global  system.  The  global  system  is  not  defined  explicitly  by 
input  cards  but  is  given  implicitly  by  the  values  of  the  orientation  parameters. 

Stiffeners  are  defined  as  structural  elements  attached  to  the  shell  surface. 
Effects  of  stiffener  cross  section  warping  or  deformation  are  neglected. 
Stiffeners  can  either  be  considered  as  discrete,  in  which  case  they  are  de- 
fined one  by  one,  or  their  contribution  to  the  shell  wall  stiffness  can  be 
"smeared"  over  the  shell  surface.  Smeared  stiffeners  are  used  for  r mvenience 
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(to  reduce  the  number  of  data  cards)  or  in  order  to  suppress  the  local  de- 
formation between  the  stiffeners.  Present  constraints  in  the  computer  program 
are  as  follows:  inelastic  analysis  can  only  be  used  with  an  initially  iso- 

tropic material.  It  cannot  be  used  with  material  properties  that  vary  with 
the  surface  coordinates.  The  yield  stress  must  be  independent  of  temperature. 
Bifurcation  buckling  and  small  amplitude  vibration  analyses  cannot  include 
plasticity. 

2.2  The  TRRID  Program 

For  computation  of  temperatures  in  the  shell  the  code  TRRID  developed  by 
Compton  and  Schultz  (Refs.  29,30)  computes  the  transient  thermal  response  of  a 
body  subject  to  arbitrary  aerodynamic  heating.  The  shell  or  plate  under  con- 
sideration is  divided  into  a number  of  layers  (called  nodes).  Individual 
layers  can  be  of  different  materials.  Temperature  varient  thermal  properties 
can  be  included  for  all  layers.  The  mid-surface  of  the  geometry  is  assumed 
to  be  adiabatic,  and  the  exterior  surface  is  subjected  to  time-varient  aero- 
dynamic heating.  Heat  conduction  into  the  body  is  one-dimensional,  and  radia- 
tion from  the  surface  to  the  surroundings  can  be  taken  into  account.  The 
finite  difference  equations  are  derivedby  conventional  energy  balance  con- 
siderations for  each  layer  using  an  implicit  difierence  scheme.  The  tempera- 
ture response  is  determined  by  a tri-diagonal  macrlx  solution  at  each  time 
step. 

The  calculation  of  the  aerodynamic  heat  transfer  coefficient  is  performed  at 
each  time  step  based  on  updated  ambient  conditions.  In  summary,  the  aero- 
thermodynamic  techniques  used  are : 

1.  Stagnation  point  - Fay-Riddel  method. 

2.  Flat  plate  - Blaslus  method  for  laminar  flow  and  the 

Spalding-Chi  method  for  turbulent  flow. 

3.  Conical  flow  - Flat-plate  methods  modified  for  conical 

flow  as  given  by  VanDriest  in  Ref.  31. 

Transition  from  laminar  to  turbulent  flow  is  assumed  to  occur  Instantaneously 
at  a boundary  layer  edge  Reynold's  number  of  500,000,  which  is  very  much  on  the 
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conservative  side  for  the  Mach  number  range  of  interest.  No  effort  is  made 
to  compute  a transition  zone  or  to  compute  a virtual  origin  for  the  turbulent 
boundary  layer.  It  is  highly  desirable  that  a code  be  developed  which  gives 
a better  estimate  of  the  temperature  field. 


2 . 3 Code  Verification 

A few  examples  follow  here  in  which  results  from  STAGS  are  compared  to  re- 
sults obtained  by  other  means.  The  first  of  these  excludes  temperature  effects. 
The  next  series  of  examples  includes  a temperature  field  defined  in  the  STAGS 
input  as  varying  parabolically  in  the  chordwise  direction.  The  last  few 
examples  are  obtained  by  use  of  STAGS  and  TRRID.  The  temperature  field  is 
determined  internally  and  thermal  degradation  of  material  properties  is  in- 
cluded. 

2,3.1  Small  Amplitude  Vibrations  of  Stress  Free  Plate 

An  important  function  of  STAGS  in  the  aeroelastic  analysis  is  to  determine 
modes  and  frequencies  of  small  amplitude  vibration.  The  capability  to  perform 
such  analyses  has  been  included  recently.  In  order  to  verify  this  new  capa- 
bility, STAGS  results  were  obtained  and  compared  with  results  from  FLUENC, 
which  is  a sub-program  of  the  MOFA  system  of  computer  programs  (Ref.  32). 

The  analysis  was  performed  for  a plate  of  rectangular  planform  with  a 10-in. 
chord  and  an  8-in.  span  as  shown  in  Fig.  2.1.  The  grid  has  11  node  points  along 
the  chord  and  5 along  the  span.  The  plate  was  clamped  over  a portion  of  the 
root  chord  as  shown  in  the  figure.  The  grid  is  rather  coarse  and  the  results 
will  not  be  very  accurate.  However,  it  has  been  used  throughout  the  program 
development  for  reasons  of  computer  economy.  It  is  sufficiently  fine  to  indi- 
cate the  gross  behavior  of  the  plate.  As  with  decreasing  grid  size  the  STAGS 
results  converge  from  below  and  the  FLUENC  results  from  above  we  can  expect 
some  discrepancy.  First  STAGS  and  FLUENC  were  applied  to  a plate  with  constant 
thickness,  0.2  in.  The  mode  shapes  from  the  two  programs  were  very  close. 

The  frequencies  corresponding  to  this  coarse  grid  are  shown  in  Table  2.1. 
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Table  2.1 


Comparison  of  Frequencies  for  Constant  Thickness  Plate 


Mode 

FREQUENCY  (Hz) 

Number 

FLUENC 

STAGS 

STAGS 

Grid  (5x11) 

Grid  (5x11) 

Refined  Grid  (15x21) 

1 

63.65 

57.12 

61.67 

2 

122.14 

106.29 

111.82 

3 

287.41 

266.84 

281.29 

4 

523.27 

466.10 

498.76 

The  relatively  good  agreement  between  the  results  Indicates  that  with  diminish- 
ing grid  size  they  would  agree  In  the  limit.  This  Impression  Is  reinforced  by 
the  fact  that  with  a refined  grid  (3  times  denser  In  the  most  critical  span- 
wise  direction)  STAGS  gives  the  results  shown  In  the  last  column  of  Table  2.1. 

In  another  example  a plate  with  variable  thickness  is  considered.  The  chord- 
wise  section  has  a "diamond  profile"  (double  wedge).  The  section  properties, 
Independent  of  the  spanwlse  coordinate,  correspond  to  0.05-ln.  thickness  at 
the  leading  and  trailing  edges  and  a maximum  thickness  of  0.45-ln.  at  mid- 
chord. The  same  grid  was  used  as  In  the  previous  case.  We  may  expect  some- 
what larger  discrepancy  In  this  case  since  the  plate  Is  very  thin  at  the 
leading  and  the  trailing  edges,  and  thus  a finer  spacing  Is  required  for  con- 
vergence. Results  are  shown  In  Table  2,2.  The  last  column  shows  results  ob- 
tained by  use  of  a grid  with  variable  spacing.  These 

results  verify  the  validity  of  the  program  and  underscore  the  usefulness  of 
a variable  spacing  -.  rid  for  this  type  of  plate.  It  shoud  be  noticed  that 
with  constant  thickness  elements  in  FLUENC  we  expect  the  frequencies  to  be 
somewhat  high. 


Table  2.2 


Mode 

Number 

FREQUENCY  (Hz) 

FLUENC 
Grid  (5x11) 

STAGS 

Grid  (5x11) 

STAGS 

Grid  (9x11) 

1 

125.99 

103.40 

117.97 

2 

272.47 

237.07 

238.87 

3 

644.58 

570.95 

604.02 

4 

1044.34 

782.05 

798.10 

2.3.2  Influence  of  Temperature  Gradients 

For  minimum  drag  many  plate-like  wing  structures  subjected  to  hypersonic  flow 
are  designed  so  that  they  are  very  thin  at  the  leading  and  trailing  edges 
(the  diamond  shaped  profile).  This  results  in  temperature  distributions  with 
markedly  higher  temperatures  at  the  edges  and  a minimum  at  mid-chord.  The 
thermal  expansion  results  in  that  case  in  spanwise  compression  at  the  edges 
and  tension  in  the  cooler  midsection.  Such  temperature  distributions  can 
have  a profound  effect  on  the  aeroelastic  behavior  of  the  wing.  Consequently, 
the  effect  has  been  studied  in  the  past  and  several  research  papers  have  been 
published  (Refs.  33-40).  These  papers  were  published  before  the  h ^h  speed 
computer  had  made  detailed  analysis  of  more  complex  configurations  possible. 

In  the  following  we  will  compare  results  obtained  from  STAGS  with  the  closed 
form  solutions  of  Ref.  34  for  constant  thickness  plates.  The  temperature  dis- 
tribution, however,  will  be  one  that  is  typical  of  a wing  with  a diamond 
shaped  profile. 

2.3.3  Thermal  Buckling 

With  the  temperature  distribution  discussed  above,  compressive  stresses  at 
the  edges  are  increasing  with  the  temperature.  At  some  temperature  level 
thermal  buckling  will  occur.  According  to  Ref.  34,  the  critical  value  of 
the  difference  between  edge  and  mid-chord  temperature  is 
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(2.5) 


AT 


crit 


18.5 


D 

E O'  t c2 


(g  + n (g  + 3) 


For  parabolic  temperature  variation,  with  a minimum  at  mid-chord,  E = 2. 

In  that  case,  a square  aluminum  plate  with  a 20  in.  side  and  0.5  in.  thick- 
ness will  buckle  at  “ 153°F.  With  a discretized  analytical  model, 

as  in  the  STAGS  code,  the  results  are  accurate  only  if  a sufficiently  fine 
grid  spacing  is  used.  Therefore,  for  comparison  with  the  analytical  results 
in  Ref.  3A  STAGS  results  were  obtained  for  various  grid  sizes.  With  a 
relatively  coarse  grid  spacing,  2.0  in.,  the  critical  temperature  difference 
is  118°F,  far  below  the  analytical  results,  153°F.  The  primary  reason  for 
the  discrepancy  is  that  the  temperature  gradient  is  large  at  leading  and 
trailing  edges  and  the  thermal  stresses  are  therefore  not  well  represented 
in  the  prebuckling  solution  with  this  coarse  grid  size.  As  the  grid  is  re- 
fined, values  of  ^T  closer  to  the  value  of  Ref.  34  are  obtained.  Thus, 

crit 

with  a grid  spacing  of  1.0  in.,  ” 142°F. 

With  a sufficiently  fine  grid  the  STAGS  results  should  be  more  accurate  than 

those  of  Ref.  34.  The  reason  for  this  is  that  Eq.  (2.5)  is  based  on  the 

assumption  that  in  the  buckling  mode  the  chord  remains  straight.  The  STAGS 

results  show  some  chordwlse  bending  in  the  buckling  mode.  If  essentially 

the  same  restriction  is  introduced  in  STAGS  (by  increasing  the  chordwise 

A o 

bending  stiffness  by  a factor  of  10  ) ^"^crit  F.  Thus, 

the  validity  of  the  STAGS  code  is  verified  for  this  type  of  application. 

The  model  with  coarse  grid  spacing,  2.0  in.,  although  somewhat  inaccurate, 
reflects  correctly  the  basic  structural  behavior.  In  the  Interest  of  computer 
economy,  subsequent  results  were  obtained  with  this  grid  size. 

The  temperature  field  under  consideration  here  is  symmetric  about  the  mid- 
plane of  the  plate.  . Consequently,  for  a perfectly  flat  plate  there  will  be 
no  out  of  plane  deformation  in  the  prebuckling  stage,  AT  < Above 

the  critical  temperature  this  equilibrium  configuration  becomes  unstable. 

With  the  temperature  difference  increasing  above  a new  deformation 

mode  gradually  develops.  This  mode  is  essentially  represented  by  pure  twist- 
ing of  the  wing. 
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The  value  of  the  twist  at  the  wing  tip  is  shown  as  a function  of  the  tempera- 
ture in  Fig.  2.2.  The  temperature  difference  is  normalized  with  respect  to 

so  that  thermal  buckling  occurs  at  X = 1.0.  The  twist  parameter  0 
is  defined  in  Ref.  34.  The  results  for  the  perfectly  flate  plate,  as  obtained 
in  Ref.  34,  are  shown  by  the  curve  marked  0^  = 0. 

If  the  plate  contains  an  imperfection  in  shape  that  is  not  orthogonal  to  the 
buckling  no  equilibrium  can  exist  without  lateral  displacements.  By 

use  of  a nonlinear  analysis  it  is  possible  to  determine  the  growth  of  this 
lateral  displacement  with  increasing  values  of  X . An  initial  imperfection 
in  the  form  of  a uniform  twiot  of  the  wing  was  introduced  in  Ref.  34.  The 
amplitude  of  the  twisting  mode  v^as  then  determined  as  a function  of  X.  Re- 
sults were  obtained  for  different  values  of  the  initial  twist  as  shown  in 
Fig.  2.2. 

For  two  values  of  the  initial  twist,  0 = 0.05  and  0.20,  results  were  ob- 
tained by  use  of  STAGS  and  are  also  shown  in  Fig.  2.2.  At  any  finite  value 
of  X the  STAGS  results  indicate  a lower  value  of  the  twist.  It  should  be 
noted  that  the  curves  with  STAGS  results  are  normalized  with  respect  to  a 
lower  value  of  STAGS  results  normalized  with  respect  to  the 

critical  load  of  Ref.  34  would  indicate  higher  values  of  the  twist  than  those 
shown  in  Ref.  34.  This  reflects  the  fact  that  the  discrete  model  is  weaker 
due  to  the  coarser  spacing  and  also  due  to  the  inclusion  of  chordwise 
bending.  The  occurrence  of  chordwise  bending  was  pointed  out  in  Ref.  35. 

The  twist,  defined  as  the  rotation  of  an  element  of  the  chord  around  a span- 
wise  axis,  is  constant  in  the  analysis  of  Ref.  34,  but  is  allowed  to  vary 
along  the  chord  in  the  STAGS  analysis.  STAGS  results  for  the  local  twist 
are  shown  in  Fig.  2 3 . 

The  variation  of  the  torsional  stiffness  of  the  lifting  surface  with  tempera- 
ture and  initial  twist  is  also  discussed  in  Ref.  34.  STAGS  results  for  the 
torsional  stiffness  are  compared  to  those  from  Ref.  34  in  Figure  2.4.  The 
agreement  between  results  from  the  two  methods  provides  additional  verifica- 
tion of  the  STAGS  code.  For  a perfect  plate  the  torsional  stiffness  (or  the 
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square  of  the  corresponding  vibration  frequency)  decreases  linearly  with 
increasing  temperature,  so  that  it  becomes  equal  to  zero  at  X 1.0.  The 
imperfect  plate  retains  a finite  stiffness.  The  explanation  of  the  fact 
that  the  imperfect  shell  retains  a finite  stiffness  at  X 1.0  is  that  the 
edges  are  stretching  as  the  wing  twists.  This  reduces  the  temperature  induced 
spanwise  compression  stress. 

Less  information  is  available  on  the  effect  of  the  temperature  distribution 
on  the  bending  stiffness  of  the  plate.  Kochanski  and  Argyris  (Ref.  36)  and 
Mansfield  (Ref.  37)  consider  the  problem  for  wings  of  infinite  aspect  ratio. 
While  those  results  do  not  apply  to  the  low  aspect  ratio  wing  considered  in 
the  STAGS  analysis,  it  is  still  possible  to  compare  certain  trends.  It  is 
explained  in  Ref.  37  how  a temperature  field  with  higher  temperatures  at  the 
leading  and  trailing  edges  reduces  the  bending  stiffness  of  the  plate.  When 
the  plate  is  subjected  to  spanwise  bending,  the  Poisson  effect  results  in 
chordwise  bending  as  well.  Due  to  this  chordwise  bending,  the  higher  tempera- 
tures are  on  one  side  of  the  neutral  axis  and  the  lower  temperatures  on  the 
other  side.  This  introduces  "thermal  bending"  of  the  plate  in  the  same  direc- 
tion as  the  applied  "mechanical  bending".  Hence,  with  increasing  temperature 
the  plate  bending  stiffness  should  be  decreasing.  This  trend  is  shown  by  the 
STAGS  results  presented  in  Figure  2.5.  Again,  the  imperfection  is  in  the 
twisting  mode. 

The  ratio  between  the  fundamental  bending  and  torsional  frequencies  is  one  of 
the  important  parameters  in  flutter  analysis.  This  ratio  can  readily  be  ob- 
tained from  the  results  of  the  stiffness  calculations  with  STAGS  discussed 
above.  Results  are  presented  in  Figure  2.6.  As  is  to  be  expected,  the  fre- 
quency ratio  changes  markedly  with  the  temperature  and  its  behavior  is 
strongly  dependent  on  the  initial  imperfection  amplitude.  For  a perfectly 
flat  plate  the  frequency  ratio  becomes  infinite  at  X = 1.0  because  the  fre- 
quency of  the  torsional  mode  vanishes. 

In  conclusion,  STAGS  numerical  results  agree  well  with  the  analytical  results 
of  Ref.  34.  Whenever  differences  occur,  these  can  be  traced  to  the  restric- 
tion imposed  in  the  analysis  of  Ref.  34  in  order  to  make  a closed  form  solu- 
tion possible. 
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2. 3. A Lifting  Surfaces  with  Diamond  Shaped  Profile 

It  was  demonstrated  above  how  the  thermal  stresses  Induced  by  aerodynamic 
heating  dramatically  reduce  the  torsional  stiffness  of  the  lifting  surface. 
The  example  discussed  Illustrates  a basic  trend,  although  It  Is  somewhat 
less  than  realistic,  since  a plate  of  uniform  thickness  Is  used  In  the 
elastic  analysis  with  a temperature  distribution  typical  of  a surface  with 
a diamond  profile. 


In  other  attempts  at  analysis  (Refs.  33,  35-40))  the  variable  thickness  was 
accounted  for  In  the  elastic  analysis.  These  publications  consider  the  effect 
on  the  torsional  stiffness  of  spanwlse  thermal  stresses  In  wings  of  rectangu- 
lar planform.  The  wing,  treated  as  a simple  beam.  Is  clamped  at  the  root 
chord  and  generally  considered  to  have  an  Infinite  aspect  ratio.  Singer 
(Ref.  39)  Is  among  those  who  Introduces  an  approximate  correction  for  finite 
aspect  ratio. 


Material  property  degradation  with  the  temperature  is  another  complication 
that  cannot  easily  be  accounted  for  in  the  analytic  solutions.  Only  through 
use  of  a numerical  solution  procedure  is  it  possible  to  obtain  reasonably 
accurate  results.  Nevertheless,  the  analytic  solutions  are  useful  since 
they  give  an  indication  of  what  trends  to  expect  in  a computer  solution  of 
the  governing  partial  differential  equations. 


With  thin  leading  and  trailing  edges  on  the  lifting  surface  It  appears  neces- 
sary also  to  consider  local  deformation  modes.  Mansfield  (Ref.  40)  presents 
a closed  form  solution  for  the  thermal  bifurcation  buckling  of  a lifting 
surface  with  a sharp  leading  edge.  The  critical  load  In  that  case  depends 
only  on  the  temperature  gradient  at  the  edge,  the  thermal  expansion  coeffi- 
cient O'  , and  the  wedge  angle  (2*3^,).  A parabolic  temperature  distribution 
discussed  earlier  is  assumed.  The  critical  temperature  at  the  sharp  edge  is 
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This  simple  closed  form  solution  cannot  take  Into  account  the  effects  of 
property  degradation,  nose  bluntness,  finite  aspect  ratio,  or  boundary  condi- 
tions at  the  wing  root.  Neither  Is  It  possible  to  evaluate  wing  behavior 
when  the  critical  temperature  Is  exceeded.  However,  the  formula  should  be 
useful  for  a first  estimate. 

The  STAGS  program  with  the  thermal  analyzer  TRRID  was  used  to  study  the  be- 
havior of  rectangular  wings  with  diamond  shaped  cross-section  under  hyper- 
sonic flight  conditions.  The  cross-section  of  the  wing  Is  shown  In  Fig.  2.2. 
The  edges  are  somewhat  blunted  so  that  the  thickness  varies  from  0.05  In.  at 
the  leading  and  trailing  edges  to  0.45  In.  at  mid-chord.  For  comparison  with 
the  analytical  solutions  (Ref.  40)  the  spatial  temperature  distribution 
should  be  held  constant  In  time  also  In  the  STAGS  analysis.  The  temperature 
distribution  shown  In  Fig.  2.7  Is  typical  for  results  obtained  by  use  of 
TRRID.  The  four  lowest  vibration  frequencies  and  corresponding  modes  were 
determined  for  different  values  of  the  amplitude  of  this  distribution.  The 
results  are  shown  In  Fig.  2.7  both  for  a wing  clamped  along  the  entire  length 
of  the  root  chord  and  for  one  In  which  the  displacements  are  constrained  only 
at  the  node  points  7,8,9. 

For  the  wing  with  complete  root  chord  constraint  the  trends  of  the  lowest  fre 
quencles  agree  with  that  of  the  analytic  solutions  at  the  lower  edge  tempera- 
tures. The  torsional  stiffness,  represented  by  the  square  of  the  frequency 
of  the  second  mode,  decreases  with  the  temperature.  Extrapolation 

of  the  results  obtained  Indicate  that  the  stiffness  would  vanish  at  a tempera 
ture  of  about  680°F.  Analytic  solutions,  valid  only  for  longer  wings,  Indi- 
cate a higher  critical  temperature  of  about  760°F.  The  lowest  frequency, 
corresponding  to  the  first  bending  mode.  Is  less  sensitive  to  the  temperature 
The  elementary  beam  theory  does  not  predict  any  Influence  of  thermal  stress 
on  the  bending  stiffness.  However,  higher  order  theories,  such  as  those  by 
Mansfield  (Ref.  37)  or  Kochanskl  and  Argyrls  (Ref.  36)  show  the  same  trend 
as  the  STAGS  solution. 

The  frequencies  of  the  wing  with  only  partial  root  chord  constraint  shows 
very  little  dependence  on  the  temperature.  This  result  Is  to  be  expected 
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since  the  thermal  stresses  at  the  leading  and  trailing  edges  will  be  less 
significant  if  points  on  the  root  chord  are  free  to  move  in  the  spanwise 
direction.  Thermal  stresses  computed  for  the  rectangular  wing  are  shown  for 
three  different  boundary  conditions  in  Fig.  2.8.  The  results  in  Fig.  2.7 
also  include  the  reduction  in  stiffness  that  is  due  to  material  property 
degradation.  Since  the  thicker  midsection  of  the  wing  stays  relatively  cool, 
the  reduction  is  moderate. 

The  frequencies  of  the  third  and  fourth  modes  are  also  shown  in  Fig.  2.7  as 
functions  of  temperature.  At  all  temperatures  at  which  results  were  obtained, 
the  mode  corresponding  to  the  four  lowest  eigenvalues  are  clearly  identified 
as  the  first  and  second  bending  modes,  and  the  first  and  second  torsion  modes. 
However,  at  a temperature  of  690°F  the  solution  indicated  in  addition  one 
negative  eigenvalue.  This,  of  course,  represents  the  mode  corresponding  to 
leading  edge  buckling.  The  critical  load  then  is  exceeded  at  T * 690°F  but 
not  at  T = 460°F  (the  second  highest  temperature  for  which  results  were  obtained) 

It  is  clear  that  the  5x11  finite  difference  grid  used  is  too  coarse  to  give 
a good  estimate  of  the  critical  load  when  the  buckling  mode  is  localized. 

With  1.0  in.  chordwise  grid  spacing  the  critical  temperature  for  bifurcation 
buckling  is  found  to  be  570°F.  However,  examination  of  the  buckling  mode 
indicates  that  a much  finer  grid  is  needed  near  the  leading  edge.  A grid 
with  variable  spacing  was  introduced  (Fig.  2.9),  and  wit;h  this  grid  the 
critical  leading  edge  temperature  was  found  to  be  370°F.  The  corresponding 
buckling  mode  is  shown  in  Fig.  2.10.  Beyond  the  5th  row,  0.75  in.  from 

the  leading  edge,  the  lateral  displacements  are  very  small.  According  to 
Eq.  (2.6)  with  P = 0.08,  and  a sharp  edge,  the  critical  temperature  is  230°F. 

The  behavior  of  the  plate  in  the  postcritical  temperature  range  was  also  de- 
termined. The  nonlinear  analysis  does  not  distinguish  between  stable  and 
unstable  equilibrium  configurations.  True  bifurcation  modes  are  orthogonal 
to  the  precritical  deformation  pattern.  Therefore,  if  buckled  configurations 
are  to  be  computed,  it  is  necessary  to  introduce  a small  disturbance  as  a 
trigger.  In  this  case  a couple  of  lateral  loads  of  one  pound  each  and  with 
opposite  direction  were  applied  at  the  leading  edge  near  the  root  chord  (rows 
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3 and  4).  It  Is  found  that  the  amplitude  of  the  buckle  grows  gradually  with 
the  temperature,  so  that  the  deformation  still  Is  very  moderate  at  tempera- 
tures well  beyond  the  critical.  Modulus  degradation  with  the  temperature  does 
not  have  a detrimental  effect.  In  fact,  the  Increase  In  the  buckle  amplitude 
Is  somewhat  retarded  through  property  degradation,  probably  because  a reduced 
shear  modulus  allows  the  leading  edge  to  expand  mrr'  freely  in  comparison  to 
the  cooler  material  In  the  Interior.  It  Is  also  found  that  a leading  edge 
buckle  has  little  effect  on  the  overall  stiffness  of  the  wing. 

As  an  additional  check  on  the  overall  behavior  of  this  wing  in  response  to 
thermal  gradients,  the  analysis  was  repeated  for  the  quadratic  chordwlse  tem- 
perature distribution.  Figure  2.11  demonstrates  the  precipitous  degradation  of 
stiffness  for  the  fully  clamped  wing.  The  detailed  behavior  is  rather  different 
from  a temperature  distribution  computed  by  use  of  TRRID.  The  Instability  In 
that  case  was  connected  with  the  localized  leading  edge  buckling,  while  for  the 
parabolic  temperature  variation  four  thermal  bifurcation  modes  are  predicted 
for  leading  edge  temperatures  of  683.7,  687.0,  869.8,  and  901. 1°F.  All  the 
associated  modes  Involve  large  overall  deformation  as  opposed  to  the  localized 
buckling  In  the  previous  case. 

Results  obtained  for  a clipped  delta  wing  show  similar  trends  as  shown  In  Fig. 

2.12. 

In  order  to  verify  that  STAGS  and  TRRID  work  properly  in  combination,  other 
cases  were  executed  In  which  the  spatial  temperature  distribution  was  allowed 
to  vary  with  time  as  prescribed  by  the  results  from  TRRID.  it  was  observed 
that  In  that  case  the  solution  of  the  nonlinear  equations  (In  STAGS)  converged 
considerably  slower.  This  Is  tentatively  attributed  to  the  fact  that  presently 
used  extrapolated  initial  estimates  are  good  only  If  the  temperatures  vary 
proportionally  with  the  aerodynamic  pressure.  This  deficiency  should  be  rela- 
tively easy  to  remedy. 

It  should  be  emphasized  that  the  temperatures  for  leading  edge  buckling  and, 
to  a lesser  degree,  the  degradation  of  the  torsional  stiffness  are  determined 
in  a conservative  way  by  the  STAGS -TRRID  combination.  The  reason  for  this  is 
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that  the  relief  through  heat  transfer  in  the  chordwise  direction  has  been 
neglected.  A two-dimensional  thermal  analysis  would  certainly  result  in 
more  realistic,  less  conservative  results.  A more  realistic  treatment  of 
the  boundary  layer  transition  will  probably  work  in  the  same  direction. 

Finally,  the  results  obtained  for  wings  with  different  support  at  the  root- 
chord  indicate  that  deformation  of  the  mlssi''e  body  may  have  a substantial 
effect  on  the  stiffness  of  a heated  wing  and  consequently  on  its  aeroelastic 
performance.  (The  aerodynamic  coupling  between  missile  body  anti  fin  may 
have  an  even  larger  effect  on  the  aeroelastic  characteristics,  as  will  be 
discussed  in  the  following  section.) 
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Figure  2.2  Tip  Twist  as  a Function  of  Temperature 
Ratio  & Initial  Imperfection  Amplitude 
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Figure  2.5  Bending  S 
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Fig.  2.8  Stress  Resultants  in  Rectangular  Planform  Diamond  Profile 
Leading  Edge  Temperature  600°F 
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Figure  2.9  Rectangular  Wing  with  Variable  Grid 
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TORSIONAL  STIFFNESS  FRACTION 
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WING  L.E.  AND  T.E.  TEMPERATURE 


Section  3 

AERODYNAMIC  ANALYSIS 


It  was  indicated  in  the  Introduction,  Section  1,  that  numerous  theoretical  methods 
exist  throu^  which  the  hypersonic  steady  and  unsteady  aerodynamics  of  simple 
geometries  can  be  determined.  However,  the  methods  are  almost  without  exception 
too  idealized  to  be  of  much  help  in  the  design  of  hypersonic  vehicles.  This  is  because 
the  following  flow  phenomena  are  not  adequately  simulated: 

1.  The  invlscid  shear  flow  or  entropy  gradients  generated  by  the  small  amount 
of  nose  bluntness  that  is  required  for  thermodynamic  survival  in  hypersonic 
flow. 

2.  The  viscid-inviscid  flow  Interactions  which  become  important  at  hypersonic 
speed,  especially  in  the  lower  part  of  the  atmosphere  where  boundary  layer 
transition  greatly  complicates  the  physical  flow  simulation. 

The  main  efforts  of  the  aerodynamic  analysis  presented  here  are  concentrated  on 
removing  these  two  deficiencies  and  provide  analytic  means  that  can  produce  realistic 
predictions  of  the  hypersonic  steady  and  unsteady  aerodynamics.  In  the  analysis  to 
be  presented  only  small  perturbations  from  zero  angle  of  attack  will  be  considered, 
i.  e. , a = 0 if  not  stated  otherwise. 

3. 1 INVISCID  FLOW  OVER  AXISYMMETRIC  BODIES 

3. 1. 1 Bluntness 'Induced  Effects 

Nose  bluntness  effects  play  an  important  role  in  aerodynamics  of  bodies  at  high  super- 
sonic and  hypersonic  speeds.  A simple  analytic  theory  has  been  developed  earlier 
which  describes  the  unsteady  aerod3mamic  characteristics  of  slender  blunted  cones  at 
hypersonic  speeds  (Refs.  22  and  23).  In  the  present  analysis  the  theory  is  examined 
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in  regard  to  existing  limitations.  It  is  found  that  when  the  Mach  number  decreases 
much  below  = 20  or  the  cone  angle  becomes  less  than  the  Mach  angle,  the  unsteady 
embedded  Newtonian  theory  (Ref.  23)  is  not  applicable  in  its  original  form.  The 
needed  modifications  are  simple,  however,  and  a more  general  embedded  Newtonian 
theory  has  been  formulated. 


.3.1.2  Embedded  Newtonian  Flow 


The  simple  Newtonian  theory  (Ref.  41)  has  provided  an  indispensible  tool  for  hyper- 
sonic vehicle  design.  It  works  well  as  long  as  the  bow  shock  follows  the  body  contour 
relatively  closely,  such  as  is  the  case  for  pointed  bodies  and  for  very  blunt  geometries 
However,  when  the  bow  shock  cannot  follow  the  body  contour,  as  is  the  case  for  blunt 
nose  cylinder-flare  bodies,  for  example,  the  Newtonian  theory  gives  prohibitively 
large  errors.  Severe  loss  of  flare  efficiency  in  the  reduced  velocity  flow  field  behind 
a strong  bow  shock  has  been  noted  (Refs.  42  and  43).  Less  emphasis  has  been  placed 
on  the  effect  of  the  entropy  gradient  created  behind  a curved  bow  shock  (Ref.  44),  and 
the  tendency  has  been  to  "average  out"  the  velocity  variation  over  the  flare  (Ref.  45). 

In  the  present  analysis  the  complete  entropy  gradient  effects  are  considered. 


In  the  embedded  Newtonian  flow  concept  first  proposed  by  Seiff  (Ref.  45),  Newtonian 
theory  is  used  in  the  non-uniform,  inviscid  flow  field  defined  by  the  bow  shock 
(Refs.  44  and  46).  The  shock  shape  is  defined  by  the  nose  bluntness  and  can  be  re- 
lated to  the  nose  drag  (Refs.  47  and  48).  The  pressure  coefficient  on  the  "embedded" 
body  element  AA  in  Figure  3. 1 can  be  written: 


C = 
P 


(3.  1) 


where  is  well  approximated  by  the  stagnation  pressure  behind  a normal  shock 

(C  =1,  8)*‘for  blunt-bodies  and  by  the  Newtonian  value  (C  ^ = 2)  for  slender 

pmax  ' ' pmax  ' 


*For  M^  2 3 this  gives  less  than  3 percent  deviation  from  the  modified  Newtonian  value 
of  Ref.  49,  which  is 
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geometries.  Vj^  is  the  local  velocity  component  normal  to  the  surface  of  the  body 
element,  and  U is  the  local  axial  velocity  component  in  the  inviscid  shear  flow,  that 
is,  Vj^A/  is  determined  geometrically  and,  for  the  body  element  in  Figure  3. 1, 
pitching  with  the  rate  of  q radians /sec,  it  is 

(X  - Xpj,  + r tan  0 ) q 

Vj^/U  = cos  « sin  0 + sin  O'  cos  0 sin  ^ + cos  0 sin  <f>  (3.2) 


The  dynamic  pressure  ratio  is  a function  of  radial  position  in  the  inviscid  shear  flow 
profile.  If  similar  profiles  are  assumed,#  it  can  be  expressed  in  the  following  form 


= f*  (R/Rsh^ 


(3.3) 
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Rgjj  is  determined  by  the  nose  drag  (Refs.  47  and  50) 
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where,  for  M^  — oe,  a Xj^,  0.8<  K^j^ll.C.  For  hemisphere-cylinders,  = 1.0 
is  found  to  correlate  experimental  data  (Ref.  47). 


R is  obtained  from  Figure  3. 1 and  is  simply 


f*  (R/Rg|^)  In  Eq.  (3. 3)  Is  of  parabolic  shape  and  can  be  written 


(3.5) 


(3.6) 


^his  should  be  a good  assumption  except  for  the  flow  region  close  to  the  shock,  e.  g. , 
the  near  nose  region  for  small  a. 

'^#Deflection  at  x relative  to  the  nose,  x , due  to  body  deformation. 

o 
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Seiff's  data  for  a hemisphere-cjdinder  (Refs.  44  and  47)  show  that  the  dynamic  pres- 
sure data  can  be  fitted  by  the  following  polynomial 


a 

q 


P 


cc 


U“ 

oc 


f*  (X*)  = 0.  17  + 2.  7f)  X*  + 4X*“ 


(:i.  8) 


The  blast  wave  pressure  is  constant  over  the  inner  half  of  a cylindrical  shock  (Refs.  4 

and  50).  The  pressure  coefficient  C in  Eq.  (3. 1)  is,  therefore,  assumed  to  be 

po 

constant  over  the  radial  extent  of  a submerged  flare.  This  assumption  is  consistent 
with  the  assumption  of  similar  flow  profiles,  Eq.  (3.3).  Kuehn  (Ref.  51)  gives  the 
follow’ing  correlation  formula  for  the  axial  distribution  of  the  blast  wave  pressure 
p^  (0)  on  the  blast  wave  centerline  (also  measured  on  the  surface  of  a hemisphere- 
cylinder).^ 


Equations  (3.  l)-(3.  9)  give  the  following  embedded-Newtonian  formulation  for  the 
aerodynamic  pressure  coefficient  on  a surface  "embedded"  in  the  "entropy  wake" 
generated  by  a blunt  nose  at  hypersonic  speed  (M  = =c,  y = 1.4). 

CO 


c 

p 


+ C f*  (X*) 

PnEWT 


(3.  10) 


^The  pressure  given  by  second  order  blast  wave  theory  is  (Ref.  52). 


Po(»)/P„  = 


k(y)/A  + 0.44 


I 


For  0=0  and  small  p and  Q , the  pressure  on  an  embedded  flare  is  defined  by 
Eq.  (3. 10)  as 


C 

P 


C f*  (X*) 

PnEWT 


X 


N 


(3. 11) 


3. 1.3  Generalized  Embedded  Newtonian  Flow 


( 


It  has  been  shown  that  through  an  appropriate  modification  the  blast-wave  parameters 

can  be  extended  to  the  correlation  of  bluntness-induced  pressures  on  cones  and  wedges 

(Refs.  53-56).  For  cones,  where  P = $ , the  result  is  (Refs.  22  and  23) 

c 


Cp/2  = i*{X) 


(3. 12) 


The  theoretical  pressure  distributions  computed  by  the  M(etbod)  0(f)  C(haracterlstics) 
(Ref.  57)  at  Moo=  20  for  spherically  blunted  5®,  10®,  and  15®  cones  are  shown  in  this 
form  in  Figure  3, 2**.  Also  plotted  are  the  blast  wave  pressures  for  hemisphere- 
cylinders  using  Kuehn's  correlation  (Ref.  51).  The  difference  between  cone  pressures 
and  blast  wave  pressures  is  a measure  of  the  dynamic  pressure  profile  in  the  "entropy 
wake"  (Refs.  22  and  23). 

fW  = q/q  = (Ap  J /(Ap  J (3.13) 

“ ^ BLUNT/  SHARP 


BLUNT  / SHA 

where  tan  = Q , 
c c 


p 

C /2  6 for  slender  cones 
p c 
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The  f*-function  defined  in  this  manner  is  shown  in  Fig.  3.3  for  2.  r>®  < 0^<  20°, 

3 < M^<  20.  The  expected  collapse  of  the  pressure  distribution  curves  occurs  only 
for  cone  angles  above  a critical  value.  Closer  examination  reveals  that  the  critical 
value  is  the  freestream  Mach  angle.  Including  the  Mach  angle  in  the  shock-angle/ 
cone-angle  chart  of  Ref.  58  indicates  that  when  the  cone  angle  "dips  below"  p the 
shock  angle  remains  more  or  less  insensitive  to  (see  Figure  3.4).  In  tlfe  limit, 
when  0^  « \i,  the  shock  angle  is  completely  independent  of  fi  . Figure  3.3  shows 
that  the  entropy  profile  is  becoming  steeper  as  6^  is  decreased.  A steepening  of  the 
profile  is  also  observed  when  going  from  blunted  cones  to  blunted  cylinder-flare  bodies 
(Ref.  23).  At  hypersonic  speeds  the  f*-function  for  blunted  cylinder  flare  bodies  is 
(Ref.  21) 

f*  = 0.17+  2.  75  X*  + 4 X*~  (3.  8) 


and  for  blunted  cones  (Ref.  22)# 

X*  s 0.49 

0.49  <X*<  0.7  (3.14) 

0.7  < X* 


u.  i / 1-  1.  J /.I  A " 

f*  =1  1 - 3.43  (0.7  - X*l' 
1 


For  O'  = 0 the  similarity  parameter  X*  can  be  written 
X*  = |(x-x„)/d, 


'N 


tan^  o/c\C^ 

N 


(3. 15) 


Where  a = for  the  cylinder-flare  configuration  and  o = 6^  for  the  blunted  cone 
(see  Fig.  3.  5).  Comparing  Figures  3.  3 and  3.  5 one  easily  arrives  at  the  picture 
sketched  in  Figure  3.  0,  which  leads  to  the  following  "analytic  curve  fit"  for  blunted 
slender  cones  (see  also  Figure  3. 7). 


#Note  that  f = 2f* 
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X.  - Y* 


: x*S  X, 


: X„  < X*  < 


! X S y * 
X 


“•’  = -X„!  X„,t  ' (I  - f*)/AI! 

(tl-y/AD,  0.0  -2^(0. 0-f^, 

*0  ^ ’‘3  ! >‘3  ' >‘1  * 


B/2  (X*  - Xj) 


» 0.1G5  + [9.65/(M^  + 8.7))®:  5 3 


1.375  - 0.5  (0.17  (10  - MJ)®:  3 < < 10 


1,375 


3.33 


M > 10 

CD 


e^/li  £ 0.  23 


3.33  - 3.03  |((?*4i)  - 0.23|  : 0. 23  < 0*4:  < 1 


1 < e^/^ 


<Xo  X,)|l  [(X- - xJ)  AbJ  / txJ-xJ)ABj| 

<*  - 'i>  I*  * <x;  - xJ)  ab/  p - (xJ  - X?  abJ  J 


6 for  a « 0 
c 


For  0 /M  Si  and  M = «o , Eq.  (3.16)  reduces  to  Eq.  (3.14),  as  it  should. 
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A comparison  between  predictions  through  Eq.  (3. 16)  and  MOC-computations  is  made 
in  Figure  3.  8 for  Moo=  5 and  ^^0=  8.  In  addition  to  the  dynamic  pressure  profile  given 
by  Eq.  (3. 16)  the  velocity  profile  is  needed.  Through  reasonir  ’•  analogous  to  that  for 
hypersonic  f.ow  conditions  (see  Refs.  22  and  23),  the  velocity  ratio  is  defined  as  follows: 


g*{X*)  = UAJ„ 


1 = . X* 


k = ] 0.304  B/|l  + (3.33  - A)/2.33| 


(3.17) 


S*  = 


1 - k X 
g crit 


When  « # 0 the  similarity  parameter  X*  is  not  given  by  Eq.  (3. 15)  but  by  Eqs.  (3,4), 
(3.  5),  and  (3.  7)  using  the  tangent  cone  approximation. 


el  = sin*^  (Vl/U) 

V I /U  = cos  a sir  6 + sin  or  cos  i9  sin  d> 
-*•  c c 


(3.18) 


Incorporating  Eqs.  (3. 16)-(3. 18)  in  the  embedded  Newtonian  analysis  of  Reference  23 
gives  ihe  results  shown  in  Figure  3.9.  It  can  be  seen  that  the  stability  derivatives  of 
slender  blunted  cones  show  sensitivity  to  Mach  number  even  for  >10.  However, 
the  effect  is  probably  too  small  to  be  discernible  from  the  data  scatter  in  dynamic 
tests.  This  is  indicated  by  the  comparison  with  experimental  data  (Refs.  59-61) 
made  in  Figure  3. 10.  The  agreement  with  experimental  data  is  very  good,  as  good 
(at  least)  as  for  the  mu  jh  more  time  consuming  numerical  flow  field  theory  of  Rie 
et  al  (Ref.  62).  Note  titat  the  opposite  effect  of  nose  bluntness  on  static  and  dynamic 
stability  derivatives  is  well  predicted.  That  the  experimentally  observed  effects  of 
Mach  number  is  well  predicted  by  the  present  theory  is  documented  by  Figure  3. 11. 
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By  representing  the  nose  bluntness  effect  as  the  fractional  change  of  sharp  cone  char- 
acteristics the  following  advantages  are  gained:  Tangent  cone  methods  can  be  used 
even  at  supersonic  Mach  numbers,  and  predictions  made  by  inviscid  theory  can  be 
compared  with  viscous  experimental  data.  The  justification  for  this  statement  is 
simply  that  in  using  this  ratio  between  blunted  and  sharp  cone  characteristics  the  in- 
viscid and  viscous  crossflow  effects  that  are  common  to  sharp  and  blunted  cones  are 
eliminated.  The  same  is  true  about  sting  interference  effects  which  often  plague 
dynamic  wind  tunnel  tests  (Ref.  63).  Thus,  realistic  comparisons  between  inviscid 
predictions  and  experimental  results  can  be  made  as  long  as  nose  bluntness  has 
negligible  influence  on  cross  flow  effects  and  sting  Interference.  This  would  certainly 
be  the  case  for  tests  with  completely  turbulent  flow.  However,  hypersonic  testing  is 
usually  done  in  laminar  flow  with  a relatively  thick  boundary  layer,  and  nose  bluntness 
may  significantly  Influence  the  viscous  effects.  Orlik-Ruckemann  (Refs.  64  and  65) 
has  shown  that  the  "boundary-layer-cushionlng"  effect  usually  is  small  on  cones,  at 
least  for  standard  C.  G.  locations.  The  effect  of  nose  bluntness  on  viscous  ^rossflow 
with  associated  convective  time  lag  effects  should  also  be  small  as  long  as  boundary 
layer  transition  does  not  occur  on  the  (aft)  body*  (Refs.  24  and  25).  The  most  difficult 
problem  by  far  is  the  effect  of  nose  bluntness  on  dynamic  support  interference  (Ref.  63), 
especially  in  presence  of  aft  body  transition  effects  (Ref.  66).  There  is  experimental 
evidence  (Ref.  67)  indicating  that  nose  bluntness  could  have  a significant  effect  on  sting 
Interference.  Fortunately,  Walchner  and  Clay  made  sure  that  dynamic  support  inter- 
ference did  not  distort  their  experimental  results  (Refs.  59  and  68),  and  also  Ward's 
data  (Ref.  61)  can  be  expected  to  be  relatively  free  of  sting  interference  effects 
(Ref.  23).  Thus,  the  comparisons  made  in  Figure  3. 10  should  be  valid  indicating  that 
the  present  analytic  inviscid  flow  theory  can  adequately  predict  the  effect  of  nose 
bluntess  on  slender  cone  stability  characteristics.  The  comparison  made  in  Fig- 
ure 3. 11  is  also  valid  as  far  as  sting  interference  goes.  However,  boundary  layer 
transition  effects  have  distorted  the  low  Mach  number  comparisons  somewhat  in  the 
manner  described  in  Ref.  22.  The  M^,  = 5 data  have  been  corrected  for  asymmetric 
transition  effects  but  not  the  results  for  M^  = 3.  At  M^  = 10  no  asymmetric  transition 
effects  were  present 


♦Experiments  indicate  the  effect  to  be  less  than  :10  percent. 
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Comparing  Eqs,  (3.  8)  and  (3. 14)  in  view  of  the  trends  shown  in  Figure  3.5  leads  to  the 
following  definition  of  the  generalized  embedded  Newtonian  flow  field  over  a blunt-nose 
cylinder. 

f*  = f*+2BX*+  (1.45  B A *)^  < 1 (3.19) 

* 

where  f^  and  B are  obtained  from  Eq.  (3.  ir»)  ami  X*  is  obtained  from  Eqs.  (3.4), 

(3.  5)  and  (3.  7)  by  allowing  for  the  spanwise  coordinate  y/d^  (see  insert  in  Figure  3. 13) 
when  computing  R/dj^ , Eq.  (3.5). 

This  formulation  is  consistent  with  the  y~2  ratio  between  shock  radii  for  cones  and 
cylinders  used  in  Ref.  69.  That  is,  it  is  assumed  that  the  dynamic  pressure  profiles 
remain  similar.  The  dynamic  pressure  profile  above  the  cylinder  surface  .six  (6)  nose 
diameters  dowTistream  of  a hemispherical  nose  tip  is  showm  in  Figure  3.  ;2  fo»-  = 6. 
Also  shown  are  the  positions  of  the  conic  frustum  for  the  two  cone  angles  between  which 
the  entropy  wake  changes  from  cylindrical  to  fully  conical.  If  one  places  fins  of  unit 
aspect  ratio  with  their  leading  edge  at  (x  - x ) 'd.,  = r>.  thev  will  be  exposed  to  this 
dynamic  pressure  distribution,  as  is  illustrated  in  Figure  3. 13.  In  one  case  the 
missile  body  is  a simple  hemisphere-cylinder,  in  the  other  case  shown  the  same 
hemispherical  nosetip  is  followed  by  a shallow  conic  frustum,  0^  = 0.4  p,  and  a 
cylindrical  aft  body.  The  figure  illustrates  how  nonlinear  the  d>mamic  pressure  pro- 
files are,  and  how  a subtle  change  of  body  geometry  can  double  the  experienced 
dynamic  pressure  gradients. 

That  the  nonlinear  aerodynamics  are  real  is  illustrated  by  the  data  in  Figure  3.  14. 

The  present  theory  predicts  the  measured  nonlinear  (O')  characteristics  up  to 
O'  = 20'’.  At  higher  angles  of  attack,  bow  shock-flare  shock  interactions  cause  a loss 
of  flare  lift  (Ref.  70).  At  some  angle  of  attack  the  windward  fin  of  the  missile  configu- 
rations in  Figure  3. 1.3  will  start  to  interact  with  the  bow  shock.  At  = 6 this  will 
happen  when  a > 9°  for  the  sphere-cylinder  geometry  and  at  0^7°  for  the  sphere- 
cone-cylinder  geometry.  At  these  higher  angles  of  attack  it  is  the  shock-induced 
heating,  the  "blow-torch-effect"  described  by  Edney  (Ref.  71),  that  is  of  over-riding 
concern  rather  than  the  aerodjuamic  effects.  It  w'ould  be  highly  advisable  to  design 
the  missile  such  that  this  often  catastrophic  shock-induced  heating  is  avoided.  This 
can  be  accomplished  by  careful  design.  It  would  be  more  difficult  to  do  amThing  about 
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shock-induced  heating  on  forward  canard  surfaces.  The  only  recourse  designwise 
will  then  probably  be  to  sweep  the  leading  edge  back  until  the  shock-induced  heating 
predicted  by  Edney  (Ref.  71)  becomes  tolerable. 


3. 2 INVISCID  TWO-DIMENSIONAL  FLOW 
3. 2. 1 Sharp  Leading  Edge 

At  hypersonic  speeds  the  unsteady  aerodynamics  of  an  airfoil  with  sharp  leading  edge  is 
often  given  with  satisfactory  accuracy  the  simple  piston  theory.  When  the  piston 
velocity  is  so  small  that  no  shocks  are  generated  the  pressure  on  an  infinitely  thin  air- 
foil is  given  by  the  "simple  wave"  equation  (Ref.  72). 

2 y 

JL  = /i  + Izl  ^ - 1 (3.20) 

' 2 a«/ 

which  in  expanded  form  becomes 

js.  = 1 + y K + yjy/  1)  + yjy  1)  +.  . . (3.21) 

p«  ~ 12  \aj 

Piston  theory,  like  the  pure  Newtonian  theory,  considers  only  the  local  instantaneous 
flow  conditions.  Thus,  both  theories  are  of  quasi -steacly  nature  and  can  be  compared  in 
regard  to  accuracy  using  static  conditons.  With  U denoting  the  surface  inclination  to 
the  free  stream,  W/a^^  is 

W/a^  = sin  u cos  U (3.22) 
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For  the  small  disturbance  assumption  on  which  the  piston  theor>’  is  based,  the 
angle  v is  small  and 


W/a  = K = u 
00  00 

with  Cp  = 2 (p  - Pg^)/y  Eq.  (3.21)  together  with  Eq,  (3.22)  gives 


111  k'2 

12  ^ 


(3.  24) 


2 

As  Cp/u  =2  according  to  pure  Newtonian  theory,  Eq.  (3.  24)  indicates  that 
piston  theory  and  Newtonian  theory  should  give  very  different  results  depending 
upon  the  Mach  number  (and  flow  deflection  angle).  Figure  3. 15  shows  that  the 
differences  indeed  are  large.  This  is  puzzling  in  view  of  the  relatively  close 
agreement  between  the  two  theories  found  by  Yates  and  Bennet  (Ref.  76)  in  their 
hypersonic  flutter  calculations.  The  reason  for  this  will  soon  become  evident. 

The  static  pressure  ratio  through  a "tangent  wedge"  compression  or  expansion 

## 

is  (See  Refs.  41  and  77). 


#The  thickness  effects  which  are  present  in  supersonic  flow  (Ref.  73)  disappear 
already  at  > 2 (Refs.  74  and  75). 

##The  expression  for  K > 0 is  equal  to  the  so-called  strong  shock  piston  theorj’ 
when  K > 1 (Ref.  78). 
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The  pressure  ratios  given  by  piston  theory  and  Newtonian  theory  are 


(-E.) 

' P® ' pT 

= 1 + YK  1 

^ Y + 1 

+ 4 K + K 

0 

; K<  0 

= 1 +< 

\Pcx)/  N 

YK^ 

; K £ 0 

The  pressure  ratios  given  by  Eqs.  (3.25)  through  (3.27)  are  shown  in  Figure  3. 16  for 
air  ( V = 1.4).  One  can  see  that  there  is  a substantial  K- range  for  which  the  slopes 
d(p/p^)/dK  given  by  Eqs.  (3.25)  - (3.27)  will  not  differ  much,  indicating  that 
stability  derivatives  and  flutter  boundaries  would  not  be  very  different  either,  in  agree 
ment  with  the  results  obtained  by  Yates  and  Bennet  (Ref.  76).  Figure  3. 16  also  shows 
that  Newtonian  and  tangent  wedge  pressure  ratios  ^roach  different  limits  for  large 
Mach  numbers  (M^— m;.  This  can  be  seen  directly  from  Eq.  3.25  which  gives 


lim 
K— 00 


1 


9 

This  agrees  with  the  Newtonian  value  1 + YK  , given  by  Eq.  (3.27),  only  for  Y = 
In  order  for  the  Newtonian  theory  to  give  the  ’’correct"  limit  also  for  other  Y - values 
the  following  "gamma-correction"  has  to  be  made. 


= 1 . nr*  1) 


#This  corresponds  to  the  value  C 


Y + 1 used  in  Ref.  76  . 


(3.26) 


(3.27) 


(3.28) 

1. 


(3.29)# 
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For  stability  and  flutter  considerations  the  pressure  derivatives  are  more  interesting. 
They  are 


V /Pi 


; - 1 < K < 1 


(3.31; 


K < 0 
K 2 0 


(3.32) 


Figure  3. 17  shows  the  derivatives  givoi  by  Eqs.  (3.30)  through  (3.32)  for  V = 1.4.  It 

i^ipears  that  piston  theory  can  be  used  for  - 1 < K < 1 and  "gamma-corrected" 

Newtonian  theory  for  K > 2.  The  h}rpersonic  similarity  results  above  can  be  extended 

/ 2 

down  to  (low)  supersonic  Mach  numbers  by  substituting  the  Prandtl-factor  v -1 
for  Mg  in  the  definition  of  K in  Eq.  (3.23)  (Refs.  41  and  79).  Eq.  (3.31)  shows  that  for 
y = 1.4  the  quadratic  term  (cubic  in  Eq.  (3.26))  increases  the  linear  slope  (obtained 
at  K = 0)  with  60  percent  at  the  limit  (K  = ±1)  for  the  piston  theory.  Thus,  as 
Lighthill  pointed  out  (Ref.  72),  it  is  important  to  include  this  last  term  of  the  piston 
theory,  and  the  so  called  third  order  piston  theoiy  should  be  used.  However,  if  the 
piston  theory  is  used  only  to  describe  the  unsteady  perturbaticm  from  the  steady  local 
flow  condition,  one  needs  only  the  linear  term,  as  long  as  the  unsteady  perturbation  is 
small,  e.g. , for(6c/U;j,)  «1  (see  Ref.  76). 
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For  a ilgid  airfoil  describing  oscillations  of  low  reduced  frequency,  ( ccu/Uoo  )^«  1, 
the  hypersonic  similarity  parameter  can  be  written 

K = M„Vj./U 


where  from  Figure  3. 18  a V|_  / U is  defined  as  follows  for  lower  and  upper  surfaces. 


I^l\  ■ *CG  + yL  ^ 

= sin(OL  ^ — cos  u. 


(H 


= sin  (Ou  - Oq  - 0)  - 


U 


(X  - + yjjtanuy)q 

U 


cos  V 


u 


V3  = tan”^  (dy/Sx) 


6 is  the  perturbation  in  pitch  aroimd  the  trim  angle  of  attack  . For 
^o  = ^ = Q = 0 this  definition  of  K agrees  with  that  given  earlier  in  Eq.  (3.23). 
The  surface  pressure  coefficient  is  determined  as  follows: 


(3.33) 


(3.34) 


(3.35) 
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In  the  case  of  the  elastically  deforming  airfoil,  Eq  (3. 34)  becomes 


sin  +|A)  + -i  eosUj_ 

sin  (U^.  - “o  ■ 'll’  ■ U If  “l- 


(3.36) 


h is  the  centerline  displacement  normal  to  the  rigid  wing  chord  (see  Figure  3. 18b).  In 
the  aeroelastic  analysis  to  be  discussed  later  the  deflection  h is  the  sum  of  the  defor- 
mations in  a large  number  of  elastic  degrees  of  freedom.  A special  case  which  has 
been  treated  extensively  is  the  so  called  binary  flutter  in  which  the  \ving  is  deforming 
in  bending  and  torsion  only;  i.  e. , the  airfoil  section  is  undeformed  (see  Figure  3. 18c). 
For  this  case  £q.  (3. 36)  becomes 
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h^  is  the  bending  deformation  of  the  airfoil  elastic  axis  (see  Figure  3. 18c). 


The  chordwise  load  distribution  is  in  coefficient  form 


(3.38) 
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At  - 0 Eqs.  (3.33)  - (3.38)  define  the  following  derivatives  for  a symmetric, 
sharp  wedge  (i)j^  = ^ ^w^ 


3cn/3{  = 2 M„Cp^  “»  e„ 


(3.39) 


^CNq/Sl  - 2M«Cp^  [4(1  + tan^0^)  - 4^]  cos  0^ 


(3.40) 


where  , obtained  from  £q.  (3.30),  is 
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For  a thin  wedge  (tan  0^  « 1 and  cos  0^  * 1 ) the  following  stability  derivatives  are 


obtained  at  a =0 
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where  Cpj^  is  given  by  Eq.  (3.41)  with  K = (or  K = '^^oo  ~ ^ 

low  supersonic  Mach  numbers). 

For  a double-wedge,  diamcmd  airfoil  shape,  the  result  at  = 0 is 

^‘^ne^F  " 1 

^‘^ne^A  ^ Sk-  / (3.43) 
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where  C „ 
pK- 

is  obtained  from  Eq.  (3.41) 

with  K 

= (or 

- - 1 ) 

Equatiixis  (3.41)  - (3.43)  show  that  M^c^^  and  are  functions 

of  K only.  This  is  illustrated  in  Figure  3. 19  for  single  and  double  wedges  Fig- 
ures 3. 19a  and  3. 19b,  respectively).  As  e:q)ected  the  tangent  wedge  values  are  well 
represented  by  third  order  piston  theory  for  low  K - values,  K < 1,  and  by  the 
gamma-corrected  Newtonian  theory  for  high  K - values,  K > 1.  Orlik-Riickemann's 
experimental  pitch  damping  data  for  a double  w'edge  (Ref.  80)  fall  well  above  the  inviscid 
prediction  (see  Figure  3.19b).  The  Helium  data  point  (y  = 5/3)  has  been  corrected 
for  y-effects  using  Eq.  (3.43).  One  feels  inclined  to  agree  with  the  authors  that  the 
observed  50  percait  overshoot  of  inviscid  predictions  cannot  be  explained  by  pure 
boundary  layer  effects  (Refs.  64  and  65),  but  that  some  other  effects,  e.g. , shock  wave- 
boundary layer  interaction  at  the  trailing  edge,  must  have  caused  it.  One  such  mecha- 
nism has  becai  demonstrated  by  Hulcher  and  Behrens  (Ref.  81).  They  showed  that  at 
Ot  = 15°  the  leeside  laminar  flow  on  a flat  plate  at  = 6 was  separated 
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from  70  percent  chord  to  the  trailing  edge  due  to  the  wake  recon^)ressioa  (shock).  Thus, 
one  can  also  expect  the  9°  double  wedge  to  have  a trailing  edge  separation 
(at  a = 0),  especially  in  view  of  the  higher  Mach  number.  It  has,  in  fact,  been  shown 
that  a biconvex  airfoil  experiences  laminar  separated  flow  on  the  leeside  aft  half  chord 
(Ref.  82).  The  loss  of  leeside  expansion  due  to  separation  will  generate  a negative, 
statically  destabilizing  tail  load.  The  dynamic  effect  will  be  stabilizing  due  to  the  con- 
vective time  lag  effects  associated  with  this  wake- induced  separation  (Refs.  83-85),  and 
could  easily  cause  the  observed  50  percent  increased  damping.  The  method  presented 
in  Reference  82  for  prediction  of  the  separation  induced  pressure  increase  at  supersonic 
Mach  numbers  can  possibly  be  implied  also  in  the  region  of  present  interest,  3 < < 6. 

Thus,  with  the  information  contained  in  References  82  through  85  as  a basis,  analytic 
methods  could  be  developed  for  prediction  of  the  separated  flow  effect  on  the  unsteady 
aerodynamic  characteristics  of  double-wedge  or  lenticular  shape  airfoils.  This  is, 
however,  beyond  the  scope  of  the  present  study. 

For  a single  wedge  no  trailing  edge  separation  is  present  and  the  deviations  between 
inviscid  predictions  and  preliminary  experimental  results  (Ref.  86)  seem  to  be  accounted 
for  by  the  attached  flow  boundary  layer  effect  described  Orlik-Ruckemann  (see  Fig. 
3.20).  However,  one  should  be  careful  not  to  interpret  the  results  as  proof  positive 
because  of  the  inherent  difficulties  in  performiirg  a two-dimensional  (itynamic  test,  as 
the  authors  of  Reference  86  point  out.  The  peculiar  frequency  effects  observed  in  the  test 
could  well  have  been  caused  by  wind  tunnel  interference  as  there  is  no  known  reason 
for  the  effects  to  be  present  in  the  free  flight  case.  Possibly  the  only  way  to  obtain  the 
desired  strip-data  (two-dimensional  airfoil  data)  would  be  to  test  a half  model  of  the 
missile-fin  configuration  varying  the  aspect  ratio  of  the  fin  (Ref.  87).  In  such  a test  the 
effect  of  missile  nose  bluntness  on  the  unsteady  fin  aerodynamics  could  also  be  in- 
vestigated to  check  out  the  prediction  techniques  described  in  Section  3. 1. 
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3. 2. 2 Bluntness  Induced  Effects 


The  blast  wave  pressure  in  two-dimensicxial  flow,  as  it  would  be  measured  on  a flat 
plate  with  blunted  leading  edge  (the  2-0  equivalent  to  Eq.  (3. 9)  ),is  as  follows 
(Refs.  50-53).# 
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+ 0.  56 
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■ “’n) 
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k'(y)  = 0.121  for  r =1.4 


(3.44) 


Thus,  the  following  two-dimensional  embedded-Newtonian  formulation  is  obtained  for 
the  aerodynamic  pressure  coefficient  on  a surface  ''embedded"  in  the  "entrc^y  wake" 
at  hypersonic  speeds  (Mco-*  ® ) 
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(3.45) 


jfSubscript  2D  or  a prime  siqjerscript  is  used  to  indicate  two-dimensional  qu  ntities 
eqidvalent  to  those  used  earlier  in  the  three-dimensional  case. 

##The  pressure  given  by  second  order  blast  wave  theory  is  (Ref.  52) 

Po<°^/Pa,  = k'(y)/A'  - 0.56 

That  is,  the  same  85  percent  "efficienty  factor"  measured  by  Kuehn  for  hemisphere- 
cylinders  (Ref.  51)  is  applied  here.  No  attempt  is  made  to  modify  the  constant  0. 56 
as  it  is  unimportant  for  the  determination  of  aerodynamic  force  derivatives. 
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For  0=0  the  similarity  parameter 
slender,  blunted  wedge  is  (see  Chemyi, 


X that  correlates  the  pressures  go  a 
Ref.  53). 


'2D 


= tan^  6 c,~  ^ 

w dN  ^ 
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(3.47) 


♦ 

That  is,  the  parameter  X 2q  at  O ^ 0 
by  Eq.  (3. 7)  for  the  3-D  case,  is 


corresponding  to  the  parameter 


X*  defined 


(3.48) 


Here  d'^  is  the  thickness  of  the  blunted  leading  edge  and  (R^^^  ^^^2D 
follows  (Ref.  52) 


'# 

It  is  assumed,  now  as  before,  that  and  = 1 . For  a slender  wedge 

with  (semi)  cylindrical  leading  edge  the  results  shown  in  Figure  3.21  are  obtaimxl, 
which  give  ^ 2D  That  is,  the  pressure  on  a wedge  with  semi-cylindrical  leading 


Using  another  value  would  give  a different  definition  of  f 2D*  ^ would,  however,  give 
the  same  Cp-value,  Eq.  (3. 45),  because  f £D  ^ defined  as  the  difference  between  the 
surface  pr^sure,  computed  M O C (Ref.  57)  and  the  blast  wave  pressure,  given 
by  Eg.  (3.44). 


edge  bluntness  is  simply  the  sum  of  the  blast  wave  pressure  and  the  pressure  measured 
on  a sharp  wedge. 
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(Cp)  gjj^p  is  given  by  Eq.  (3.35) 


(3.50) 


That  there  is  no  velocity  deficit  effect  on  the  wedge  pressures,  only  on  the  cone  pres- 
sures, is  explained  by  the  differences  in  flow  fields.  For  the  wedge  all  flow  turning 

occurs  at  the  shock,  whereas  substantial  flow  turning  occurs  downstream  of  the  shock 

2 

on  the  c(me.  For  a slender  blunted  wedge  (tan  0^ « 1 and  cos  0^  = 1 ) the  deriva- 
tives correspxxiding  to  the  sharp  wedge  derivatives  given  in  Eq.  (3.42)  are 
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For  a (semi)  cylindrical  leading  edge 
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In  Figure  3. 22  these  inviscid  predictions,  Eqs.  (3.  51)  and  (3. 52),  are  compared  with 
East's  experimental  data  (Ref.  78).  The  static  data  seem  to  agree  with  predicti«is, 
considering  the  data  scatter,  whereas  the  dynamic  experimental  data  do  not.  For 
completeness,  sharp  wedge  data  are  also  shown.  East  shows  how  the  wedge  shock 
caused  substantial  separation  of  the  sidewall  boundary  layer,  as  would  be  expected  (see 
Figure  3. 23  and  Ref.  88).  It  can  be  seen  that  for  = 0 the  side  plate  boundary 

ULr 

layer  sqiarates  well  upstream  of  the  wedge,  and  as  much  as  35  percent  of  the  wedge 
surface  lies  within  the  influence  of  the  sidewaU  shock  according  to  East.  For 
^CG  ~ 1*0  ^6  corresponding  value  is  15  percent.  Thus,  the  wall  interference  varies 
with  It  will  also  be  different  for  the  sharp  and  blunted  wedge.  Already  at  low 

speeof  these  sidewall  effects  cause  dynamic  simulation  problems  (Kei.  89),  problems 
which  will  become  even  greater  at  hypersaiic  speeds  with  its  complicated  "comer-flow” 
(see  Figure  3.24  and  Refs.  27,  90,  91,  and  92).  Thus,  the  deviations  between  present 
predictions  and  East's  experimental  data  could  be  caused  ty  this  side  plate  interference. 
East  used  experimental  wedge  data  to  obtain  a dependence  of  the  blast  wave  pressure, 
in  Eq.  (3. 50),  on  the  "tangent  wedge"  angle.  As  the  inclusion  of  this  wedge  slope 
dependence  of  did  not  give  any  more  signiflcant  effect  of  nose  bluntness  on  the 
predicted  stability  derivatives  than  the  present  simpler  formulation,  no  refinements  of 
Eq.  (3. 50)  are  contemplated  at  this  time.  It  should  probably  be  pointed  out  that  the 
viscous  interacticm  effects,  which  will  be  discussed  later,  are  negligible  in  East's  test, 
mainly  because  of  the  relatively  cool  wall,  T /T.  = 0.33. 

W I 


Although  East's  data  are  somewhat  difficult  to  use  when  checking  two-dimensional 
theoretical  estimates  because  of  the  large  interference  effects  from  the  side  plate 
boimdaiy  layer,  they  are  very  valuable  in  that  they  point  out  what  the  effect  from  the 
missile  body  boundary  layer  could  be  on  the  fln  stability  derivatives.  If  one  wants  to 
use  all-movable  fins  with  some  margin  of  static  stability  one  could  end  iq}  with  an  axis 
position  in  the  range  0.2  < 4 < 0. 4;  i.  e. , the  very  range  in  which  East  measured 

negative  aerodynamic  damping.  Viscous  crossflow  affects  on  the  missile  body  could 
possibly  aggravate  the  situation.  Thus,  one  will  have  to  stucty  this  side  wall  effect 
further  if  one  wants  to  design  a finned,  hypersonic  missile.  If  one,  instead  of  movable 
fins,  decides  to  use  flap  controls,  one  will  encounter  other  problems.  Although  the 
velocity  deficit  effect  is  absent  for  a blunted  wedge,  f = 1 in  Eq.  (3,45),  a deflected 
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flap  downstream  of  a blunted  leading  edge  will  experience  entropy -wake-induced 
ncm -uniform  flow  effects  as  illustrated  by  the  curved  flap  shock  on  a lifting  reentrj' 
vehicle  (Refs.  93,  94,  and  Figure  3.23).  At  high  angles  of  attack  bow  shock  - flap  shock 
interaction  presents  additional  problems  (Ref.  70). 


3.3  VISCOUS  FLOW  EFFECTS 

The  boundarj’  layer  on  the  fins  of  a tactical  missile  operating  in  the  sf>eed  range 
3 ^ vai7  from  being  completely  laminar  at  the  high  Mach  number  - high 

altitude  end  of  the  operating  range  to  being  main!’'  turbulent  at  low  Mach  numbers  iuid 
altitudes.  That  is,  at  some  intermediate  region  of  the  operating  range  boundary  layer 
transition  will  have  a dominating  influence  on  the  aerothermodj’namics  of  the  fins.  In 
regard  to  the  boundary  layer  on  the  missile  body  transition  will  play  a role  throughout 
the  whole  operating  riuige,  and  will  significantly  affect  the  unsteady  aerodjTiamics 
(Refs.  24  and  25),  effects  that  may  become  adverse  if  the  missile  tip  has  even  a small 
amount  of  nose  bluntness  (Refs.  95  and  96). 

As  the  characteristic  time  of  the  thermal  response  usually  will  be  orders  of  magnitude 
larger  than  the  periods  of  structural  vibrations  one  can  negiect  any  coupling  between 
thermodynamic  and  elastic  response.  When  considering  the  rigid  body  response 
thermodynamic  coupling  effects  may  not  be  negligible.  However,  this  problem  will  not 
be  addressed  in  the  present  study  which  only  is  concerned  with  the  aero-thermo-elastic 
response  oi  the  missile  fins.  Thus,  only  the  viscous  effects  on  the  unsteady  aero- 
dynamics of  the  fins  will  be  considered  here. 

3.3.1  Two-Dimensional  Viscous  Flow 

Viscous  flow  effects  in  low  density  hypersonic  flow  have  been  widely  investigated  by 
both  theoreticians  and  experimentalists.  This  is  at  least  true  in  regard  to  static 
aerodjTiamics  w'here  the  so-called  viscid-inviscid  interaction  effects  are  well  established 
(see  Reference  97  for  a summary^  of  the  present  state  of  the  art).  The  dynamic  effects 
of  this  hypers  nic  viscid-inviscid  interaction  have  received  much  less  attention,  the 


3-24 


only  well  known  work  being  that  by  Orlik-Riickeman  et  al  (Refs.  64,  65,  80,  86,  and 
98).  In  the  present  study  a simplified  approach  is  used,  tailored  to  the  particular  flow 
conditions  of  interest  for  a tactical  missile  operating  in  the  low  hypersonic  flow  regime. 

3. 3. 1. 1 Static  Aerodynamics.  The  boundary  layer  displacement  surface  will  increase 
the  local  inviscid  wedge  angle  6 by  an  amount  d5  /dx.  Figure  3. 19  shows  how 
this  will  increase  the  normal  force  derivative. 

The  displacement  thickness  for  a laminar  boundary  layer  is  (Ref.  41) 
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temperatures  in  the  range  200°R  < T*  < 1000“R, 
where  T*/T^  = (1  + 3 T j 6 (Ref.  97) 


(3. 53a) 


(3.53b) 


Alternatively  Eq.  (3.  53)  can  be  written  as  follows: 


(3.54a) 


(3.54b) 


For  the  isolated  wall  T = yp^ 
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The  thickness  slope  is  simply 


/dx  = 6 /2x 


(3.55) 


The  slope  ^6  /dx  will  be  very  small  compared  to  the  angular  changes  of  the 
inviscid  flow,  at  least  in  the  present  study,  and  the  effect  of  o6  /dx  can  be 
simulated  by  a small  wedge  embedded  in  the  flow  aft  of  the  bow  shock.  Thus, 
Eq.  (3. 26)  can  be  used  to  give  for  » 1 

AP  /p  = = TM  d6*/dx 

^ *^v  *^e  V e 


Combining  Eqs.  (3.53)  through  (3.56)  gives 


Ap/Pg  = 

K = M deVdx  = BX 
v e 

B = Bj  + Bj/M^ 

yc'/v/R^ 

where  2B  can  be  identified  with  the  bracket  in  the  equations  describing  6*.  Eqs.  (3. 53) 
and  (3. 54). 

For  » 1 Eq.  (3.57)  simplifies  to 
Ap/Pe  = VB^X 


(3.57) 


(3.58) 


Predictions  of  p /p  = 1 + Ap  /p  using  Eq.  (3.  58)  agree  well  \vith  experimental 
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wedge  data  for  = 10  (Ref.  97  and  Figure  3. 26  ).The  agreement  persists  up  to  large 
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X -values  (outside  the  range  of  interest  for  the  present  study).  It  should  be  pointed 
out,  however,  that  when  the  wedge  angle  goes  to  zero,  and  all  the  flow  change  is  of 
viscous  nature,  Eq.  (3. 58)  no  longer  predicts  the  measured  pressures.  In  this  case 
viscid-inviscid  coiq)ling  effects  have  to  be  included  (sc  e Refs.  41  and  97).  However, 
on  the  finite  chord  airfoil  of  consideration  in  the  present  study,  trailing  edge  and  sepa- 
rated flow  effects  will  often  dominate  when  these  higher  order  viscous  effects  should 

be  considered;  e.  g. , ata  0 for  finite  wedge  angles  (0  > 0).  Eq.  (3.57),  or 

2 ^ 

Eq.  (3. 58)  for  » 1,  is  especially  conveniait  v^4ien  using  the  local  linearization 

^qjproach  in  flutter  analysis  as  was  done  by  Yates  and  Rennet  (Ref.  76).  As  the  self- 

induced  pressure  does  not  change  the  boundary  layer  slope  d 6 /dx  up  to  second  order 

accuracy  (Ref.  41),  the  viscous  induced  effect  on  the  surface  pressure  derivative  is 

solely  a function  of  inviscid  flow  parameters.  The  following  simple  result  is  obtained. 
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Pg/P-e  is  obtained  from  Eq.  (3.25).  The  derivative  ^ [log  (ApyPg)|  /s(Pg/p^)  is 
obtained  from  Eq.  (3. 57)  using  the  oblique  shock  relations  given  in  Ref.  58.  The 
result  is  as  follows: 


3-28 


Combining  (Eqs.  (3.59i  and  (3.60)  gives 


(3.61) 


and  B = are  defined  in  Eq.  (3.57),  and  p^/p^  is  given  by 

Eq.  (3.25),  substituting  K with  K / 1 - at  lower  supersonic  Mach  numbers 

(Ref.  79).  The  relative  change  of  inviscid  aerodynamic  derivatives  due  to  viscous  first 

order  perturbation  effects  can  be  obtained  for  arbitrary  thin  airfoils  using  Eq.  (3.G1). 

For  a wedge  the  result  becomes  particularly  simple,  as  M = constant  and 
— — — 1 /2  ® 

X = Xte  ^ • Integration  of  Eq.  (3.61)  gives  then  for  upper  and  lower  surfaces 

the  following  relationship. 
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Noticing  that  (C  = (c  ^ , Eq.  (3.62)  gives  the  viscous  induced  increment  of 

p 0 e nt7  1 

upper  and  lower  surface  inviscid  aerodynamic  derivatives.  At  Q = 0 it  gives  the  total 
result  for  both  surfaces.  That  is 
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Within  the  linear  Q -range  Eq.  (3. 63)  also  expresses  the  ratio ac^^/c^..  In  Fig- 
ures 3.27  and  3.  28  predictions  through  Eq.  (3.  63)  are  compared  with  available  experi- 
mental data  (Refs.  99  luid  100).  The  good  agreement  for  the  5°  wedge 
(Fig.  3. 27)  was  expected,  but  the  good  agreement  with  the  flat  plate  data  (Fig.  3. 28) 
was  a surprise  in  view  of  the  failure  of  Eq.  (3. 58)  to  predict  Sliski’s  flat  plate  data 
(Ref.  97).  The  experimental  data  were  underpredicted  by  30  percent.  This  may  also 
be  the  lift  loss  on  a finite  chord  flat  plate  due  to  the  trailing  edge  effect  discussed  by 
the  authors  of  Ref.  100.  (See  inset  in  Figure  3.28).  Thus,  fortuitously  Eq.  (3,63)  will 
give  the  cori’ect  viscous  effect  on  the  lift  even  for  the  limiting  case  when  0 ^0. 

However,  the  moment  ccxitribution  will  not  be  predicted  as  well,  and  neither  will  the 
viscous  effect  on  unsteady  aerodynamics.  However,  as  was  discussed  before,  a very 
accurate  prediction  of  the  pure  viscous  interaction  for  zero  inviscid  flow  inclination  to 
the  surface  is  not  essential  for  prediction  of  the  airfoil  aerodynamics.  Judging  by  the 
good  agreement  between  predicted  and  measured  viscous  influence  on  the  moment  of  the 
5°  wedge,  it  does  not  seem  to  experience  much  of  this  trailing  edge  effect  (see 
Figure  3. 27).  This  is  probably  to  a large  extent  due  to  the  modest  value  (X-j.^  = 1. 3) 
of  the  hypersonic  viscous  similarity  parameter.  As  this  parameter  X.j,j,  is  increased, 
the  strong  interaction  region  near  the  leading  edge  and  the  near  wake  effect  at  the  trailing 

edge  both  become  more  important.  For  a finite  base  height,  e.  g. , for  0 0,  the 

w 

wake  recompression  effects  are  fed  upstream  through  the  near  wake  adding  complexitj’ 
to  the  trailing  edge  effect  (Refs.  84  and  85).  For  the  X.j.j,  -range  of  interest  in  the 
present  study  these  effects  are,  however,  unlikely  to  become  important.  The 
viscous  effects  are  likely  to  be  small,  as  is  illustrated  by  Figure  3. 29,  at 
least  for  steady  state  flight  conditions,  T^^/T^  = 1.0.  When  the  tactical 
missile  is  maneuvering,  it  may  be  possible  for  the  fins  to  reach  soak  tem- 
peratures as  high  as  those  expected  on  the  straight  wing  space  shuttle  during 
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reentry;  i.e.,  =2.5  (Ref.  101).  It  is  interesting  to  note  that  for  = 6 

this  condition  comes  close  to  the  adiabatic  wall  case  existing  in  most  wind  tunnel 
tests.  Even  in  this  case,  however,  the  viscous  induced  increment  of  the  inviscid 
lift  will  not  exceed  20  percent.  Consequently,  inviscid  aerocfynamic  characteristics 
have  been  used  in  the  present  aeroelastic  analysis. 

3. 3. 1.2  Unsteady  Aerodynamics.  The  change  of  effective  airfoil  geometry  produced 
by  the  boimdary  layer  displacement  surface  and  the  corresponding  effect  on  static 
stability  derivatives  were  considered  in  the  previous  section.  There  is,  of  course,  a 
similar  effect  on  the  dynamic  stability  derivative.  Equation  (3.61)  defines  the  following 
strip  load  derivative  distribution . 


di 

dd  J 

L 

de 

(3.64) 


And  the  running  load  increment  due  to  viscous  wedge  thickening  is 

''l 
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_l_nW 
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(3.65) 


For  a wedge  describing  rigid  body  oscillations  in  pitch  around  x = Eqs.(3.34), 


(3.61),  (3.  64)  and  (3. 65)  give  for  = 0 
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TE 
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Integrating  Eq.  (3. 66)  and  its  moment  contribution  gives 

“ (3  “ ^cg) 

A + (i  . t 
45  \3  ^OCl 


(3.67) 


(3.68) 
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The  ratio  given  by  Eq.  (3.61),  as  was  discussed  earlier  (See 

Figures  3, 27  - 3. 29).  Equation  (3. 68)  shows  that  the  viscous  induced  thickening  of  the 
wedge  increases  the  aerodynamic  damping,  the  effect  being  minimum  for  £„_  = 1/3. 
Orlik-Rlickeman  has  shown  that  in  addition  to  this  static  interaction,  it  is  important  to 
consider  the  "dynamic  interaction"  resulting  from  the  fact  that  the  boimdar>’  layer  dis- 
placement thickness  does  not  remain  constant  during  the  oscillation  (Ref.  98).  The 
boundary  layer  becomes  thinner  on  the  windward  side.  On  the  oscillating  airfoil  this 
boundary  layer  "cushioning"  decreases  the  "piston  effectiveness"  by  changing  the 
inviscid  velocity  ratio  (Vj^  ^^e^i  ^ amount  (d  6 /dt)/U^.  This  local  reaction  can 

be  assumed  to  occur  instantaneously  for  the  structural  and  rigid  body  frequencies  of 
practical  interest.  At  large  V one  may  have  to  consider  downstream  effects  of 
changes  occurring  at  an  earlier  time  instant  in  the  strong  interaction  region  near  the 
leading  edge,  as  well  as  upstream  effects  of  changes  in  the  near  wake  recirculatory 

flow  (also  occurring  at  an  earlier  time  instant).  These  effects  are  not  considered  in 
the  present  analysis. 

For  most  applications  it  can  be  assumed  that  the  missile  velocity  changes  slowly 

* 

compared  to  structural  and  rigid  body  frequencies.  Thus,  the  time  derivative  d6  /dt 
is  solely  the  result  of  the  perturbations  of  the  relative  piston  velocity  V 
defined  by  Eqs.  (3.34),  (3. 36),  or  (3. 37),  That  is 


d6*  _d  ) 

dt  dt  \u^/ 


(3.69) 


Whereas  the  boundary  layer  slope,  d5  /dx  , is  relatively  insensitive  to  the  ambient 
pressure  gradient,  the  boundary  layer  thickness  itself  is  not,  but  becomes  thinner  for 
increasing  pressures  (Ref.  41).  Bertram  et  al  (Ref.  102)  have  shown  that  this  fact 
can  be  represented  by  assuming  the  boundary  layer  thickness  to  be  inversely  proportional 
to  \/"p  /p  . Consistent  with  the  perturbation  approach  used  here  p /p  is 

V OO  V CO 

represented  by  the  inviscid  pressure  ratio  Pg/p^  • Thus,  Eq.  (3.  55)  is  modified  as 
follows: 


(3.  55u) 


When  p /p  — 1 the  flat  plate  value,  Eq.  (3.  55),  is  recovered. 

C CO 
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From  Eqs.  (3.55a),  (3.57),  and  (3.60)  is  obtained 
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The  change  of  piston  effectiveness  ^ Vu^  can  be  written 
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Again,  using  the  embedded  formulation,  Eq.  (3. 56),  gives 

AP^(6*)/Pg  = ‘>'Mg6*/Ue 


Combining  Eqs.  (3.70,  (3.71),  and  (3. 72)  gives  finally 
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In  general  M and  AP  /p  are  functions  of  However,  for  a wedge 
Eq.  (3. 73)  simplifies  to 


TE  ^1/2 
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Integrating  Eq.  (3.74)  gives  the  following  result  for  pitch  oscillations  around 

\ 


X = X cG  a = 0 . 
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(3. 75) 
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Combining  Eqs.  (3.61),  (3.62),  (3.63),  (3.68),  and  (3. 75)  with  Eq.  (3.42)  gives 
the  damping  ratios 
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and  the  following  viscous-induced  forward  shift  of  the  aerodynamic  center, 


(3.  78) 


Gg  varies  according  to  the  wall  temperature  ratio.  Figure  3. 30  illustrates  the 

results  for  a 5 degree  wedge  at  = 6.  With  (Ac  / (c  from  Figure  3.29 

^ nt7  V i 

one  notes  that  although  the  viscous  induced  effect  on  the  lift  will  not  exceed  20  percent, 
the  viscous  effect  on  imsteady  aerodynamics  can  be  as  high  as  50  percent  for  aft  eg 
or  elastic  axis  locations. 


As  AP„/P.,  is  known,  the  estimate  of  (Ac_^)  in  Eq.  (3. 77)  can  be  improved  by 

^ ® -1/2  V -1/2 

substituting  (Pg/p^  ) in  the  expression  for  with  (p^/p^)  , 


where 


Seeing  these  large  effects  of  pure  viscous  interaction  one  must  be  concerned  about 
what  the  effects  of  transition  might  be.  It  is  not  within  the  scope  of  the  present  study 
to  investigate  this,  but  it  is  a problem  that  needs  future  attention,  especially  in  view 
of  the  associated  thermal  effects.  The  present  predictions  are  compared  with  Orlik- 
I^’ckeman's  theoretical  and  experimental  data  (Ref.  86)  in  Figure  3.31.  The  agree- 
ment between  the  present  first  order  theory  and  Orllk-Ruckeman's  higher  order  theory 
is  very  satisfying  and  assures  that  the  present  embedded  viscous  perturbation  method 
will  be  sufficiently  accurate  for  the  tactical  missile  imder  study.  Figure  3. 32  shows 
a more  detailed  comparison  between  Orlik-Ruckeman's  theory  (Ref.  98)  and  the 
present  perturbation  method.  The  difference  in  invlscid  characteristics  is  the 
difference  between  piston  theory  used  by  Orllk-Ruckemann  and  the  tangent  wedge 
formulation  used  in  the  present  analysis. 


3. 3. 1. 3 Non-Vlscous  Applications.  The  perturbation  analysis  performed  for 
viscous  flow  can  easily  be  applied  to  other  surface  perturbations,  e.g. , those 
caused  by  elastic  deflections. 
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Mandl  (Refs.  lO:’.— lOH)  has  investigated  numerically  the  effect  of  one  particular  type 
of  surface  deformation  (see  inset  in  Figure  3. 33).  In  this  case  tlie  following  equivalence 
holds. 


d6*/dx  = K6  i (3.79) 

w ' 

2 

where  K = c/R.6  ; k « 1.  K > 0 for  concave  and  K < 0 for  convex  surface 

L W 

deformation.  Again  using  the  embedded  flow  approach,  the  equation  corresponding 
to  Eq.  (3.56),  giving  the  deformation-induced  relative  pressure  increase,  is 


AP„/Pe  = YM^  k C. 


(3.  80) 


and  the  following  equation  corresponding  to  Eqs.  (3.59)  — (3. 61)  is  obtained. 
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The  equation  corresponding  to  Eqs.  (3. 62)  and  (3. 63)  is 
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The  effect  on  the  dynamic  derivatives  is  obtained  equally  simply,  the  equation  cor- 
responding to  Eq.  (3. 66)  being 
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which  after  integration  gives 
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As  the  piston  is  solid  in  this  case  (Ac^J)^  = (Ac^^)j  =0  and  the  total  effect  on 
damping  is  the  - change.  That  is 
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(3.  86) 


Instead  of  moving  the  aerodynamic  center  forward,  as  in  Eq.  (3, 78),  an  equally 
large  aft  moment  is  obtained. 


(3.  87) 


In  Figure  3.33  the  K-derivatives  obtained  from  Eqs.  (3. 80) -{3. 87)  are  compared 
with  Mandl’s  numerical  results  (Refs.  103-105).  The  agreement  is  satisfactory  for 
axis  locations  of  interest,  < 0.75,  and  indicates  that  the  present  perturbation 
method  is  sufficiently  accurate  for  the  small  elastic  deformation  to  be  considered  in 
the  present  aeroelastic  analysis. 
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3.4  CKNTHIPTGAL  PHFSSl'HF  CORRFCTION 


When  the  surface  has  a longitudinal  curvature,  the  pressure  formulas  derived 
earlier  have  to  lx?  corrected  for  the  curvature  effect  (Ref.  41».  With  Itusemann's 
centrifu^jal  pressure  correction  included,  the  N«*\\1onian  pressure  is  (Ref.  Hi 
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e is  the  densitx  ratio  through  the  (U)\\  > shock  and  lt|j(x)  is  the  (longitudinal) 
radius  of  curvature.  For  the  present  analysis  of  thin  airfoils  d - e i ^ 2 ( ; * li 

O 

and  1 R|^(.\i^  ~ ' b Fittaire  3.  Is  is  obtained 

( V|^)j  = h (s)  - h (0)  ± (V  (3.  s!i) 

F I' 

For  a thin  airfoil  the  pressure  ccK'nicient  {riven  bv  Fq.  (3.  ‘'S{  is 

(‘'pInB  (^s)c  \ 

(^p)rs  = (>'* 


Vj_/l'  is  {riven  by  Eq.  (3.  30)  or  Fq.  (3.37).  This  centrifu{ral  pressure  correc- 
tion (AC  ) should  also  lx>  added  to  tlie  tangent  wedge  pressure,  Eq.  (3.35). 

P ^ 
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3.  4. 1 Effect  of  Airfoil  Camber, 


An  infinitely  thin  airfoil  with  circular  arc  camber  is  sketched  in  Figure  3.34.  The 
following  relationships  between  geometric  parameters  hold  for 


According  to  Newlonlan  theor\’  only  the  windward  side  has  non-zero 
pressures.  For  a small  angle  of  attack  Eq.  (3.  89)  becomes 
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And  Eq.  (3.  81)  gives 
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The  0- derivative  is 
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Finalh,  the  effeet  of  camber  on  the  o-derivative  at  0=0  is 
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Thus,  Uie  centrifufjal  pressure  correction  increases  the  Newtonian  pressure 
cun’ature  derivatives  b\  the  factor  |l  -►  (>  +l)  In  the  cas<*  of  the  pure 

N'ewtonian  theoiT,  7-^1  2 , and  Busemann's  eorrection  increases  the 

Newtonian  pressure  bv  50  percent  (see  Ref.  41). 


3.  4.  2 Kffect  of  Surface  Curvature 


If  one  turns  Figure  3.  34  around  180  degrees  and  reverses  the  free  stream  flow’ 
direction,  the  figure  can  represent  the  wlndw-ard  side  of  a circular  arc  (bicusped) 
airfoil.  In  that  case,  Eqs.  (3.  93)  through  (3.  95)  will  represent  the  pressure  deriva- 


tives if  the  sign  for  the  (AC  ) -term  is  reversed.  That  is,  the  strip  nor- 

P ^ I _ 1 


pc  I 

mal  force  derivative  /d^,  is  decreased  by  the  ratio  1 1 - (1+7) 

the  Newtonian  value.  Note  that  2h^/x  represents  the  ( 
between  the  leading  edge  and  station  x (see  Eq.  (3. 91)), 


from 
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Mandl  (Refs.  103-105)  investigated  the  effect  of  constant  surface  curvature  on  the 
aerodynamic  characteristics  of  a wedge-shaped  airfoil  (see  inset  in  Figure  3.33 
for  a definition  of  geometn  ).  In  this  case  Eqs.  (3.  93) -(3.  95)  become  (for  the 
windward  side) 
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(3.  98) 


Again  the  derivatives  at  a = 0 are  increased  by  the  factor  |1  + (y  + 1)  j 
above  the  Newtonian  value  due  to  concave  curvature  effect,  (Convex  curvature 
has,  of  course,  the  opposite  effect,  the  factor  then  being  jl  - (y  + 

For  the  standard  Newtonian  theory  with  Busemann's  centrifugal  correction 
Eq.  (3, 98)  becomes  (Ref.  41). 
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(3.  9 Ha) 


That  is,  tiu'  centrifugal  pressure  correction  increases  the  Newtonian  value  In 
30  percent. 

3.  X'iscous  C'unature  Kffects 

For  the  wedjte  ;mal\  zed  earlier  in  Section  3.3  the  viscous  curxature  parameters 
are  as  follows. 


(3.99) 
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Kquution  90)  gives  with  the  above  Eq.  (3.  99) 
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one  obtains 
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The  magnitude  of  the  first  term  is  always  less  than  6^/VM^  and  the  magnitude  of 
the  second  term  never  exceeds  |2/(yM^j^|  ^Ap^/p^J  (l  + Ap^/p^)  . 

In  the  present  analysis  it  has  been  shown  that  viscld-invlscid  interaction  is  not  of 

9 2 

any  real  concern  unless  M“  « 1 and  6 « 1.  For  6 < 0.2  (11.4°)  and 

e w w ' 

TM  > 5 the  first  term  is  always  less  than  4 percent,  and  the  second  term  would 
not  exceed  6 percent  for  Ap^/p^  < 1.  Thus,  it  appears  that  the  viscous  curvature 
effect  can  be  neglected,  at  least  in  a first  order  analysis  such  as  the  perturbation 
method  developed  earlier. 
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}.  } Klastk-  IVrturlxitions 


Foi-  a thin,  *;\  m metric  airfoil  tlu-  elastic  ilcflcction  is  represented  In  the  defoi  ma- 
tion  ol  the  center  choitl  line  (st>c  I'ipir**  In  this  case  the  followini;  wpiival- 

«‘ncc  with  the  \iscous  analvsis  is  obtaini**! 
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To  this  should  lx‘  added  the  cumatiire  effect  from  l'().  i.'b '>()).  That  is 
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For  a wed^e  describinjj  paraliolic  clastic  deformations  one  can  write  h (.\)  as 
follows 
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The  equations  corresponding  to  Eqs.  <3. 104)  and  (3. 105)  become 
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The  K- derivatives  arc  in  angular  measure 
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That  is,  for  0 = 0 and  K — 0 the  centrifugal  pressure  correction  goes  to  zero. 
Consequently,  for  the  elastic  analysis  in  the  present  study  the  formuiatiou  expressed 
through  Eqs.  (3.34),  (3.35),  and  (3.38)  is  sufficient. 


3.4.5  Surface  Curvature  Perturbations 


For  a wedge  with  a parabolic  surface  curvature  superimposed  isee  inset  in 
Figure  3. 33)  the  following  equivalence  holds 

deVdx  = (3.109) 

2 

K « 1 ; X > 0 for  concave  and  K < 0 for  convex  surface  curvatiu:e. 

The  relative  pressure  increase  given  by  Eq.  (3. 56)  is  for  the  windward  side 

(AP^/PJ  =->'M^Ke„4  (3.110) 
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Throujjh  the  steps  correspondint:  to  those  leading  to  Kq.  (:{.  01)  one  obtains 
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To  this  should  lx*  added  tlie  eentrifugal  pressure  correction  front  Kq.  (3.07). 
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From  Kqs.  (3.  110)  and  (3.  Ill)  one  obtains 
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From  tlie  earlier  discussion  of  Eq.  (3.08)  it  is  ob\ious  that  for  0=0  modifving 
Kq.  (3. 110)  as  follows  will  include  the  centrifugal  pressure  term 
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Integrating  Eq.  (3. 114)  gives  the  normal  force  and  pitching  moment  derivatives 
at  a = 0 
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As  there  is  no  "piston  cushioning"  in  this  case,  the  total  effect  on  the  damping  is 
the  - change.  The  equation  corresponding  to  Eq.  (3.66)  is  as  follows 


2y(y-t-  2) 

y 1 


(• 


(3.116) 


Integration  gives 


(3. 117) 
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In  Fi^re  the  results  from  Kqs.  (3.  115)  :ind  (3.  117)  are  showii  for 

•) 

0^^.  = 1.25  and  M“  1 . Also  shown  arc  the  predictions  throujih  Newtonian, 
Newtonian- liusemann,  and  Mandl's  theory.  The  ajireement  l>etween  present  pre- 
diction and  Newtonian- lUisemann  theory  is  \ ery  );ood,  which,  of  course,  it  should  In' 
for  this  case,  i.c.  , for  — or.  The  centrifugal  pressure  correction  increases 

the  Newtonian  \alue  b\  at  most  20  percent.  The  fipire  shows  that  concaye  surface 
curvature  increases  the  d\  namic  stal)ilit\  for  all  cjx  locations  and  increases  the  static 
stabilitv  for  cp  locations  forward  of  =2  3.  The  effect  of  convex  curvature 

is  exacth  the  opjKOsite.  Whereas  the  apreement  w itJi  experimental  daUi  has  Ix'en 
found  to  Ix'  pood  for  the  Husemann  correction  of  concave  surface  curvature  effects 
(It(  f.  10b),  IK)  experimental  \ erification  of  the  convex  I'urvature  effects  have  been 
obtaiiuHl.  rhus,  in  the  present  anahsis  which  considers  only  convex  airfoil  shapi-s, 
the  centrifupal  pressure  correction  is  neplected. 

This  is  true  also  in  repard  to  the  effect  of  camlxr-inducinp  elastic  deflections, 

W'e  believe  tliat  lx*fore  includinp  surface  curvature  effects  in  the  acroeleastic  analysis, 
efforts  should  lx-  made  to  account  for  the  effects  of  Iioundarv  layer  transition  luid  How 

ff 

separation.  t)nl\  when  the  more  fundamenUil  viscous  flow  effects  have  been  included, 
can  one  justify  the  inclusion  of  higher  order  inviscid  flow  effects,  such  as  the  curva- 
ture effect  tliscussed  alxive  or  the  lx)w  shock  interaction  considered  by  McIntosh 
(Kef.  loT),  Ixtli  being  of  comparable  magnitudes  (sec  Fig.  3.35:mdHef.  lOs^. 


# Regular  viscid-inviscid  interaction  effects  are  insignificant  for  tlie  operating  range 
of  the  tactical  missile  under  study. 
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Figure  3.2  Pressure  distribution  o\er  blunted  slender  iKnlies  of  resolution  ;it  M 


Figure  3. 3 Dynamic  pressure  functicxi  for  blunted  caies  at  = 20 
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Figure  3. 5 Dynamic  Pressure  Functions  for  Blunted  Cones  and  Cylinder- 

Flare  Bodies 
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Figure  i . 6 Dynamic  Pressure  Function  for  Shallow  Blunted  Cones 
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Figure  3. 8 Comparison  between  dynamic  pressure  fimctions  for  blunted  slender 
cones  as  determined  by  analytic  approximation  and  metli^d  of  chax  acteristics 
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Figure  3. 11  Compariscci  between  predicted  and  measured  effects  of 
nose  bluntness  and  mach  number  on  the  dynamic  stability  of  a 7**  cone 

(sheet  2 of  2) 
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Figure  3. 12  Dynamic  pressure  ratio  six  nose  diameters  downstream 
of  the  tip  of  a hemisphere-cylinder  at  = 6 
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Figure  3. 14  Aerodynamic  characteristics  of  a blunt  cylinder-flare 

body  at  » 14 


Figure  3. 15  Wedge  pressure  coefficient  predicted  by  piston  theor>’ 
compared  to  Newtonian  prediction 
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Figure  3. 17  Comparison  of  wedge  pressure  derivatives  predicted  by  tangent  wedge  method 
3rd  order  piston  theory,  and  ganuna-corrected  Newtcxiian  theory 


Figure  3. 18  Definition  of  oscillatory  parameters 
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a.  View  normal  to  surface 
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Figure  3.23  Wedge-side  wall  boundary  layer  interference  (Ref.  78) 
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Figure  3.24  Wedge-wall  comer  flow  (Ref.  90) 
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a.  Viscous  Induced  Forward  Movement  of  tiie  Aerodynamic  Center 


Figure  3. 30  Viscous  effects  on  static  and  dynamic  stability  derivatives 

of  a 5”  wedge  at  = 6 (Sheet  1 of  2) 
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b.  Normalized  Viscous  Induced  Dynamic  Effects 

Figure  3.30  Viscous  effects  <»  static  and  dynamic  stability  derivatives 

of  a 5°  wedge  at  = 6 (Sheet  2 of  2) 
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Figure  3. 31  Comparison  betwem  predicted  and  measured  viscous  effects  on  static  and  dynamic 
stability  derivatives  a 3*  wedge  for  heated  wall  conditions 
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Section  4 

AEROELASTIC  ANALYSIS 


4.1  Approach 

Over  the  past  decade  many  studies  have  been  made  of  the  aeroelastlc  behavior 
of  a wide  range  of  vehicle  configurations  operating  at  hypersonic  speed.  Some 
of  these  are  referred  to  In  Section  1 which  gives  a brief  simnary  of  the  state 
of  the  art  In  the  field.  Here  an  aeroelastlc  analysis  system  Is  presented  for 
lifting  surfaces  on  tactical  missiles.  Previously  existing  codes  for  aero- 
dynamic, thermal,  and  elastic  analysis  with  modifications  as  necessary  have 
been  utilized  as  much  as  possible.  The  system  of  computer  codes  is  intended 
to  provide  a well  tested  analytical  method  for  solving  a wide  range  of  hyper- 
sonic aeroelastlc  problems  for  lifting  surfaces  In  the  Mach  number  range  of 
roughly  3 to  6. 

The  equations  governing  a vehicle  under  flight  conditions  can,  after  reduction 
to  a finite  number  of  degrees  of  freedom,  be  written  in  the  form 

[M]  fu}  + [D]  {u}  + P (u)  + K(u)  + Q(u)  - fF}  (4.1) 


Here  the  vector  u represents  the  "discretized”  displacement  field,  [M]  and 
[D]  represent  structural  Inertia  and  damping  matrices.  The  operators  P 
and  Q determine  the  components  of  aerodynamic  forces  caused  by  structural 
deformation.  In  the  linear  case  the  operator  K corresponds  to  the  stiffness 
matrix.  It  determines  internal  forces  in  terms  of  the  displacement  unknowns. 
Finally,  [F}  is  a vector  representing  other  than  displacement  dependent 
forces  on  the  structure,  including  the  effects  of  aerodynamic  heating.  The 
components  of  fF}  generally  vary  with  time.  Even  when  stationary  flight 
conditions  are  reached,  temperatures  will  continue  to  vary  until  thermal 
equilibrium  is  obtained. 
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A straightforward  numerical  Integration  of  Eq.  (4.1)  over  the  trajectory  or 
over  pertinent  parts  of  It,  Is  In  principle  feasible.  Such  a solution  would 
yield  the  answers  to  all  questions  about  structural  Integrity,  Including 
possible  aeroelastlc  Instability.  However,  the  time  period  over  which  the 
Integration  must  be  carried  out  Is  large  In  comparison  to  the  period  of  free 
vibrations  of  a typical  wing.  Therefore,  the  procedure  will  In  general  not 
be  economically  feasible,  and  even  In  the  presence  of  our  days'  high-speed 
digital  computers,  the  classical  methods  of  aeroelastlc  analysis  survive  as 
viable  options. 

The  force  vector  [F]  Is  generally  a function  that  varies  slowly  with  time 
In  comparison  with  the  free  structural  vibration.  Consequently,  unless  a 
vibration  Is  excited,  the  first  three  terms  In  Eq.  (4.1)  can  be  discarded, 
leaving  us  with  the  equation: 

K(u)  + Q(u)  = {F}  (4.2) 

Here  K and  Q are  nonlinear  algebraic  operators,  containing  terms  of  first 
through  third  degree  In  u.  The  application  of  these  operators  to  a function 
results  In  a vector  with  the  number  of  elements  N equal  to  that  of  the  load 
vector  {F}.  Consequently,  we  can  write 

K^(u)  + q Q^(u)  - F^  1 - 1,  N (4.3) 

For  reasons  that  will  be  clear  from  the  subsequent  discussion,  we  make  the 
substitution 


{“)  * {“q}  {"i) 

in  which  we  assume  that  [u^}  represents  an  equilibrium  configuration  and 
that  {Uj^3  small  that  all  terms  of  higher  than  the  first  order  in 

u^  may  be  discarded. 
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After  substitution  of  Eq.(4,4)  into  Eq.(4.2),we  can  subtract  the  equation: 


K(u  ) + q Q(u  ) - fF) 

O O'-'' 


(A.  5) 


Then  Eq. (4.2) takes  the  form 

[K]  {Ui)  + q LQ]  - 0 (4.6) 

The  elements  of  the  matrices  [K]  and  [Q]  are  in  general  functions  of  the 
displacement  vector  That  is, 

ic,,  - OK./dxJ  at  u - u 
ij  1 J o 

and 

Qij  * OQj^/SXj)  at  u “ u^  . 

For  Eq.(4.7)  to  have  nonzero  solutions 

Det  f[K]  + [Q])  - 0 . 

If  {F}  « 0 (no  thermal  effects,  symmetric  wing  profile,  and  zero  angle  of 
attack),  then  [K]  is  represented  by  the  linear  stiffness  matrix  [K]  , and 
Eq.  (4.8)  represents  the  classical  fornulation  of  the  divergence  problem 
(static  aeroelastlcity) . With  a symmetric  wing  profile  thermal  effects 
cause  only  In-plane  displacement  components.  The  eigenvector  obtained  from 
Eq.  (4.8)  contains  normal  displaceF;'>nt  components  only.  In  that  case  the 
solution  represents  a bifurcation  in  the  equilibrium  path.  The  matrix  [R] 
represents  the  linear  stiffnes  matrix  for  the  prestressed  structure.  That  is, 
it  includes  the  effects  of  thermal  stresses  on  the  stiffness  properties. 

If  the  fundamental  solution  [u]  ■ [u^}  of  Eq.  (4.2)  contains  normal  dis- 
placement components,  the  vanishing  of  the  determinant  (Eq.  (4.8))  generally 
represents  a vertical  tangent  in  a displacement/velocity  diagram  rather  than 
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a bifurcation  point.  In  that  case  the  eigenvalue  approach  to  the  divergence 
problem  is  not  rigorous.  It  is  possible  to  linearize  Eq.  (4.2)  for  computa- 
tion of  fu^}  but  the  solution  obtained  from  Eq.  (4.8)  represents  a more 
or  less  unconservative  estimate  of  the  critical  velocity.  We  refer  to  the 
solution  of  Eq.  (4.8)  as  linear  and  the  direction  solution  of  Eq.  (4.2)  as 
nonlinear  static  aeroelasticity . The  choice  between  the  methods  is  discussed 
in  subsection  4.4. 

Dynamic  aeroelastic  instability  is  said  to  occur  if  the  flight  conditions  are 
such  that  the  amplitude  of  a vibration  will  Increase  with  time.  Solution  to 
that  problem,  of  course,  requires  that  the  "dynamic"  terms  in  Eq.  (4.1)  arj 
included.  The  structural  damping  is  often  discarded  since  its  magnitude  *.s 
not  well  known  and  its  exclusion  results  in  only  slightly  conservative  results 
for  lifting  surfaces  such  as  the  missile  fins  considered  in  the  present 
analysis . 

The  time  step  in  the  Integration  procedure  must  be  relatively  small  in  compari- 
son to  the  period  for  important  vibration  modes.  Since  the  duration  of  the 
trajectory  is  very  large  in  comparison  to  a vibration  period,  the  direct  ap- 
proach is  impractical.  It  is  possible  to  "compress"  the  trajectory  into  a 
shorter  time  period  as  long  as  q and  [F}  remain  slowly  varying.  This 
leads  to  some  difficulties  with  the  thermal  terms.  Better  efficiency  may  be 
achieved  if  the  quaslstatlc  solution  first  is  obtained  (Eq.  (4.2)).  The 
analysis  can  then  be  restarted  with  excitation  of  a small  amplitude  vibration 
at  appropriately  selected  points  along  the  trajectory.  The  trajectory  is 
safe  if  the  vibration  amplitude  is  decaying  at  all  points  along  the  trajectory. 
The  direct  integration  of  Eq.  (4.1)  will  be  referred  to  as  nonlinear  dynamic 
aeroelasticity. 

A linearization  of  the  dynamic  solution  can  be  obtained  in  the  following  way. 

In  Eq.  (4.1)  is  substituted 

{u}  « [u^}  + a^  fUj^j  slniu^t  (4.9) 
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where  fu^)  represents  the  quasistatic  solution  and  a^,  {u^^}  and  fou^) 

represent  amplitude,  deformation  pattern,  and  frequency  of  the  i-th  vibration 
mode.  It  is  assumed  that  the  a^  are  so  small  that  terms  of  higher  than  first 
order  in  a^  can  be  discarded.  After  subtraction  of  the  terms  representing 
quasistatic  equilibrium,  we  obtain  a homogeneous  equation  system  from  which 
frequencies  and  modes  can  be  determined  as  eigenvalues  and  eigenvectors. 

Frequencies  and  generalized  masses  corresponding  to  a finite  nuinber  of  modes 
can  be  determined  at  appropriately  chosen  points  along  the  trajectory.  These 
results  can  be  used  In  a classical  flutter  analysis  as  described  in  numerous 
text  books.  We  refer  to  this  procedure  as  linear  dynamic  aeroelasticity 

analysis . 

The  information  obtained  from  the  linear  elgensolutlon  to  extract  the  flutter 
speed  is  somewhat  limited  since  it  does  not  indicate  the  rate  of  growth  of 
the  vibration  amplitude.  However,  it  does  not  seem  likely  that  a vehicle  can 
be  subjected  to  supercritical  condition  for  a significant  time  interval  with- 
out structural  failur>2. 

It  is  true  also  that  the  linearization  of  the  analysis  precludes  effects  of 
stiffening  of  the  structure  with  finite  deformation.  For  the  cantilevered 
plates  considered  such  structural  stiffening  is  insignificant  within  the 
range  of  tolerable  deformation.  In  most  cases  the  stiffening  of  the  "aero- 
dynamic spring"  will  also  be  Insignificant.  Thus,  it  appears  that  the  linear 
dynamic  aeroelasticity  analysis  is  quite  satisfactory  for  the  present  purpose. 

All  options  discussed  above  have  been  retained  as  capabilities  in  the  com- 
puter program  for  aeroelastlc  analysis,  STAGS-MOFA.  In  the  case  of  linear 
divergence  analysis  the  problem  defined  by  Eq.  (4.8)  is  not  solved  directly. 

The  nonconservative  nature  of  the  aerodynamic  matrices  results  in  a matrix 
that  is  not  symmetric.  Eigenvalue  routines  for  such  problems  are  not  avail- 
able for  problems  of  the  size  that  may  be  required  for  an  accurate  solution 
of  the  lifting  surface  with  thin  leading  and/or  trailing  edges.  As  the 
development  of  such  procedures  could  not  be  contained  within  the  scope  of 
the  present  program,  the  aeroelastlc  analysis  was  based  on  the  modal  approach. 
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That  is,  nonlinear  terms  in  Eq.  (4.1)  are  discarded.  Linearized  equations 
gives  quasi  static  solutions  that  define  the  prestress  for  a free  vibration 
analysis.  The  vibration  modes  and  frequencies  are  used  in  a modal  analysis  of 
linear  divergence. 


The  building  blocks  for  the  aeroelastic  analysis  then  are 


TRRID  for  thermal  analysis 

STAGS  for  elastic  analysis 

FLAM  for  computation  of  aerodynamic  forces 

MODIV  for  linear  static  aeroelastic  analysis 

MOFA  for  linear  dynamic  aeroelastic  analysis 


The  manner  in  which  those  programs  work  together  is  indicated  by  the  Flow 
Chart  in  Fig.  4.1.  Auxiliary  storage  is  used  for  intermediate  results,  such 
as  temperatures  from  TRRID.  However,  by  proper  choice  of  control  cards  it  is 
possible  to  exercise  the  progriims  in  sequence  during  the  same  computer  run. 


The  temperatures  are  f:  ^ t computed  in  TRRID  at  a number  of  time  steps  over 
the  trajectory.  The  results  are  stored  in  a file  and  accessible  so  that  the 
temperature  distribution  at  .ny  fixed  value  of  time  can  be  obtained  by  inter- 
polation. Subsequently,  the  values  of  the  aero,  ^namlc  force  matrices  are  de- 
termined in  FLAM  as  discussed  in  subsection  4 The  information  is  available 
to  STAGS  as  needed  during  a nonlinear  quasi  static  or  nonlinear  transient 
analysis.  The  quasi  static  analysis  gives  the  solution  to  the  nonlinear  di- 
vergence problem  as  well  as  a complete  picture  of  stresses  and  strains  in 
the  structure  during  the  trajectory. 

For  linear  static  aeroelasticity  analysis  the  user  will  ask  STAGS  for  a linear 
analysis.  Modes  and  frequencies  will  be  computed  and  stored  so  that  for  each 
preselected  time  step  a computed  critical  velocity  can  be  obtained  by  use  of 
the  MODIV  program.  The  trajectory  is  safe  if  the  computed  critical  velocity 
at  all  points  exceeds  the  actual  velocity. 
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Unless  the  STAGS  analysis  consists  of  direct  Integration  of  the  nonlinear 
form  of  Eq.  (4.1),  a separate  dynamic  aeroelastlc  analysis  Is  also  needed. 
For  this  purpose  the  user  can  select  to  obtain  either  a linear  or  a nonlinear 
quasi  static  solution.  If  the  user  chooses  the  option  of  nonlinear  dynamic 
aeroeleastlclty  analysis,  the  transient  analysis  will  be  restarted  at  pre> 
selected  steps  with  excitation  of  the  vibration  modes.  If  he  chooses  linear 
dynamic  aeroelastlclty  analysis,  modes  will  be  stored  on  file  and  utilized 
by  MOFA  In  a classical  flutter  analysis.  If  the  computed  flutter  velocity 
exceeds  the  actui<l  velocity  at  all  time  steps,  the  trajectory  is  safe. 

Some  results  obtained  by  the  analysis  system  discussed  here  are  presented 
In  subsection  4.4.  Based  on  those  results  recommendations  for  an  optimum 
procedure  In  using  the  program  will  be  made  In  Section  5. 


I 
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4.2  The  Aerodynamic  Input 

The  present  study  represents  a first  step  toward  the  development  of  a complete, 
integrated  system  for  aeroelastic  analysis  of  flight  vehicles  with  thin  lift- 
ing surfaces  (wings  and  fins)  operating  in  the  lower  hypersonic  regime.  It 
is  limited  to  small  angles  of  attack.  Thus,  much  of  the  basic  theory  is  al- 
ready available  since  only  the  fin  strip  aerodynamics  is  considered.  However, 
the  moderate  nose  tip  bluntness  usually  necessary  in  a hypersonic  missile 
design  generates  a nonuniform  flow  field  over  the  fins.  Analytic  means  have 
been  developed  to  account  for  these  effects,  as  was  discussed  in  the  previous 
section.  The  capability  is  presently  limited  to  hemispherical  nose  tips. 

However,  the  nonuniform  dynamic  pressure  distribution  generated  over  the  fins 
by  such  a nose  tip  can  be  computed  and  Is  Included  in  the  program. 

A significant  simplification  is  introduced  in  invlscid  hypersonic  theory  by 
the  fact  that  the  pressure  at  a point  is  a function  of  the  local  angle  of 
attack  only.  This  reduces  the  aerodynamic  calculations  to  consideration  of 
chordwise  strips.  It  was  shown  in  Section  3 that  for  the  very  small  angles 
of  attack  considered  in  the  present  study.  Third  Order  Piston  Theory  will  pro- 
vide a simple  means  for  calculation  of  the  inviscid  aerodynamic  influence  co- 
efficients with  sufficient  accuracy. 

4.2.1  Calculation  of  Third  Order  Piston  Theory  Aerodynfunic  Influence 
Coefficients 

The  geometry  of  a typical  two-dimensional  strip  will  be  discussed  with  reference 
to  Fig.  3.18  in  which  the  important  variables  are  illustrated.  The  derivation 
follows  closely  that  of  Ref.  109.  It  is  noted  that  the  slope  at  any  point  is 
the  sum  of  five  components;  namely,  initial  angle  of  attack  of  the  airfoil, 
camber,  control  surface  Increment,  variation  in  thickness,  and  structural  de- 
formation. Since  this  study  excludes  rigid  body  motions,  the  first  four  com- 
ponents can  be  determined  at  the  start  of  the  analysis.  All  changes  of  the 
effective  angle  of  attack  are  due  to  flexibility  and  are  the  thermoelastic 
perturbations  in  the  analysis  from  the  rigid  body  steady  state  condition.  With 
the  usual  small  angle  assumption,  Eq.  (3.36)can  be  written 
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V“’b  ■ "u  ■ + ‘f) 


(Vj^/D)^  - + to  + Sj) 


(4.10) 


Here  *nd  are  upper  and  lower  aurface  alopea  with  reipect  to  the 
velocity  vector,  or  ■ + 0 la  tha  rigid  body  angle  of  attack,  and 

6^  “ dh/dx  la  the  change  of  aurface  alope  due  to  atructural  flexibility 
for  a thin  airfoil.^  The  preaaure  difference  normal  to  a lifting  aurface 
is 


( 


^ - <Sl  ■ V 


(A. 11) 


In  which  the  preaaure  coefficient  la  expreaaed  In  terma  of  the  total 

slope  at  a point  as 


(4.12) 


On  substituting  Eq.  (4.12)  into  Eq.  (4.11)  and  distinguishing  between  upper, 
and  lower  aurface  slopes  one  obtains 


f-5 

• CD 

Inserting  Eq.  (4.10)  into  Eq.  (4.13)  and  expressing  the  result  In  ascending 
powers  of  the  slope  (pt  + fi^)  gives 

^ - Qq  + Qi  (a  + 6j)  + Q2  (a  + 6^)^  + Q3  (or  + 6^)^  (4.14) 

in  which  the  are  the  Aerodynamic  Influence  Coefficients  (A.I.C.)  as 

function  of  the  total  slope. 


«f>H 


- <‘f>t 
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% 

^2 

^3 


2 / N 4.  J ^ / 2 2 

00 

^ + (Y  + IHv^  + u„)  + <4^ 

00 

”-  (“l  - “u> 

(V  + 1)  M /3 
00 


) + (uj  - uj) 


Charts  of  these  coefficients  are  given  In  Ref.  109  for  many  profiles.  Expand* 
Ing  Eq.  (4.14)  to  obtain  an  expression  for  the  lifting  pressure  In  terms  of 
the  perturbation  or  flexibility  components  of  the  streamwlse  surface  slope 
gives 


f ■ «0  ^ *£  * ‘f  * 

In  which  the  A.I.C.s  as  functions  of  the  rigid  body  surface  slopes  are 

Qj  - + 2QJ  a + 3QJ  or^ 

Q2  - Q2  + 3Q3  a 

Q3-Q3 

The  are  discussed  In  some  detail  and  given  graphically  for  zero  angle 

of  attack  In  Ref.  109,  showing  that  a significant  Increase  In  complexity  re- 
sults for  a nonsymmetrlc  section  or  nonzero  Initial  angle  of  attack.  In  the 
nonsymmetrlc  case  ^ H terms  In  Eq.  (4.13)  must  be  conq>uted. 

For  a symnetrlc  section,  ^0  " ^2  " ^ * *****  there  Is  no  justi- 

fication for  retaining  Q^.  Thus,  the  aerodynamic  formulation  reduces  to  a 
linear  one. 
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In  what  has  been  said  up  to  now,  no  distinction  has  been  made  between 
dynamic  (unsteady)  and  static  behavior.  All  invlscld  pressures  are  known 
In  terms  of  the  local  slope  and  the  distinction  between  the  two  lies  in  a 
different  Interpretation  of  "slope”.  For  static  flow,  the  chordwise  slope 
is  simply  that  due  to  static  structural  flexibility. 


6 


f 


5x 


(4.16) 


where  h represents  lateral  displacement  and  the  x the  distance  along  the 
chord.  Lateral  motion  Introduces  an  additional  component  equal  to  the 
lateral  velocity  divided  by  the  freestream  velocity,  so  that  the  effective 
slope  becomes 


eff 


dx  u at 

00 


(4.17) 


A discussion  of  generalized  aerodynamic  forces  is  included  here  for  complete- 
ness. 


4.2.2  Computation  of  Generalized  Aerodynamic  Forces 

The  aeroelastic  analysis  as  discussed  above  includes  the  use  of  generalized 
aerodyn£imlc  forces.  It  is  assumed  that  the  lateral  displacement  h , of  a 
thin  platelike  lifting  surface  can  be  defined  in  terms  of  its  first  L 
normal  modes. 

L 

h (x,y,t)  - ^ ^ 0^  (x,y)  (t)  (4.18) 

X-1 
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0^(x,y)  is  le  mode  shape  and  qjj(t)  Is  the  deflection  amplitude  at  a 
reference  point.  Then,  as  is  well  known,  the  equations  of  motion  can  be 
written 

M {?}  + D [5]  + M {§}  * {G}  (4.19) 

in  which  the  x are  the  natural  frequencies  with  M the  generalized 

a* 

mass,  and  the  structural  damping.  The  components  of  the  generalized 

force  fG^j  are 

where  S is  the  area  of  the  lifting  surface  and  F is  the  normal  force 

z 

for  the  case  of  a thin  surface 


F (x,y,t)  - (p  - p ) « Ap  (x,y,t)  (4.20) 

z L U 

where  4p  Is  the  pressure  difference  between  upper  and  lower  surfaces. 

In  a linear  analysis  at  a * 0,  the  only  term  contributing  to  the  lifting 
pressure  is  that  containing  the  first  power  of  6^  , i.e., 

4p  (x,y,t)  * qJJj  (x,y)  6^  (x,y,t)  . (4.21) 


The  streamwise  effective  slope  is  given  by 


6^  (x,y,t) 


+ i 5h 

3x  5t 

00 


(4.22) 


For  simple  harmonic  motion  of  frequency  x , Eq.  (4.22)  can  be  written  in 
terms  of  the  reduced  frequency  k « cx/2U^  . 

. T iiut  /^.  Oil  h.  ixt 

6f(t)  = e = (^  + 2ik  -)  e 

Introducing  matrix  notation  and  expressing  the  displacements  and  streamwise 
slopes  in  terms  of  their  respective  modal  matrices  gives 
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(4.23) 


(6j)  - (’J 

in  which  [Df]  and  [*]  are  the  modal  slope  and  modal  deflection  matrices, 
respectively.  The  former  is  usually  written  as  an  operator  matrix  times  the 
modal  deflection  matrix 


[Dt]  - [D][$]  . 


We  can  now  write  the  aerodynamic  force  vector  as  follows  for  a ^ 0 


(4.24) 


A is  the  diagonal  matrix  of  element  planform  areas.  The  generalized  forces 
are  defined  by  the  vector  [G}  * [Q]  [q]  where 


[Q] 


LQi  I LaI  [m[  + ^ [»] 

LqJ  LaI  L[D]  + ^ "I  [*]  • 


(4.25) 


Using  the  definitions  of  dynamic  pressure  and  reduced  frequencies,  we  have 


7 


2 2 
P C (D 

8r2 


(4.26) 


On  substituting  Eq.  (4.26)  into  Eq.  (4.25)  the  final  matrix  [Q  ] needed  to 
represent  the  generalized  aerodynamic  forces  in  the  flutter  program  can  be 
written  in  a form  consistent  with  that  in  Ref.  32,  namely 


[^l]  - ^ [ [BR]  + ^ [BI]  ] 
8k 


(4.27) 


in  which 


[BR] 


m [»] 
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and 


[BI]  = [f]^  LQj'^  LaI  [»1  . 


These  are  the  two  components  of  the  [Q^]  matrix  found  in  Eq.  (4.1.12)  of 
Ref.  32.  Once  they  are  computed,  along  with  the  generalized  mass  matrix  and 
the  natural  frequencies  of  the  lifting  surface,  these  data  are  fed  into  the 
MOFA  program  (Ref.  32)  to  obtain  the  flutter  speeds  according  to  the  linearized 
formulation. 

4.3  Program  Capabilities 

A computer  program  ICAM  has  been  derived  in  which  aerodynamic  influence  co- 
efficients and  generalized  aerodynamic  forces  are  determined  for  use  in  the 
STAGS  and  MOFA  programs.  The  programs  include  the  chordwise  thickness 
variation  of  wedge,  diamond  and  parabolic  sections,  and  it  is  very  easy  to 
modify  a subroutine  PR0FIL  to  account  for  any  other  thickness  distribution. 

The  Mach  number  is  assumed  to  be  sufficiently  high  so  the  correction  for 
sweepback  (R.jf.  32)  is  negligible. 

4.4  Numerical  Results 

4.4.1  Linear  Static  Aeroelasticity 

Some  classic  divergence  results  have  been  obtained  by  use  of  STAGS  and  MODIV 
in  combination.  With  the  prestress  from  a linear  analysis  STAGS  computes  vibra- 
tion modes  for  use  in  MODIV. 

The  effective^*^  dynamic  pressure  as  a function  of  the  distance  from  the  root  chord 
was  determined  by  use  of  the  method  discussed  in  Section  3.  The  dynamic  pressure 
for  a straight  wing  with  the  leading  edge  42.0  inches  behind  the  spherical  nose 
of  the  missile  body  is  shown  in  Figure  4.2.  In  Figure  4.2a  results  are  shown  for 
one  body  with  a 2.0  inch  and  one  with  a 6.0  inch  diameter.  Figure  4.2b  shows  the 
effect  of  increasing  angle  of  attack  from  zero  to  a * 6°  for  the  6.0  inch  body 
diameter.  In  all  cases  the  effect  of  the  body  reduces  considerably  the  dynamic 
pressure  over  a large  portion  of  the  fin.  Thus,  the  presence  of  the  missile  body 
should  result  in  an  increase  of  the  free  strefun  dynamic  pressure  at  which  aero- 
elastic  instability  occurs. 

# The  distribution  at  the  quarter  chord  of  the  diamond  airfoil  section. 
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The  results  presented  in  Table  4.1  apply  to  a rectangular  wing  with  a 10.0  in. 
chord  and  8.0  in.  span.  The  cross-section  of  the  wing  is  "diamond-shaped"  as 
shown  in  Fig.  2.1.  In  Table  4.1A  are  shown  the  results  for  a wing  that  has 
only  partial  support  along  the  root  chord,  i.e.,  all  displacements  are  con- 
strained at  nodes  7,  8,  9 but  the  other  nodes  are  free  (see  Fig.  2.1).  It 
is  seen  that  the  critical  velocity  increases,  as  expected,  when  the  effect 
of  the  missile  body  is  taken  into  account.  The  slight  decrease  in  critical 
velocity  with  angle  of  attack  is  due  to  the  changes  in  the  generalized  aero- 
dynamic forces.  For  the  wing  alone  this  simply  reflects  the  nonlinear  aerodynamics 
described  by  piston  theory.  In  the  presence  of  the  missile  body  the  effect  is 
reinforced  by  the  Increase  with  y of  the  local  dynamic  pressure^^(see  Fig. 
3.13). 

If  the  same  wing  is  clamped  along  the  entire  root  chord,  the  critical  free 
stream  pressure  increases  dramatically,  as  can  be  seen  from  Table  4. IB. 

Again,  the  presence  of  the  missile  body  results  in  an  increase  of  the  critical 
pressure.  Most  of  the  results  presented  in  Table  4. IB  include  the  effect  of 
aerodynamic  heating.  They  correspond  to  a parabolic  temperature  distribution 
with  450°F  at  the  leading  and  trailing  edges  and  0°F  at  midchord. 

4.4.2  Nonlinear  Static  Aeroelastlcity 

It  was  indicated  in  subsection  4.1  that  the  linear  static  aeroelastic  analysis 
presents  a rigorous  solution  only  if  no  component  of  the  divergence  mode 
(eigenfunction)  is  contained  in  the  precritical  deformation  pattern.  This 
generally  is  the  case  only  if  the  precritical  mode  is  free  from  lateral  dis- 
placement, that  is,  if  the  wing  profile  is  symmetric  and  the  angle  of  attack 
is  zero.  If  a component  of  the  divergence  mode  is  present  at  subcritical 
velocity,  then  this  component  will  grow  nonlinearly  from  the  outset.  Rigorously 
correct  deformations  and  stresses  in  the  structure  can  only  be  determined  from 
the  nonlinear  quaslstatlc  solution.  The  traditional  linear  solution  is  ade- 
quate for  small  dynamic  pressures  but  becomes  progressively  less  credible  as 
dynamic  pressure  increases. 

The  shear  flow  Effect  at  a ^ 0 (Profile  A-A  in  Fig.  1-2)  , which  is  not 

Included,  could  amplify  this  a-effect  significantly. 
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TABLE  4.1 

LINEAR  DIVERGENCE  RESULTS 

RECTANGUUR  PLANFORM  WING  10.  IN.  CHORD,  8. IN.  SPAN.  4.57o  T/C  DIAMOND  SECTION 


A)  3 Node  Root  Support.  Room  Temperature. 


BODY 

ANGLE  OF  ATTACK 

(DEG)  DIVERGENCE 

DYNAMIC  PRESSURE  (LBS/IN^) 

No 

ne 

0. 

410. 

2, 

404. 

4. 

391. 

6. 

370. 

6 

In 

Dia. 

0. 

566. 

2. 

553. 

4. 

510. 

5 

6. 

447. 

B) 

Same  Wing  but  Clamped  Along  the  Root 

• 

BODY 

ANGLE  OF  ATTACK 

L.E.& 

T.E.  TEMP. 

DIVERGENCE  DYNAMIC 

(DEG) 

(DEG.  F.) 

PRESSURE  (LBS/IN^) 

None 

0. 

0. 

2400. 

6 

In 

. Dia. 

0. 

0. 

3550. 

None 

0. 

4! 

50. 

2200. 

2. 

2180. 

4. 

2150. 

6. 

2110. 

6 

In.  Dia. 

0. 

2670. 

2. 

2620. 

4. 

2490. 

6. 

1 

’ 

2260. 
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The  shortcomings  of  the  linear  static  aeroelastic  analysis  were  first  noted 
by  Padlog,  Donato,  and  Batt  (Ref.  19).  With  a simple  but  still  representa- 
tive model  of  a wing  susceptible  to  divergence  it  is  shown  that  the  classical 
linear  analysis  may  give  nonconservative  results.  The  reason  for  this  is 
that  the  nonlinear  contribution  to  the  aerodynamic  forces  is  such  that  the 
twisting  moment  grows  with  the  angle  of  twist  at  a rate  that  exceeds  the 
linear  rate.  Consequently,  the  curve  representing  the  secondary  equilibrium 
form  bends  back  from  the  maximum  at  the  critl  1 dynamic  pressure.  The 
equilibr  t on  this  secondary  path  is  unstable  so  the  situation  is  analogous 
to  the  well  known  and  much  discussed  problem  with  imperfection  sensitivity  in 
elastic  stability  analysis.  Of  course,  the  nonlinear  structural  terms  should 
also  be  accounted  for.  If  the  structure  were  stiffening  under  increasing  de- 
formation, the  situation  could  be  reversed,  resulting  in  a stable  secondary 
equilibrium  path.  However,  at  moderate  deformation  structural  nonlinearities 
are  insignificant  for  the  cantilever  plate.  If  a listurbance  is  present  that 
^ will  cause  a lateral  displacement,  this  displacement  will  grow  at  an  accelerated 

rate,  until  at  some  point  the  slope  of  the  curve  is  infinite.  At  that  point 
the  equilibrium  becomes  unstable  and  the  existence  of  a vertical  slope  may  be 
considered  to  be  a nonlinear  divergence  condition.  The  dynamic  pressure  corre- 
sponding to  aeroelastic  instability  is  below  that  determined  by  the  linear 
theory  and  decreases  with  increasing  size  of  the  disturbance. 

It  should  be  noted  that  the  vertical  tangency  point  may  correspond  to  very 
large  deformations.  Excessive  stresses  may  confine  flight  to  a narrower  range. 
Therefore,  it  seems  that  in  the  majority  of  applications  it  is  necessary  to 
carry  out  a nonlinear  analysis.  In  addition  to  the  need  to  assess  structural 
integrity,  there  is  a need  to  consider  the  effects  of  wing  deformation  on 
flight  mechanics.  In  particular,  the  aerodynamic  stability  derivatives  are 
affected  by  such  deformations.  It  is  demonstrated  in  Ref.  14  that  relatively 
small  aerodynamic  nonlinearities  can  sometimes  produce  dramatic  changes  in 
the  stability  derivatives  at  high  values  of  the  dynamic  pressure. 

For  code  verification  and  to  illustrate  the  relationship  between  linear  and 
, nonlinear  static  aeroelastic  analysis  STAGS  was  again  applied  to  the  rectangular 
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planform  wing  with  diamond  shaped  profile  (Fig.  2.1).  The  same  finite  dif- 
ference grid  was  used  and  the  wing  was  attached  only  at  the  grid  points  7, 

8,  and  9.  No  thermal  or  missile  body  aerodynamic  effects  were  included.  The 
linear  divergence  velocity  was  determined  by  use  of  STAGS  - computer  vibration 
nodes  in  MODIV.  The  critical  dynamic  pressure  "o  computed  is  shown  in  Fig. 
4,j.  The  imperfection  causing  lateral  displacements  was  introduced  in  the 
form  of  a small  lateral  force  at  the  tip  of  the  leading  edge  or  by  flight  with 
a nonzero  angle  of  attack.  Three  values  of  this  perturbation  force  were  con- 
sidered as  shown  in  the  figure.  For  these  three  cases  the  angle  of  attack 
was  zero.  In  a fourth  case  no  imperfections  were  introduced,  but  the  initial 
angle  of  attack  was  2.0  degrees. 

The  stiffness  properties  of  aluminum  was  used  in  the  analysis.  Even  for  a 
high  performance  alloy,  the  stresses  will  be  well  into  the  inelastic  range 
when  the  wing  displacement  reaches  a value  of  0.5  in.  Consequently,  at 

flight  with  an  initial  angle  of  attack  of  2.0°  ^ the  critical  dynamic 
pressure  is  less  than  one-third  of  that  computed  from  the  linear  theory. 

It  may  be  observed  here  that  the  effect  of  the  missile  body  is  to  reduce  the 
dynamic  pressure  on  a large  part  of  the  wing.  For  a given  maneuver  this  may 
have  to  be  compensated  for  by  an  increase  in  the  angle  of  attack.  In  view  of 
decreasing  body  effect  with  increasing  angle  of  attack  and  the  additional 
shear  flow  effect  at  a 0,  the  presence  of  the  missile  body  may  on  the 
balance  have  an  unfavorable  effect  on  the  static  aeroelastic  stability  of  the 
wing  even  in  inviscid  flow.  When  the  viscous  wing-body  interaction  is  con- 
sidered, the  adverse  effect  of  the  presence  of  the  missile  body  will  be  much 
greater,  as  was  discussed  earlier  in  connection  with  Fig.  3.22. 

4.4.3  Dynamic  Aeroelasticity 

The  flutter  speed  was  computed  for  a number  of  cases  by  use  of  the  STAGS-MOFA 
combination.  The  flutter  analysis  in  MOFA  was  based  on  the  three  lowest  fre- 
quencies and  corresponding  modes,  as  computed  by  use  of  STAGS.  As  a first 
example  the  rectangular  wing  with  a 10.0  in.  chord  and  8.0  in.  span  was  con- 
sidered (see  Fig.  2.1).  The  wing  was  clamped  along  the  entire  root  chord.  It 
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was  assumed  to  have  Its  leading  edge  A2.0  in.  aft  of  the  missile  body  nose 
tip.  The  flight  was  assumed  to  be  at  10,000  ft.  altitude  and  at  a Mach 
number  of  5.0.  Results  are  shown  in  Table  A. 2.  As  in  the  static  divergence 
case,  the  ci.ltical  velocities  are  seen  to  Increase  significantly  due  to  the 
missile  body  effect,  and  to  decrease  slowly  with  increasing  angle  of  attack. 

To  assess  the  importance  of  thickness  on  the  aerodynamic  forces  and  the  re- 
sulting flutter  condition  results  were  also  obtained  using  the  "Mach  Box" 

2 

theory  of  Ref.  32.  The  resulting  flutter  dynamic  pressure  of  1600  Ibs/in 

2 

is  markedly  unconservative  relative  to  the  900  Ibs/in  computed  accounting 
for  thickness  effects.  This  evidence  of  the  destabilizing  influence  of  thick- 
ness is  consistent  with  the  majority  of  available  literature,  such  as  Refs. 

110,  111,  112,  and  113. 

Table  A. 3 shows  some  results  for  a clipped  delta  wing.  Planform  and  cross- 
section  of  the  wing  are  shown  in  Fig.  A. A.  The  results  shown  in  Table  A.3A 
do  not  Include  effects  of  aerodynamic  heating.  The  effects  of  the  missile 
body  and  of  angle  of  attack  are  similar  to  those  for  the  rectangular  wing. 

The  slight  decrease  in  flutter  speed  with  angle  of  attack  agrees  with  Refs. 

76,  llA,  115  and  116.  Table  A.3B  shows  results  corresponding  to  a parabolic 
temperature  distribution  with  600°F  at  leading  and  trailing  edges  and  0*^F  at 
midchord. 

The  body  effect  is  greater  for  the  thermally  stressed  clipped  delta  wing  at 
isothermal  conditions.  This  is  due  to  the  fact  that  at  high  temperatures  the 
structural  stiffnesses  are  being  reduced,  tending  to  amplify  the  relative  im- 
portance of  the  aerodynamic  terms  which  are  temperature  invariant.  Thus  we 
nave  reason  to  believe  that  in  general,  the  hotter  the  structure  is,  the  more 
important  do  the  aerodynamic  terms  become. 

As  was  discussed  above,  solution  of  the  nonlinear  equations  of  motions  including 
displacement  dependent  aerodynamic  loads  would  yield  a complece  solution  to  the 
aeroelastic  problem.  While  the  flutter  analysis  only  Indicates  the  velocity  at 
which  a small  vibration  will  grow  with  time,  the  integration  of  Eq.  (A.l)  will 
also  show  the  rate  of  growth  of  the  oscillation.  It  may  show,  for  example, 
that  a maneuver  that  includes  a short  time  at  a speed  slight  above  the  flutter 
velocity  still  is  safe. 
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TABLE  4.2A 

BINARY  FLUTTER  RESULTS -RECTANGULAR  PLANFORM  WING 

10. IN.  CHORD,  8. IN.  SPAN,  4.5%  T/C  DIAMOND  SECTION.  ROOM  TEMPERATURE. 

CLAMPED  ROOT,  MACH  NO.  = 5.0  AT  10,000  FT.  ALTITUDE 
WING  L.E.  AT  42.0  IN.  AFT  OF  BODY  NOSE 


BODY 

ANGLE  OF  ATTACK  (DEG) 

FLUTTER  DYNAMIC  PRESSURE  (LBS/IN^) 

None 

0. 

900. 

2. 

850. 

4. 

800. 

6. 

750. 

6 In. 

Dia. 

0. 

1200. 

2. 

1100. 

4. 

950. 

6. 

850. 
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TABLE  4.3 


BINARY  FLUTTER  RESULTS-CLIPPED  DELTA  PLANFORM  WING 

10.  IN.  ROOT  CHORD,  2.5  IN.  TIP  CHORD,  8.0  IN.  SPAN,  4.5%  T/C  DIAMOND  SECTION 

MACH  NO.  5.0  AT  10,000  FT.  ALTITUDE, 

WING  L.E.  AT  ROOT  CHORD  AT  42.0  IN.  AFT  OF  BODY  NOSE 

a)  Room  Temperature 

BODY  ANGLE  OF  ATTACK  (DEG)  FLUTTER  DYNAMIC  PRESSURE  (LBS/IN^) 


None 

0. 

1900. 

2. 

1880. 

4. 

1840. 

6. 

1810. 

6 In. 

Dia. 

0. 

2840. 

2. 

2730. 

4. 

2530. 

1 

6. 

2340. 

B) 

Leading  and  Trailing  Edges 

at  600  DEGREES  FAHRENHEIT 

BODY 

ANGLE 

OF  ATTACK  (DEG) 

FLUTTER  DYNAMIC  PRESSURE  (LBS/IN^) 

None 

0. 

1020. 

2. 

980. 

4. 

930. 

6. 

850. 

6 In. 

Dia. 

0. 

1720. 

2. 

1620. 

4. 

1500. 

6. 

1300. 
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Data  may  be  saved  from  a quasi  static  analysis  and  STAGS  can  be  restarted  at 
some  points  along  the  trajectory  with  dynamic  terms  Included  together  with 
an  excitation  of  a small  amplitude  vibration.  For  demonstration  of  such  a 
solution  the  rectangular  wing  with  a diamond  profile  was  selected.  At  room 
temperature,  zero  angle  of  attack  and  body  effect  on  the  dynamic  pressure 
excluded,  the  critical  velocity  is  found  to  be  155,000  in/sec.  The  transient 
analysis  with  STAGS  was  started  with  a vibration  excited  by  application  of  a 
load  of  0.1  lbs  at  the  tip  and  the  leading  edge  of  the  wing.  This  load  was 
applied  Instantaneously  and  then  held  constant.  Restarts  were  made  at 
12,000,  140,000,  and  160,000  In/sec  (Fig.  4.5).  It  Is  seen  that  dynamic 

instability  occurs  somewhere  between  140,000  and  160,000  In/sec.  The  oscilla- 
tion Is  not  centered  around  zero  since  the  load  at  the  wing  tip  Is  retained 
and  produces  a nonsymmetrlc  effect  that  grows  with  the  velocity  of  the  free 
air  stream  (see  Fig.  4.3). 
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Figure  4.1  (.hart  Dpmonstratlng  Suggested  Strategy  for  the  Digital  Program. 
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LATERAL  ELASTIC  DISPLACEMENT  AT  LEADING  EDGE  TIP, 


Mach  No.  = 5.0  „ 

Air  Density  ■ 0.75x10  Ibs/soc  /in 


0 100  200  300  400  500  600 

AERODYNAMIC  PRESSURE,  q.,  (psi) 


Figure  4.3  Static  Aeroelastic  Behavior  of  the  Rectangular  Planform, 
Diamond  Section  Wing  Supported  at  Nodes  7,  8 and  9 along 
Root  Chord 
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CORNER  DISPLACEMENT 


Figure  U.b  Transient  Response  Analysis  for  Rectangular  Wing 
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Section  5 

RECOMMENDED  PROCEDURE 


It  is  not  claimed  that  the  work  presented  in  the  preceding  sections  provides 
a completely  satisfactory  solution  of  the  problem  of  aeroelastic  stability 
of  lifting  surfaces  in  hypersonic  flight.  However,  it  is  believed  to  advance 
the  state  of  the  art  in  such  analysis  considerably.  The  computer  program  de- 
veloped under  the  contract  is  somewhat  limited,  and  more  development  is  needed 
before  accurate  analysis  of  a tactical  missile  is  possible.  The  code  contains 
many  options  with  respect  to  the  choice  of  strategy  in  the  aeroelastic  analy- 
sis. These  are  all  left  open  to  the  user.  However,  some  recommendations  arc 
made  for  the  choice  of  procedures.  The  recommendations  are  primarily  based 
on  the  experience  acquired  during  the  development  and  limited  application  of 
the  code. 

The  first  question  that  arises  is  whether  the  linearized  options  may  be  used 
or  whether  it  is  necessary  to  solve  the  complete  nonlinear  equations  (quasi- 
static or  transient).  A few  points  relevant  to  this  question  are  listed  here; 

(1)  It  has  long  been  realized  that  for  thin  lifting  surfaces  a rigorous  defi- 
nition of  stresses  and  deformations  can  be  obtained  only  if  such  nonlinear 
effects  are  included  as  deformation  dependent  loads.  It  has  also  been  shown 
that  nonlinearities  in  the  aerodynamic  formulation  are  important  within  a 
practical  range  of  flight  parameters.  In  a linear  analysis,  the  term  Q(u) 
(Eq.  (4.2))  degenerates  into  ] (u).  It  is  sometimes  assumed 

that  such  an  analysis  gives  a good  estimate  of  deformations  and  stresses  as 
long  as  the  flight  velocity  is  somewhat  below  the  classical  divergence  velo- 
city. For  the  lifting  surfaces  and  flight  conditions  typical  for  hypersonic 
tactical  missiles  this  does  not  appear  to  be  the  case.  Therefore,  unless 
special  circumstances  dictate  differently,  it  is  necessary  to  obtain  the  non- 
linear quasistatic  solution  (Eq.  4.2)).  In  that  case  there  is  no  need  for  a 
solution  to  the  linearized  aeroelastic  problem  as  divergence  will  be  indi- 
cated by  a rapidly  increasing  wing  tip  deflection. 
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(2)  The  nonlinear  transient  analysis  is  somewhat  more  accurate  than  the 
linearized  flutter  analysis  since  it  retains  a much  larger  number  of  degrees 
of  freedom.  In  addition,  it  provides  some  additional  information,  such  as 
the  rate  of  growth  of  the  vibration  modes  in  the  post  critical  condition. 

For  structural  configurations  for  which  the  stiffness  grows  with  increasing 
deformation  the  linearized  analysis  may  be  unduly  conservative.  However, 
these  objections  appear  to  be  slight  in  comparison  to  the  computational  ad- 
vantages that  are  offered  by  the  linearized  flutter  analysis. 

(3)  The  linearized  dynamic  aeroelastlcity  analysis  is  executed  by  MOFA,  using 
frequencies  and  virtual  masses  computed  by  STAGS  and  generalized  aerodynamic 
forces  computed  by  FLAM.  It  is  clear  th.it  within  the  range  of  flight  pa- 
rameters of  interest  thermal  effects  are  important.  It  is  not  equally  clear, 
however,  how  Important  the  effects  of  structural  deformation  on  the  aero- 
dynamic forces  are  (Q(u)  in  Eq.  (4.2))  or  if  aerodynamic  and  geometric  non- 
linearities  need  be  included  in  the  STAGS  analysis.  However,  in  the  general 
case  it  is  necessary,  as  stated  above,  to  obtain  the  nonlinear  quasistatic 
solution.  In  that  case  vibration  modes  are,  of  course,  obtained  for  use  in 
the  flutter  analysis  at  the  time  when  this  solution  is  established.  In  fact, 
the  whole  aeroelastic  problem  can  then  be  solved  in  one  computer  run.  Like- 
wise, whether  or  not  geometrically  nonlinear  terms  ir.  the  prestress  analysis 
need  to  be  retained  is  not  clear  at  this  point.  Such  terms  can  easily  be 
suppressed  by  use  of  the  standard  input  cards  in  STAGS.  Suppression  of  these 
terms  may  lead  to  better  convergence  with  little  sacrifice  in  accuracy.  The 
most  important  effect  of  structural  nonlinearity  is  probably  the  coupling  be- 
tween in-plane  stress  (thermal)  and  torsional  stiffness  of  the  plate. 

(4)  The  capability  to  evaluate  the  linear  divergence  speed  may  be  useful  as 

a first  estimate.  Some  idea  of  when  divergence  is  going  to  occur  is  needed  for 
selection  of  the  time  steps  at  which  vibration  analysis  is  to  be  carried  out. 
The  nonlinear  flutter  option  can  be  used  after  a preliminary  linear  flutter 
analysis.  A restart  with  excitation  of  small  amplitude  vibrations  can  be 
executed  at  a time  step  previous  to  that  at  which  flutter  was  indicated. 
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Based  on  these  arguments,  it  seems  that  the  following  approach  can  be  recom- 
cianded : 

o Perform  a complete  nonlinear  quasistatic  analysis 
covering  a pertinent  part  of  the  trajectory. 

o At  a number  of  intervals  obtain  vibration  modes  and 
frequencies  from  STAGS  for  use  in  a linearized 
flutter  analysis  with  MOFA. 

o If  the  aeroelastic  characteristics  of  the  lifting 
surface  is  relatively  well  known,  this  may  be 
achieved  in  one  run.  If  not,  it  is  recommended 
that  some  preliminary  analysis  is  made,  such  as  a 
linear  divergence  or  flutter  analysis  based  on 
linear  prestress  and  an  approximate  temperature  field. 

o Before  an  extensive  study  (parameter  variation)  is 
performed,  it  may  be  advisable  to  perform  some  com- 
parative analyses  to  determine  if  structural  and/or 
aerodynamic  nonlinearities  may  be  neglected  in  the 
analysis. 

The  problem  of  leading  or  trailing  edge  buckling  deserves  some  consideration. 
Presently  it  seems  that  the  beat  recommendation  is  that  it  be  treated  with 
benign  neglect.  Buckling  caused  by  controlled  displacement,  as  is  the  case 
if  stresses  are  thermally  Induced,  is  generally  a stable  phenomenon.  The 
gradually  increasing  amplitude  of  the  buckle  relieves  the  compressive  stress. 
In  Section  2 it  was  observed  that  the  development  of  leading  edge  buckling 
had  little  effect  on  the  basic  stiffness  of  the  plate.  It  seems  likely  that 
this  behavior  is  typical  and  that  the  observation  can  be  extended  to  similar 
configurations.  Two  different  procedures  for  assessing  thermal  buckling 
effects  are  available.  One  is  to  introduce  a small  perturbation.  In  that 
case  the  quasistatic  solution  corresponds  to  stable  equilibrium  for  all  levels 
of  the  temperature.  In  some  temperature  range  the  lowest  vibration  modes  will 
include  the  very  localized  modes  corresponding  to  leading  edge  buckling.  In 
view  of  the  results  discussed  above,  it  is  not  likely  that  the  presence  of 
these  modes  will  significantly  affect  the  flutter  velocity. 
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If  no  perturbation  is  introduced  to  trigger  thermal  buckling,  convergence  will 
be  considerably  easier.  In  that  case  the  quasistatic  solution  obtained  beyond 
the  point  of  thermal  buckling  will  represent  unstable  equilibrium.  The  coeffi- 
cient matrix  used  in  the  vibration  analysis  has  at  least  one  negative  root,  and 
the  corresponding  eigenvalue  is  negative.  Such  solutions  to  the  vibration 
analysis  are  ignored  in  the  aeroelastic  analysis.  If,  as  expected,  leading 
edge  buckling  has  little  effect  on  the  flutter  characteristics,  the  two  ap- 
proaches should  give  very  similar  results.  There  is  presently  insufficient 
experience  available  on  which  to  base  more  definitive  recommendations  regarding 
the  problem  of  leading  edge  buckling. 


Section  6 
CONCLUSIONS 


A study  of  the  aerothermoelastic  characteristics  of  a tactical  missile 
operating  at  near  zero  angle  of  attack  in  the  Mach  number  range  3 ‘ - 6 

has  produced  the  following  results: 

o An  integrated  digital  program  has  been  developed  for  the  static  and 
dynamic  linear  and  nonlinear  aeroelastic  response  of  sharp  edge.aero- 
dynamically  heated,  cantilevered  lifting  surfaces  representative  of 
a wide  range  of  tactical  missile  wings  and  fins.  The  program  utilizes 
a one-dimensional  "through  the  thickness"  flat  plate  solution  for  com- 
putation of  the  structural  temperature  effects,  third  order  piston 
theory  aerodynamics  modified  to  account  for  the  nonuniform  flow  gene- 
rated by  a hemispherical ly  blunted  missile  body,  and  the  STAGS  non- 
linear finite  difference  program  for  the  structural  (response)  calcu- 
lations . 

o In  general  it  is  found  to  be  desirable  to  run  a full  nonlinear  quasi- 
static analysis  for  each  trajectory  till  its  completion  or  the  onset 
of  divergence.  Flutter  analyses,  linear  or  nonlinear,  should  be  per- 
formed only  at  selected  discrete  times  within  each  trajectory  using 
the  appropriate  prestress  conditions  from  the  quasistatic  solution. 

The  nonlinear  flutter  solution  provides  a direct  numerical  solution 
of  the  coupled  nonlinear  equations  but  is  in  its  present  form  markedly 
less  efficient  than  the  eigenvalue  solution  obtained  through  standard 
type  linear  flutter  analyses. 

Important  results  obtained  in  the  analyses  are: 
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o The  most rue tural  analyses  conducted  to  compare  with  existing  (closed 
fom)  solutions  show  excellent  agreement  where  such  could  be  expected 
and  the  numerical  results  follow  well  established  trends  where  the 
problems  become  too  complex  for  closed  form  solutions.  Computed  re- 
sults illustrate  the  themally  Induced  stiffness  changes  in  wings  and 
demonstrate  the  Importance  of  boundary  conditions,  temperature  sensi- 
tive elastic  moduli,  Initial  imperfections,  planform  and  thickness 
variations . 

o A simple  analytic  theory  has  been  developed  that  can  predict  the  non- 

uniform  flow  characteristics  produced  by  spherical  nose  bluntness  on 

slender  axisymmetric  bodies  at  supersonic  and  hypersonic  speeds, 

M ^3.  It  is  shown  that  nose  bluntness  has  adverse  effects  on  the 
00 

rigid  body  vehicle  dynamics.  The  theoretical  predictions  are  in 
excellent  agreement  with  available  experimental  data. 

o Missile  nose  bluntness  generates  a large  spanwise  dynamic  pressure 
gradient  over  the  fin,  and  will  at  nonzero  angle  of  attack  also 
create  an  inviscid  shear  flow  which  makes  the  aerodynamic  forces  on  the  fin 
displacement  dependent.  Cylindrical  leading  edge  roundness,  however, 
is  found  to  have  insignificant  influence  on  the  two-dimensional  un- 
steady fin  aerodynamics. 

o A simple  viscous  perturbation  theory  has  been  developed  which  gives 
results  that  agree  well  with  available  numerical  and  experimental  re- 
sults. It  is  found  that  for  flight  conditions  of  present  interest  viscld- 
inviscid  interaction  effects  on  the  fin  are  small  in  the  absence  of 
boundary  layer  transition  and  flow  separation. 

o Although  boundary  layer  transition  and  flow  separation  will  have  sig- 
nificant effects  on  the  two-dimensional  aerothermoelastic  charac- 
tistics  of  the  fin  (strips),  it  is  the  three-dimensional  flow  separation 
caused  by  the  interaction  between  the  fin  and  the  fuselage  boundary 
layer  that  is  of  most  concern.  It  can  cause  negative  aerodynamic 
damping  of  movable  fins  for  typical  axis  locations,  and  will  in  general 
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destroy  the  applicability  of  the  simple  piston  theory,  as  the  aerodynamic 
forces  no  longer  are  dependent  only  upon  local  deflections.  This  should 
be  considered  when  interpreting  results  from  the  aeroelastic  analysis  with 
the  present  version  of  the  program. 
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Section  7 


RECCMIENDATIONS  FOR  FUTURE  'ORK 


The  present  study  represents  a first  step  toward  development  of  the  needed 

analytic  capability.  A tactical  missile  has  to  perform  high  g maneuvers. 

Thus,  further  analytic  efforts  are  needed. 

In  the  field  of  unsteady  aerodynamics  the  following  analytic  efforts  would 

appear  to  be  most  urgent: 

o Define  the  aerodynamic  environment  of  the  fins  in  the  angle  of  attack 
range  of  interest  for  a maneuvering  missile.  The  highly  nonlinear 
aerodynamic  characteristics  can  be  defined  by  further  extension  of 
the  embedded  Newtonian  analysis  technique.  The  effects  of  elastic 
perturbations  of  small  amplitude  can  often  be  determined  by  local 
linearization  techniques,  such  as  were  demonstrated  in  the  present 
study.  However,  based  upon  the  results  in  Ref.  14  the  nonlinear  aero- 
structural  coupling  effects  can  be  very  large  and  a full  nonlinear 
analysis  may  be  necessary. 

o Extend  present  analytic  capability  to  include  realistic  missile  body 
geometries  by  considering  biconlc  nose  tips  and  ogival  forebody  geo- 
metries in  addition  to  the  hemisphere-cylinder  and  hemisphere -cone 
geometric  investigated  in  the  present  study. 

o Extend  the  developed  analytic  theory  for  viscid-inviscid  interaction 
to  include  the  effects  of  boundary  layer  transition.  Asymmetric 
transition  can  have  a significant  effect  on  the  fin  aerodynamics  even 
when  thermal  effects  are  ignored,  and  can  be  expected  to  have  a 
greater  influence  when  the  full  aerothermoelastic  coupling  is 
inc luded. 
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o Develop  means  for  determining  how  trailing  edge  flow  conditions  and 
leeside  flow  separation  affect  the  unsteady  two-dimensional  airfoil 
aerodynamics  of  a finite  chord  fin. 

o Conduct  an  investigation  of  the  interaction  between  the  fin  and  the 
fuselage  boundary  layer  which  will  lead  to  a better  unstanding  of 
the  "real  life"  unsteady  fin  aerodynamics.  The  goal  is,  of  course, 
to  develop  means  by  which  the  separated  flow  effects,  which  always 
will  be  present  for  realistic  fin  thicknesses  (Ref.  117)  can  be  in- 
cluded in  an  aerothermoelastic  analysis.  Due  to  the  complexity 
of  the  flow  field,  this  can  be  accomplished  only  through  a powerful 
combination  of  analytic  and  experimental  techniques.  An  experimental 
approach  to  this  problem,  which  also  could  define  the  nose  bluntness 
generated  entropy  wake  effects,  is  discussed  in  Ref.  87. 

o Develop  means  for  determination  of  the  rigid  body  unsteady  aerodynamics 
in  the  angle  of  attack  range  of  interest  for  a maneuvering  tactical 
missile.  In  the  inviscid  flow  field  the  difficulties  are  caused  by 
the  large  nonlinear  effects  of  nose  bluntness.  It  is  believed  that 
the  embedded  Newtonian  theory  can  be  extended  to  describe  these 
effects  for  realistic  missile  geometries.  In  addition  to  this  invis- 
cid flow  investigation  it  is  important  that  some  means  be  developed  to 
assess  the  potentially  large  effects  of  viscous  cross  flow  on  the 
vehicle  dynamics,  in  particular  the  cross  flow  effects  on  boundary 
layer  transition  and  (near)  base  flow  recirculation. 

o It  may  be  desired  that  all  the  fin  problems  discussed  above  are  con- 
sidered also  for  canard  surfaces.  Because  of  the  close  proximity  to 
the  nose  tip,  the  canard  surfaces  will  have  their  own  unique  aero- 
thermoelastic  problems.  In  addition,  their  presence  will  change  the 
flow  environment  experienced  by  the  fins.  Although  the  full  coupling 
between  the  flow  fields  cannot  be  considered  for  some  time,  some  pre- 
liminary, exploratory  investigation  probably  should  be  made  at  a rather 
early  stage  to  get  an  assessment  of  how  large  these  coupling  effects  are. 
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In  the  field  of  structural  and  aeroelastic  analysis  a number  of  tasks  appear 
to  be  of  considerable  interest.  Among  the  more  Important  are  the  following: 

o Exercise  the  developed  capability  in  the  field  of  hypersonic  aero- 
elasticity  to  Improve  our  general  understanding  of  the  problems  en- 
countered and  give  some  direction  for  continued  research  in  the  area. 

The  following  items  may  be  mentioned  as  worthy  of  special  attention: 

. Effects  of  leading  edge  buckling 

. The  accuracy  of  the  aeroelastic  analysis  as  a function  of  the 
number  of  modes  Included  in  Che  analysis 

. Effect  of  structural  nonlinearity  (in  kinematic  relations) 

. Effect  of  wing  plan  form  and  profile  on  the  aeroelastic 
characteristics , 

o In  the  completed  study  the  lifting  surface  was  considered  fixed  at  the 
root  in  one  way  or  another  so  that  all  rigid  body  motions  are  con- 
strained. Ultimately  it  may  be  deemed  necessary  Co  Include  in  the 
analysis  the  rigid  body  modes  of  the  vehicle.  Before  a substantial 
effort  along  chose  lines  is  planned,  it  seems  prudent  to  acquire 
more  experience  by  exercising  presently  available  computer  codes  and 
consulting  with  experts  in  pertinent  areas. 

o Modification  of  the  extrapolation  method  for  initial  estimates  of  the 
unknowns  can  lead  to  considerable  improvements  in  computer  time  for 
the  quasi-static  solution.  By  reorganizing  the  procedures  for  com- 
puting partial  derivatives  of  the  strain  energy  with  respect  to  the 
discrete  variables,  considerable  reduction  in  computer  time  can  be 
achieved.  It  is  expected  that  the  computer  time  for  a typical  tran- 
sient response  analysis  may  be  reduced  by  a factor  from  two  to  five, 
depending  on  Che  problem  parameters. 
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o The  thermal  gradients  of  the  leading  and  trailing  edges  of  a typical  wing 
structure,  as  based  on  solutions  of  the  one-dimensional  heat  equation, 
are  extremely  sharp.  It  appears  that  a much  more  favorable  distribution 
would  be  obtained  if  chordwise  heat  conduction  were  accounted  for.  There- 
fore, a computer  code  should  be  developel  and  integrated  in  STAGS  that  is 
based  on  the  more  realistic  two-dimensional  solution. 

o Solving  the  transient  nonlinear  response  requires  a large  (3  times  the 
number  of  nodes)  order  system  of  coupled  differential  equations  and  is 
expensive  in  computer  time.  Alternative,  more  economic  methods  should 
be  developed  to  handle  the  nonlinear  flutter  problem.  An  obvious  con- 
tender for  this  task  is  the  method  of  modal izing  the  equations.  The  ap- 
proach here  will  naturally  be  more  complex  than  the  linear  procedure 
used  in  the  past.  This  is  due  to  the  coupling  of  the  equations  introduced 
by  the  nonlinear  terms,  which  also  causes  the  number  of  differential  equa- 
tions to  be  solved  to  exceed  the  number  of  normal  modes  chosen.  Neverthe- 
less, for  practical  wing  structures  the  reduction  of  the  number  of  equations 
will  be  many-fold  and  the  resulting  savings  in  computer  costs  is  signifi- 
cant. It  is  also  important  to  note  that  the  linear  eigenvalue  method  is 
subject  to  some  approximations. 

. The  displacement  pattern  is  restricted  so  th'^t  it  can  be  de- 
fined as  the  sum  of  displacement  modes.  Generally,  the 
flutter  analysis  is  restricted  to  four  modes  or  less. 

. The  solution  gives  the  point  in  time  at  which  an  infinitesimal 
disturbance  will  begin  to  grow. 

. The  analysis  does  not  indicate  the  rate  of  growth  of  the 
amplitude  of  the  vibration. 

If  a modal  transient  analysis  is  developed  in  which  the  nonlinear  terms  are 
retained,  the  last  two  approximations  are  eliminated.  This  is  not  believed 
in  itself  to  be  too  important.  However,  if  the  nonlinear  transient  analysis 
is  reduced  by  the  modal  approach  to  entail  only  a few  degrees  of  freedom, 
the  direct  solution  of  the  equations  of  motion  (Eq.  4.1)  will  probably  be 
more  efficient  than  the  linearized  approach.  Hence,  no  price  is  paid  for 
the  additional  information  yielded. 


7-4 


o In  Its  present  form  the  computer  program  does  not  provide  a complete  output 
of  all  the  quantities,  e.g.,  the  aerodynamic  derivatives  of  concern  to 
aeroelasticians . It  is  highly  desirable  that  the  program  be  modified  to 
post-process  the  data  in  order  to  extract  the  maximum  information  from 
each  analysis. 
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Appendix  A 
NOMENCLATURE 


{1} 

a 

a. 

1 

a,  b 
[Al 
Aj,  A^ 

AA 

A,  B 

B,  Bj,  B2 
[BI],  (BR] 
c 

C* 

D 

(D) 


d 


N 


d 


B 


unit  vector 
speed  of  sound 

til 

deformation  amplitude  in  i*"  vibration  mode 
parameters  defined  in  Eq.  (3.16) 
surface  area  matrix  (Eq.  (4.24)) 
constants  defining  f*  (Eq.  (3.6)) 
surface  area  element 
parameters  in  Eq.  (3.16) 
viscous  parameters  defined  in  Eq.  (3.57) 
matrices  defined  in  Eq.  (4.27) 

reference  length,  cylinder  caliber  or  base  dirmeter  for  bodies, 
chord  length  for  air  foils 

Chapmen-Rubesin  parameter  based  on  reference  temperature 
plate  stiffness,Et  (12  (1  ~V^)] 
damping  matrix 
nose  drag 

N 

nose  drag  of  airfoil  strip 
base  diameter 


A-1 


N 


[D*] 

E 


i,  F 
f* 


F 

z 


{F} 

{^zl 

R* 

G^,  C 
h 


i 

k 

kr  and  k 
f R 

K,  K 


K 


, k (7) 


f 

M 

[M] 


diameter  of  axisjanmetric  nosebluntness 

l)ase  height  of  wedge 

diameter  of  cylindrical  leading  edge 

modal  slo|)e  matrix 

Voung's  Modulus 

pressure  correlation  functions 

2 2 

d\Tiamic  pressure  ratio,  pU 
lift  component  (Eq.  (-4.20)) 

total  force  v'ector 
aerodynamic  force  vector 

velocity  ratio,  V 

parameters  defined  in  Eqs.  (3.61)  and  (3.77) 

surface  displacemcat 
imaginary  unit,  i 2 = _i 
reduced  frequency,  cw  2U^ 

parameters  in  Eq.  (3.16)  and  Eq.  (3.17)  respectively 

nonlinear  algebraic  stiffness  o|xjrators  (Section  4) 
hypersonic  similarity  parameters  (Eqs.  (3.23)  and  (3.33)) 
shock  parameters  (Eqs.  (3.4)  and  (3.9)) 

sharp  cone  body  length 
Mach  Number,  U/a 
mass  matrix 


A-2 


m 

P 

M X N 
N 


n 


P 


pitching  moment 

Cm  = Mp/(p^U«2/2)  Sc 

pitching  moment  for  an  airfoil  strip 

c = m /(P^U„^/2)  c^ 

size  of  matrix 
normal  force 

Cn  = /2)  S 

normal  force  for  an  airfoil  strip 

"=n  = « c 

static  pressure 

Cp  = (P  -P,)  / lP^vJ/2) 


A-3 


blast  wave  pressure 

Cp  = <p„  - U^/2) 

P aerodynamic  operator 

)p[  pressure  force  vector 


Pr 

q 

q 

Q,  Q 

II 

n 

% 

It 


s 

t 

t 

T 


■p* 


AT 


crit 


Prandtl  number 

rigid  body  pitch  rate  (Section  3) 

deflection  amplitude  in  {.*th  mode 

dymamic  pressure,  p U^/2 
aerodynamic  operators 

aerodynamic  influence  coefficients  (Eq.  (4. 14)) 
body  radius  (Fig.  3. 1) 

radial  distance  from  bow  shock  centerline  (Fig.  3. 1) 
thermal  forces 

longitudianl  body  radius  of  curvature 
Ueynolds  number,  Ue^  = x U^/i  ^ and 

Re  = X U^/i' 

reference  area, it  c^/4  for  bodies  and  projected  wing  area  for  air  foils 
time 

plate  thickness  (Eq  (2. 5)) 
temperature 
reference  temperature 
critical  temperature  difference 


A-4 


u 

U 

u 

V 

w 

X 

X,  y,  z 


^X 

X,  Y 
y/^N 

Xf  S'!] 


displacement  field  (Section  4) 
deformation  pattern  in  1th  vibration  mode 

strain  energy  (Eq  (2. 1)) 
axial  velocity 

total  potential  energy  (Eq  (2. 1)) 
velocity  normal  to  body  surface 

airfoil  sink  velocity,  normal  to  free  stream 
lateral  load  applied  at  the  wing  tip 
axial  coordinate  (Fig.  3. 1) 
cartesian  coordinates  (Eq  (2.4)) 
global  coordinates  in  STAGS 

centroid  location  of  equivalent  nose  bluntness  (x^  = 0 for  spherical 
or  cylindrical  nose  bluntness) 

cone  center  of  gravity  location  forward  of  base 

displacement  component  (Eq  (2. 1)) 

surface  coordinates  (Eq  (2. 4)) 
spanwiso  fin  coordinate  (Fig.  3. 13) 

distance  from  chord  centerline  to  lower  and  upper  airioil  surface 
(Fig.  3. 18  a) 

wind-fixed  surface  coordinate  (Fig.  3. 34) 
translatory  coordinate  (Fig.  3. 1) 


A-F 


(V 

a 

■y 

•> 

6 

6f 

C* 


V 

i 

1. 


thermal  expansion  coefficient 

angle  of  atUick 

trim  angle-of-attack 

angular  body  coordinate  (Fig.  3.1) 

constant  in  Eg  (4.12) 

ratio  of  specific  heats  ( T =1.4  for  air) 

local  flow  inclination  to  airfoil  surface.  For  small 
inclination  angles  6 = Vj_/U 

change  in  6 due  to  deformation 

boundary  layer  displacememt  thickness 

lx?nding  deformation 

densih  ratio  across  benv  shock 

body  perturlxition  in  pitch 

cone  half  angle 

flare  imgle 

wedge  half  angle 


dimensionless  surface  curvature,  c/R  jj 


0 


w 


K for  airfoil  center  chord 


nose  bluntness  parameter  (Eq.  (3.9)) 

M:ich  angle 
Poisson’s  ratio 

kinematic  viscosity  of  air  (Section  3) 

viscous  hypersonic  similarity  parameter,  P = X 

dimensionless  x-coordinatc,  ^ =x/c 


temperature  parameter  in  Eg(2.5) 
coordinates  for  grid  generation  in  STAGS 


A -6 


shape  parameter  for  temperature  distribution  (Eg  (2.5)) 
amplitude  of  vibration  mode  (Eg  (4.19)) 
air  density 

angulcr  surface  coordinate  (Fig  3.5) 
surface  inclination  to  free  stream  velocity 
azimuth  (Fig.  3.1) 
mode  shape  of  f'th  mode, 
modal  displacement  natrix 

inviscid  hypersonic  similarity  parameters  defined 
in  Eg,  (3.7  and  (3.11)  respectively. 

\iscous  h.Npcrsonic  similarity  parameter  defined  in  Eg.  (3.57) 
natural  Wbration  frequency  of  i’th  mode 
work  done  by  external  forces  (Eg.  (2.1)) 


V. 


I 


Subscripts 


A 

AC 

b 

n 

B 

c 

crit 

CG 

d 


e 


K 


f 

F 

FP 

L 

LK 

i 

i 

M 

N 

N and  NEWT 

7N 

PT 


aft  half 

aerodynamic  center 
body  surface 
base 

Buscmann  correction 
centrifugal  pressure  corrrection 
critical 

center  of  graviW  or  oscillation  center 

deformation 

boundary  layer  edge 

total  elastic  deformation  effect 

flexibility  effect 

forward  half 

flat  plate 

lower  surface 

leading  edge 

inviscid  flow  (Seetton  3) 

initial  imperfection  (Section  2) 

main  airfoil 

nose 

Newtonian  value 
y - corrected  Newtonian  value 
piston  theory 


A-8 


sh 


bow  shock 


SHARP  sharp  leading  edge  or  pointed  nose 

t total  or  stagnation  value 

TE  trailing  edge 

TW  tangent  wedge 

U upper  surface 

V viscositj’  effect 

w wall 

o initial  value  or  constant  parameter  value 

0,  1,  2.  3 parameter  values  bounding  different  regions  ( Eq.  (3.  IG)) 

® free  stream  conditions 

2D  two-dimensional  value 

Superscript 

' prime  superscript  indicates  two-dimensional  value 

- barred  sjTnbols  indicate  absolute  \’alue;  c.g.,  h - |h| 

Derivative  S\-mbols 


e 

ac 



Cv 

■ da  ' 
dC 

m 

c 

m^ 

1 II 

n ^ 

m 

q 

' 

L/ 

m^ 

c 

= C + 

"'q 

C 

c 

Pk 

dK 
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APPENDIX  B 


UTILIZATION  OF  THE  STAGS  HYPERSONIC 
AEROELASTICITY  PACKAGE 


B.  1 Program  Organization 

In  order  to  perform  a complete  hypersonic  aeroelastic  analysis 
of  an  "aerodynamic"  structure,  such  as  a wing  or  a wing-fuselage 
assembly,  two  approaches  are  currently  available.  The  first  approach 
is  to  perform  a STAGS  nonlinear  "quasi- static"  analysis  including  the 
effects  of  both  aerodynamic  heating  and  transverse  forcing  functions, 
while  at  the  same  time  computing  and  saving  vibration  modes  at  selected 
pre-s.tress  states  along  the  flight  trajectory.  In  this  way  a determination 
of  divergence  and/or  thermal  buckling  conditions  may  be  obtained  directly, 
while  a prediction  of  flutter  behavior,  i.e.,  critical  speed,  may  be  obtained 
in  a subsequent  modal  analysis  with  a post-processor.  The  second  approach 
is  to  include  inertial  effects  at  the  outset  by  performing  transient  (time 
integration)  analysis.  This  of  course  obvial  the  need  for  the  modal  flutter 
post-processor  but  at  a significantly  greater  expense  to  the  user.  In  either 
case,  t%vo  external  programs,  HYPATH  (Hypersonic  Aerodynamic  Thermal 
Analysis)  and  HYPALM  (Hypersonic  Aerodynamic  Load  Matrices)  should 
be  executed  in  preparation  for  STAGS,  while  the  optional  "post"-processor; 
FLAM  (Flutter  Analysis  Modal)  may  be  executed  when  using  the  first 
approach.  The  integration  of  these  programs,  which  is  based  on  a 
peripheral  data  structure,  as  well  as  their  individual  input/output  require- 
ments and  those  adjustments  additionally  required  by  STAGS  for  such  an 
analysis  are  described  in  the  following  sections.  The  general  organization 
is  summarized  in  Figure  B.  1. 
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B.2  Program  HYPALM  (Hypersonic  Aerodynamic  Load  Matrices) 


B.2.1  Function 

Prepares  for  a STAGS  aeroelastic  analysis  by  creating  an 
aerodynamic  load  matrix  file,  QTO  , QTl,  QT2,  QT3,  for  each 
"aerodynamically  active"  branch.  These  matrices , in  addition  to  the 
dynamic  pressure  and  elastic  rotations  of  the  structure  (i.e.,  wing) 
are  incorporated  by  STAGS  in  the  formation  of  the  aerodynamically 
equivalent  forcing  function. 


B.  2.  2 Data  Card  Input 

GENERAL  CONTROL  CARDS. 
A-1  Case  Title  Card 


Variable 


Label 


Format 


12A6 


Columns 


1-5 


B-1  I/O  Summary  Card 


Variable 


1ST  EE  R 


Format 


15 


Columns 


1-5 


1ST  AGS 


15 


6-10 


Description 
Case  title. 


Description 

0 - Thickness /chord 

ratios  at  span -wise 
stations  to  be  input 
from  cards. 

1 - Thickness /chord 

ratios  to  be  defined  in 
user  written  subroutine 
PSTEER. 

0 - Print  results  only. 

1 - Save  QT  matrices  on 

tape  for  STAGS 
analysis. 


B-2 


Variable 

Format 

Columns 

Description 

NBRNCH 

15 

11-15 

Totsd.  number  of  branches 
in  STAGS  model. 

B-2  Active  Branch  Summary  Card 

Variable 

Format 

Coltunns 

Description 

LAB  (I) 

(1=1,  NBRNCH) 

1615 

1-80 

List  of  aerodynamically 
active  branches. 

-1  - inactive,  do  not 
perform  analysis  on  branch. 
+ 1 - active,  perform  analys: 
on  branch. 

THE  REMAINING  CARDS  ARE  REPEATED  ONCE  FOR  EACH  "ACTIVE" 
BRANCH  AS  DEFINED  ON  THE  B-2  CARD. 

GEOMETRY 

C-1  Shell  Type  Definition  Card 

Presently,  the  shell  geometry  for  each  active  branch  is  restricted 
to  that  of  a quadrilateral  plate,  as  representative  of  a wing.  Due  to  the 
nature  of  the  analysis  the  quadrilatercil  is  further  restricted  such  that  two 
of  its  sides  must  be  parallel  to  the  local  x axis  (see  description  of 
coordinate  systems  in  STAGS  user's  manual).  This  ensures  a unique 
"span-wise"  station  for  each  column  in  the  branch. 


Variable 

Format 

Columns 

Description 

NSHELL 

15 

1-5 

Basic  shell  geometry. 
3 - Plate. 

NCART 

15 

6-10 

0 - Rectangular  grid. 

1 - Isoparametric  mapping 

to  a quadrilateral. 
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u 


C-i  Surface  Constants  Card  #1 


Corresponds  to  rectangular  dimensions  Lj^  and  L.^  in  STAGS 
manual  for  NSHELL=3,  if  NCART  = 0.  If  NCART  = 1,  and  represent 

maximum  "natural"  coordinates. 


Variable 


Format  Columns 


Description 


PR0P1 


E10.6  1-10 


Length  parallel  to  local  x 
axis.  Use  1 . 0 if  NCART=1. 


PR0P2 


E10.6  11-20 


Length  parallel  to  local  y 
axis.  Use  1.  0 if  NCART  = 1. 


C-3  Surface  Constants  Card  #2 


This  card  is  included  only  if  NCART=1.  It  represents  local  coordinates 
of  the  quadrilateral  branch  corner  points  corresponding  to  the  system  established 
in  the  STAGS  user's  manual. 


Variable 


Format  Columns 


Description 


X(I),  Y(J)  8E10. 6 1-80 

(J-1,  4) 


Local  cartesian  coordinates 
of  branch  corner  points. 
Must  be  defined  in  clockwise 
order  starting  with  the 
coordinates  (0. , 0. ), 


DESCRETIZATION 
D-1  Grid  Definition  Card 


Variable  Format  Columns 

NR  15  1-5 


Description 

Number  of  rows.  NR  positive, 
spacing  is  constant.  NR 
negative,  spacing  is  variable, 
read  in  X or  S values  on 
card  D-2.  ^ 
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Variable 

F ormat 

Columns 

Description 

NC 

15 

6-10 

Number  of  columns.  NC 
positive,  spacing  is  constant. 
NC  negative,  spacing  is 
variable,  read  in  Y or  Tl 
values  on  card  D-  3. 

D-2  X - or  % - Coordinate  Card 

Include  this  card  only  if  NR  < 0 on  the  D-1 

card. 

Variable 

Format 

Columns 

Description 

X(I) 

(1=1,  NR) 

8E10.6 

1-80 

X 

or  1 value  for  row  I. 

D-3  Y-  or 

Tl-  Coordinate  Card 

Include  this  card  only  if  NC  < 0 on  the  D-1 

card. 

Variable 

Format 

Columns 

Description 

V(J) 

(J=l,  NC) 

8E10.6 

1-80 

Y 

or  T)  value  for  Column  J. 

AERODYNAMIC  PARAMETERS 

E-1  Aerodynamic  Analysis  Type  Summary 

Variable 

Format 

Columns 

Description 

NTHK 

15 

1-5 

1 

- Include  thickness  effects. 

NCAM 

15 

6-10 

1 

-Include  camber  effects. 

NCON 

15 

11-15 

1 

- Include  control  surface 

effects . 

NDYN 

15 

16-20 

1 

- Include  fuselage  nose  effects 

MODAL 

15 

21-25 

1 

- Save  input  file  for  sub- 

sequent modal  flutter 
analysis. 

£-2  Flight  Parameter  Card 

Variable 

Format 

Columns 

Description 

ALPHA 

ElO.  6 

1-10 

Angle  of  attack  (deg. ) 

B-5 


Variable 

Format 

Columns 

Description 

PR 

ElO.  6 

11-20 

Ambient  pressure  (psi) 

RHO 

£10.  6 

21-30 

Ambient  density  j 

AMACH 

£10.  6 

31-40 

Representative  Mach  no. 

GAM 

ElO.  6 

41-50 

Ratio  of  specific  ’leats  (=’1.4) 

BREF 

£10.  6 

51-60 

Ref.  semi-chord  (in. ) 

E-2  Section  Profile  Card 

Variable 

Format 

Columns 

Description 

NSYM 

15 

1-5 

I'  - Sy..iirietric  profile 

15 

1 - Nonsymmetrlc  profile 

NPROF 

6-10 

1 - Diamond  profile 

2 - Wedge  profile 

3 - Parabolic  profile 

E-4  Thickness 

/Chord  Ratios, 

Upper  Profile 

Include  this  card  only  if  NTHK=1  on  the  E-1 

card  and  ISTEER^O  on  the 

B-1  card. 

Variable 

Format 

Columns 

Description 

TU(J) 

8E10.  6 

1-80 

Thickness /chord  ratio 

(J=l,  NC) 

at  span- wise  station  J, 
upper  profile. 

E-5  Thickness /Chord  Ratios, 

Lower  Profile 

Include  this  card  only  if 

1)  ISTEER=0  on  the  B-1  card, 

2)  NTHK=1  on  the  E-1  card,  and 

3)  NSYM=0  on  the  E-3  card. 


B-6 


Variable 

Format 

Columns 

Description 

TL(J) 

8E10.  6 

1-80 

Thickness /chord  ratio 

CJ=1,  NC) 

at  span-wise  station  J, 

lower  profile. 

E-6  Local  Chord  Lengths 


Include  this  card  only  if  NCAM=1  on  the  E-1  card. 


Variable 


Format 


Columns 


Description 


CB(J)  8E10. 6 1-80  Local  chord  lengths. 

(J=l,  NC) 


E-7  Control  Surface  Parameter  Card  #1 


Include  this  card  only  if  NCON=l  on  the  E-1  card. 


Variable  Format  Columns 

XH(J)  8E10. 6 1-80 


Description 

Distance  of  leading  edge 
from  hinge -line. 


E-8  Control  Surface  Parameter  Card  #2 

Include  this  card  only  if  NCON=l  on  the 


Variable 

Format 

Columns 

BETA(J) 
(J=l,  NC) 

8E10.6 

1-80 

E-1  card. 

Description 

Control  deflection 
in  radians. 


E-9  Fuselage  Nose  Parameter  Card 


Variable 


Format 


Columns 


Description 


DN 

XN 

CDN 


E10.6 

E10.6 

E10.6 


1-10  Fuselage  nose  diameter 

11 -ZO  Distance  from  leading  edge  to 

Zl-30  Nose  drag  coefficient  body 


B-7 


B.  2.  3 User  Written  Subroutines 


At  present,  a user  written  subroutine  is  available  for  defining  the 
maximum  thickness /chord  ratios  (i.e.,  profile)  as  a function  of  span-wise 
station.  It  is  called  when  ISTEER=1  and  is  referred  to  as  subroutine 
PSTEER: 

SUBROUTINE  PSTEER  (J,  Y,  C,  TU,  TL)  PSTEER.  2 

END  PSTEER.  3 


where: 


{J  = Span-wise  station  (i.e. , column). 

Y = Coordinate  (local)  of  span-wise  station. 
C = Chord  length. 


Returned 


TU  = Maximum  thickness/C  (upper  profile). 
TL  = Maximum  thickness/C  (lower  profile). 


Note:  TL  may  be  neglected  when  the 
profile  is  symmetric  (NSYM=1). 


B.2.4  Output  File  Creation 

Two  peripheral  files  are  created  by  HYPALM: 

1)  The  QTj  matrices  are  written  to  TAPE3  in 
preparation  for  a STAGS  aeroelastic  analysis 
(if  ISTAGS=1  on  the  B - Icard). 
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2)  The  entire  input  file  is  written  to  TAPES 


in  preparation  for  FLAM,  the  post-STAGS 
flutter  analyzer  (if  MODAL=l  on  the  E-1  card). 

It  is  therefore  necessary  to  save  these  files  either  on  disk,  or 
a physical  tape  unit  should  the  user  decide  to  execute  the  program  package 
in  stages  rather  than  within  the  same  runstream. 
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B.3 


Program  HYPATH  (Hypersonic  Aerodynamic  Thermal  Analysis) 


! 

i 

i 

fi 


i 


B.  3.  1 Fiinction 


HYPATH  prepares  for  a STAGS  analysis  by  creating  a temperature 
history  file,  based  on  aerodynamic  heat  convection,  radiation  and  one 
dimensional  conduction  theory,  for  a given  flight  profile.  The  analysis 
includes  all  "aerodynamically  active"  STAGS  branches. 


B.  3.  2 Data  Card  Input 

GENERAL  CONTROL  CARDS 

A-1  Case  Title  Card 

Variable  Format 

Columns 

Title  12A6 

1-72 

B-1  Execution  Summary  Card 


Variable 

Format 

Columns 

ISTEER 

15 

1-5 

ISTAGS 

15 

6-10 

NBRNCH 

15 

11-15 

NCHK 

15 

16-20 

Description 
Case  Title 


Description 

0 - Nodal  properties  to  be 

input  on  cards. 

1 - Nodal  properties  to  be 

defined  in  user  written 
subroutine  STEER. 

0 - Suppress  temperature 

history  file. 

1 - Create  temperature 

history  file. 

Total  number  of  branches  in 
STAGS  model. 

0 - Perform  analysis. 

1 - Verify  input  data  only. 


i 
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Variable 

Format 

Colunuis 

Description 

NROB(I), 

1515 

1-75 

Number  of  rows  in  branch  I, 

IROC(I) 

number  of  columns  in  branch  I, 

(1=1,  NBRNCH) 

and  row-or-column  only  switch: 

0 - Every  node  in  branch  is  to 

be  analyzed. 

1 - Only  one  node  per  row  need 

be  analyzed. 

2 - Only  one  node  per  column 

need  be  analyzed. 

3 - Only  one  node  is  necessary 

to  represent  entire  branch. 
1 - Ignore  branch  in  analysis 
(i. e.,  "inactive"). 


JFMAX 

15 

6-10 

lATM 

15 

11-15 

Number  of  layers  thru  thickness 
to  be  considered.  Maximum  = 6. 
Note:  NLAY  + 1 temperatures 
are  computed.  If  symmetry  is 
later  specified,  2 * NLAY  + 1 
temperatures  are  implicitly 
obtained. 


Number  of  data  cards  in  flight 
profile  deck. 


Atmospheric  data  specification. 

0 - Atmospheric  profile  deck 

will  be  read. 

1 - 1966  U.S.  standard^ 

Atmospheric  data  will  be  used 

2 - Mil.  Std.  210- A Std.  Hot  day 

data  will  be  used. 

3 - Mil.  Std.  210-A  std.  Cold 

day  data  will  be  used. 
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Variable 


Format 


1 


I ) 


JAMAX 


IMORV 


ITIME 


IGRAD 


IPLOT 


Colvimns 

Description 

16-20 

Number  of  data  cards  in 
atmospheric  profile  deck 
(only  necessary  if  IATM=0). 

21-25 

0 - Flight  profile  velocity 

defined  in  teirns  of  Mach 
number. 

1 - Flight  profile  velocity 

defined  in  ft. /sec. 

26-30 

0 - Time  in  flight  profile 

given  in  seconds. 

1 - Time  in  flight  profile 

given  in  minutes. 

31-35 

0 - Initial  temperi.tures  are 

uniform  thru  thickness. 

1 - Initial  temperatures  vary 

thru  thickness. 

36-40 

IPLOT  negative,  temperature 
vs.  thickness  plots  are  required. 
IPLOT  zero  or  positive,  no  plots 
are  required.  (Note:  plots  are 
only  allowed  when  ISTAGS=0  on 
B-1  card. 

B-4  Time  Step  and  Computation  Summary  Card 


Variable 

Format 

Columns 

Description 

TIME  1 

ElO.  6 

1-10 

Initial  time  to  be  considered  in 
analysis.  (Must  lie  within  range 
specified  in  flight  profile  deck, 
card). 

TIME  2 

ElO.  6 

11-20 

Final  time  to  be  considered  in 
analysis.  (Must  lie  within  range 
specified  in  flight  profile  deck). 

NSTEPA 

15 

21-25 

Number  of  time  steps  to  be 
computed  in  finite  difference 
analysis.  The  maximum  allowed 
is  100. 

I 
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Variable 


Format 


Columns 


NSTEPW  15  26-30 


ICOM  15  31-35 


Description 

Nxunber  of  time  steps  to 
be  saved  on  STAGS  temperature 
history  file.  The  maximum 
allowed  is  50.  (In  addition  the 
initial  values  are  always  saved.  ) 

Equation  solution  type. 

0 - Tri-diagonal  solution. 

1 - Cholesky  solution. 


C-1  Plot  Specification  Card 

Include  this  card  only  if  IPLOT  < 0 on  the  card  and  ISTEER=0 
on  the  card. 

Variable  Format  C olumns 

Arg.  XTLvlE(I)  8E10.6  1-80 

(1=1.5) 

D-1  Flight  Profile  Deck  Header  Card 

Variable  Format  Columns 

FLIGHT  12AG  1-72 

D-2  Flight  Profile  Deck 

The  flight  profile  deck  consists  of  JFMAX  cards  (as  specified 
on  the  B-3  card),  each  card  prescribing  the  flight  velocity  and  altitude 
at  a given  time  along  the  trajectory.  Linear  interpolation  is  used  at 
computation  steps  falling  between  the  prescribed  times. 


Description 

Number  of  times  at  which  plots 
are  desired  and  their  values, 
respectively.  The  maximum 
allowed  is  5. 


Description 
Flight  profile. 


Variable  Format 

Columns 

Description 

TF(J),  CF(J),  3E10.6 

1-30 

Time,  velocity  and  eiltitude. 

FF(J) 

respectively,  at  various  points 

(J  = l,  JFMAX) 

along  the  flight  trajectory. 
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E-1  Atmospheric  Profile  Deck  Header  Card 


L 


Include  this  card  only  if  IATM=0  on  the 

Variable  Format  Columns 

ATMOS  12A6  1-72 


B-2  card. 

Description 

Atmospheric  profile  title. 


E-2  Atmospheric  Profile  Deck 

Include  tliese  cards  only  if  IATM=0  on  the  B-2  card.  The 
atmospheric  profile  deck  prescribes  the  ambient  prtjssure  and  temperature 
as  a discrete  function  of  altitude.  JAMAX  cards  (as  specified  on  B-2  card) 
each  representing  an  individual  altitude  are  used.  As  with  the  flight  profile 
deck,  the  independent  variable  (i.  e. , altitude)  must  be  monotonically  increas- 
ing as  linear  interpolation  is  used  for  intermediate  stations. 


Variable 

Format 

Columns 

Description 

TA(J),  CA(J), 

3E10.  fS 

1-30 

Altitude  (ft.  ),  pressure 

FA(J) 

and  temperature  (DEG. 

(J  = l,  JAMAX) 

respectively,  at  JAMAX 
separate  altitudes. 

THE  REMAINING  INPUT  SECTION  IS  REPEATED  ONCE  FOR  EACH 
"ACTIVE"  (I.E.  , IROC  > 0)  BRANCH 


GEOMETRY 

F-1  Shell  Type  Definition  Card 


Variable 

Format 

Columns 

Description 

NSHELL 

15 

1-5 

Shell  type. 

3 - Plate 
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1 Variable 

Format 

Columns 

Description 

j NCART 

t 

15 

6-10 

0 - Rectangular  grid. 

1 - Quadrilateral  grid, 

isoparametric  mapping. 

F-2  Aerodynamic  Flow  Description  Card 

This  card  prescribes  the  branch  configuration  in  terms  of  the 
aerodynamic  problem  class. 


Variable 

Format 

Columns 

Description 

i ISHAPE 

15 

1-5 

1 - Flat  plate 

2 - Cylinder  (not  operational) 

3 - Sphere  (not  operational) 

i IFLCN 

15 

6-10 

1 - Flat  plate,  perfect  gas. 

2 - Flat  plate,  cabrically 

imperfect  gas. 

3 - Stagnation  point,  perfect 

gas. 

4 - Wedge  reattachment 

(cubic  solution) 

5-17  "Reserved  for 
future  use. " 

DIA  ElO.  6 11-20  Diameter  for  cylindrical  or 

spherical  geometry  (reserved 
for  future  use). 


F-3  Surface  Constants  Card  #1 

Corresponds  to  rectangular  dimensions  and  in  STAGS 
manu'l  for  NSHELLi=3,  if  NCART=0.  If  NCART=1,  Lj  and  represent 
maximum  "natural"  coordinates. 
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Variable 


Foirmat 


Columni 


Deecription 


PRQ&Pl 


£10.6  1-10 


Length  parallel  to  local  x 
axle.  Use  1.  0 if  NCART-1. 


PROP2 


£10.6  11-20 


Length  parallel  to  local  y 
axis.  Use  1. 0 if  NCART=1. 


F-4  Surface  Constants  Card  //2 

This  card  is  included  ozily  if  NCART»1.  It  represents  local  coordinates 
of  the  quadrilateral  branch  comer  points  corresponding  to  the  system  established 
in  the  STAGS  user's  manual. 


Variable  Format  Columns 

X(I).  Y(J)  8£10.6  1-80 


DISCRETIZAOOIl' 
G-1  Grid  Definition  Card 


Variable  Format  Columns 

NR  15  1-5 


NC 


15  6-10 


Description 

Local  cartesian  coordinates 
cf  branch  corner  points. 
Must  be  defined  in  clockwise 
order  starting  with  the 
coordinates  (0. , 0. ). 


Description 

Number  of  rows.  NR  positive, 
spacing  is  constant.  NR 
negative,  spacing  is  variable, 
read  in  X or  | values  on 
card  G- 2 . 

Number  of  colximns.  NC 
positive,  spacing  is  constant. 
NC  negative,  spacing  is 
variable,  read  in  Y or  ‘I] 
Values  on  cardC-3  . 
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G-Z  X-  or  Coordinate  Card 


Include  this  card  only  if  NR  < 0 on  the  G-1  card. 


Variable 

Format 

Columns 

Description 

X(I) 

(1=1,  NR) 

8E10.  6 

1-80 

X or  ? value  for  row  I. 

G-3  Y-  or 

T|-  Coordinate  Card 

Include  this  card  only  if  NC  < 0 

on  the  G-1  card. 

Variable 

Format 

Columns 

Description 

Y(J) 

(J  = l,  NC) 

8E10.  6 

1-80 

Y or  T]  value  for  Column  J. 

‘ In  the  event  that  the  steering  subroutine  is  not  used  (i.e. , ISTEER=0 

on  the  B-1 

card,  the  following  cards  are 

required  for  each  node  being  con- 

sidered  within  an  "active" 

branch. 

H-1  Nodal  Information  Card  §1 

Variable 

Format 

CcTumns 

Description 

NSN 

15 

1-5 

Number  of  node  being  considerec 
(must  be  sequential). 

TT 

E10.6 

6-15 

Total  thickness  (in.)  character- 
istic re  no.  distance  (in. ). 

EM 

E10.6 

26-35 

Emissivity 

DELTA 

E10.6 

36-45 

Sweep  angle  (deg.  ) 

ALPHA 

E10.6 

46-55 

Heat  transfer  coefficient 
at  midplane. 

0.  - symmetric  section  only 
one-half  considered. 
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H-Z  Nodal  Information  Card  #2 


Include  this  card  only  if  IGRAD=0  on  the  B-3  card. 


Variable 

TO 


Format  Columns  Description 

£10.6  1-10  Initial  temperature,  uniform 

thru  thickness  (DEG.  F). 


H-3  Nodal  Information  Card  #3 

Include  this  card  only  if  IGRAD=1  on  the  B-3  card. 


Variable 


F o rm  at  C olumns 


Description 


TZER(I)  8E10. 6 1-80 

(1=1,  NLAY  + 1) 


Initial  temperature  distribution 
thru  thickness  (or  to  mid-surface 
if  ALPHA=0.  ontheH-1  card). 


H-4  Nodal  Information  Card  #4 


Variable  Format  Columns 

NMAT(I)  1615  1-80 

(1=1,  NLAY) 


U 

Description 

Material  type  of  each  layer  thru 
thickness  (or  half  of  thickness 
if  ALPHA=0.  on  the  H-1  card). 


1 - 20Z4T  A1 

2 - 701  A1 

3 - Beryllium 

4 - Copper 

5-  SAE  1010  Steel 

6 - SAE  1020  Steel 

7 - SAE  4130  Steel 

8 - 300  Mar  aging  Steel 

9 - 301  Stainless  Steel 

10  - 403  Stainless  Steel 

11  - 410  Stainless  Steel 
12-  G-18-8  Stainless  Steel 

13  - 17-7  PH  Stainless  Steel 

14  - Iconel  - X 

15  - Magnesium  (DCW-HK-31XA-F) 

16  - Tungsten 

17  - Fiberglass 

18  - Avcoat  No  . 1 

19  - Insulcork 

20  - Titanium 

21  - TMC  Propellant 


U 
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B.  3.  3 User  Written  Subroutine 


( 


( 


The  user  has  the  option  of  specifying  all  of  the  nodal  parameters 
appearing  in  data  cards  H-1  to  H-4  within  a single  user  written  subroutine, 
STEER,  In  this  way,  functional  relationships  with  branch,  row  and  column 
numbers,  transverse  station  and  spacial  coordinates  as  independent  variables 
may  be  conveniently  established. 

SUBROUTINE  STEER  (NSN,  IB,  X,  Y,  ETA,  ZETA,  IF  STEER.  2 

1,  TZER,  T,  NMAT,  ALE,  DELTA,  EM,  ALPHA)  STEER.  3 

DIMENSION  TZER(7),  NMAT(6)  STEER.  4 

(INSERT  STEERING  CODE)  STEER.  5 

END 


WHERE: 

NSN  = Number  of  structural  node.  This 

number  is  incremented  consecutively 
through  the  entire  set  of  "active" 
branches.  Depending  upon  the  value 
of  IROC(I)  (see  B-2  card)  it  may  be 
incremented  only  once  per  row,  once 
per  col. , etc. 


Given: 


ilB  = Branch  number 

X = Local  X coordinate  of  node 

y = Local  y coordinate  of  node 

ETA  = "Natural"  x coordinate  (meaningful  only 
when  NCART=1) 
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ZETA  = "Natural"  y coordinate  (meaningful  only 
for  NCART=1). 

IP  = Number  of  points  defined  thru  the  thickness 
(=NLAY+1) 


P^eturned: 


TZER(I) 
(1=1.  IP) 

See  H-1  card 
for  definition 

T 

Equivalent  to 
H-1  card 

NMAT(I) 
(1=1.  IP-1) 

See  H-1  card 

ALE 

See  H-1  card 

DELTA 

See  H-1  card 

EM 

See  H-1  card 

ALPHA 

See  H-1  card 

TT  on 


B,  3.  4 Output  File  Creation 

A thermal  history  file  containing,  essentially,  T=T  (x,  y,  t)  and  q = q(t) 
is  written  to  TAPE 3 when  ISTAGS=1  on  the  B-1  card.  Due  to  the  data  re- 
organization necessary  to  produce  this  file,  it  is  advised  that  a preliminary 
run  with  NCHK=1  is  made  for  model  verification  before  attempting  the  file 
creation  run  with  ISTAGS=1  and  NCHK=0.  As  with  program  HYPALM,  it 
will  also  be  necessary  to  save  TAPE3  on  a peripheral  device  if  the  STAGS 
analysis  is  to  be  performed  in  a separate  run. 
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B.4  The  "STAGS"  Interface 


B.4.1  Data  F ile  Input 


Assuming  the  aeroelastic  pre-processors,  programs  HYPALM  and 
HYPATH,  have  been  executed,  the  following  data  files  will  have  been  created 


1)  Aerodynamic  Load  Matrices 
(written  to  TAPES  by  HYPALM) 

2)  Thermal  History  * 

(written  to  TAPES  by  HYPATH) 

3)  Aerodynamic  Model  Description 
(written  to  TAPES  by  HYPALM) 


As  shown  in  Figure  B.  1,  two  of  these  files  are  usually  required  for 

A 

I a STAGS  aeroelastic  analysis  and  therefore  must  be  copied  to  TAPE18 

(in  the  order  depicted)  before  execution.  In  the  case  of  a "restart"  analysis 
a third  file  must  also  be  copied  to  TAPE18.  This  refers  to  the  "restart" 
file  which  is  automatically  written  to  TAPE18  in  a previous  analysis  and 
which,  of  course,  must  be  "saved"  at  the  time  and  then  copied  to  TAPE18 
(in  the  position  depicted  in  Fig.  B.  1)  before  executing  a restart  run.  Note 
that  files  1)  and  2)  above  do  not  have  to  include  the  entire  flight  trajectory 
at  the  outset,  as  the  files  may  be  re-created  with  programs  HYPALM  and 
HYPATH  to  account  for  different  intervals  in  the  trajectory  prior  to 
successive  STAGS  "restart"  runs. 


It  is  worth  emphasizing  that  TAPE3  refers  to  a local  file.  If  HYPALM 
and  HYPATH  are  to  be  executed  consecutively  within  the  same  ninstream, 
the  contents  of  TAPE 3 must  therefore  be  copied  to  an  incermediate  file 
and  TAPES  "xmloaded"  before  executing  the  second  program  and  "re- 
creating" a new  TAPES. 
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B.4.2  Data  Card  Modifications 

The  data  card  requirements  for  the  various  possible  types  of  STAGS 
analyses  are  described  in  the  STAGS  user's  manual.  Those  particular  entries 
which  are  affected  in  a STAGS  aeroelastic  analysis  are: 


1)  Two  additional  parameters  have  been  included  on  the  B-1  ( ard  - 


Variable 

Format 

Columns 

lAERO 

15 

31-35 

ITHERM 

15 

36-40 

Description 

1 - Aeroelastic  analysis. 

Aerodynamic  load  matrices 
are  to  be  read  from  TAPE18. 
0 - Not  an  aeroelastic  analysis. 


1 - Hypersonic  thermal  — 

analysis  has  been  performed. 
Thermal  history  is  to  be  read 
from  TAPE  18. 


B-2  Aerodynamic  Load  Matrix  Selection  Card 

Include  this  card  only  if  IAERO=l  on  the  B-1  card. 


Variable 

Format 

Columns 

Description 

IQTO 

15 

1-5 

1 - Read  QTO  matrix. 

IQTl 

15 

6-10 

1 - Read  QTl  matrix. 

IQT2 

15 

11-15 

1 - Read  QT2  matrix. 

IQT3 

15 

16-20 

1 - Read  QT3  matrix. 

IPQT 

15 

21-25 

1 - Print  QT  matrices. 

RHOINF 

E10.6 

26-35 

Ambient  air  density. 

B-22 


3)  Note  that  when  performing  a quasi-static  STAGS  aeroelastic 
analysis  INDIC  should  be  set  equal  to  3 for  a basic  nonlinear 
analysis,  and  to  5 if  a subsequent  modal  flutter  analysis  is 
anticipated.  In  the  latter  case  a number  of  vibration  modes 
about  various  pre-stress  states  are  computed  and  saved  as 
described  under  the  F-2  and  F-3  data  cards  in  the  STAGS 
user's  manual. 

4)  As  long  as  the  thermal  history  file  is  to  be  considered  (i.  e.  , 
ITH£RM=1)  the  dynamic  pressure  as  a function  of  time  is 
directly  accessible.  Therefore,  when  performing  a quasi- 
static nonlinear  analysis  the  load  multipliers  on  the  F-1  card 
(system  A)  represexit  actual  values  of  the  time  (in  seconds) 
and  correspond  to  the  flight  profile  deck  used  as  input  to 
program  HYPATH.  When  performing  a transient  response 
analysis  the  flight  profile  time  is  then  identical  with  the 
times  specified  on  the  H-1  data  card.  One  potential 
complication  arises  when  the  hypersonic  thermal  effects 

are  omitted  and  an  aeroelastic  analysis  is  attempted  with 
IFLU=1  and  1TH£RM=0  on  the  B-1  card.  In  this  case  the 
consequences  are  as  follows: 

a)  For  a nonlinear  static  analysis,  since  the  dynamic 
pressure  history  is  no  longer  accessible,  the  load 
multipliers  on  the  F-1  card  represent  actual  values 
and  increments  of  the  dynamic  pressure  in  psi. 
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b)  For  a transient  analysis,  the  time  parameter  is 
specified  as  usual,  but  the  flight  velocity  must 
be  specified  as  a function  of  time  in  subroutine 
VINF. 

B.4.  3 Output  File  Creation 

Potentially,  twc  data  files  are  created  during  a STAGS  aeroelastic 
analysis: 

1)  A restart  file  is  written  to  TAPE  18  in  all  cases. 

2)  A vibration  mode  file  containing  displacement  mode  vectors, 
natural  frequencies  and  generalized  masses  for  each  of  the 
selected  pre-stress  states  is  written  to  TAPE9  if  INDIC=5  and 
IAERO=i  ontheB-1  card. 

These  files,  of  course,  must  be  saved  on  a peripheral  device  for 
use  in  a subsequent  execution  of  STAGS  and/or  the  flutter  post-processor 
(FLAM). 

B.  5 Program  FLAM  (Flutter  Analysis  Modal) 

B.5.1  F unction 

FLAM  performs  a linearized  modal  flutter  analysis  about  a series 
of  nonlinear  pre- stress  states  along  a flight  trajectory,  incorporating  the 
vibration  modes  computed  by  STAGS  at  these  prescribed  states.  Internally, 
program  GALA  employs  the  displacement  mode  vectors  to  obtain  generalized 

The  da.  for  each  pre-stress  state  is  separated  by  an  END-OF-FILE 
mark,  which  must  be  considered  when  copying  TAPE9  from  one  logical 
file  to  another. 
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aerodynamic  forces  which  are  utilized  (in  addition  to  the  natural  frequencies, 
generalized  masses  and  selected  "reduced"  frequencies)  by  program  MOFA  to 
solve  a reduced  (complex)  eigenvalue  problem  and  hence  predict-  the  critical, 
i.  e.  , flutter,  speed  in  an  iterative  manner. 

B.5.  Z Data  File  Input 

Two  data  files  are  prerequisite  to  performing  a FLAM  analysis,  in 
addition  to  a minimal  number  of  data  cards: 

1)  "Model  description"  file  created  by  program  HYPALM  previously 
on  TAPES.  This  file  must  be  copied  (or  attached)  to  TAPEl 
before  FLAM  execution. 

2)  "Vibration  mode"  file  created  by  STAGS  in  a previous  analysis 
and  containing  eigenvectors,  eigenvalues  and  generalized  masses 
for  a number  of  modes  evaluated  at  a number  of  pre-stress  states 
(i.e.  , times  along  flight  trajectory).  This  file,  originally  saved 
on  TAPE9  must  be  copied  (or  attached)  to  TAPE2  before  FLAM 
execution.  Note:  Each  pre-stress  state  constitutes  an  individual 
sub-file  on  TAPE9  and  hence  more  than  one  "copy"  command 
may  be  required  in  transferring  the  contents  to  TAPE2. 

B.5.3  Data  Card  Input 

A-1  Case  Title  Card 

Variable  Format  Columns  Description 


TITLE 


12A6 


1-72 


Case  title. 


aerodynamic  forces  which  are  utilized  (in  addition  to  the  natural  frequencies, 
generalized  masses  and  selected  "reduced"  frequencies)  by  program  MOFA  to 
solve  a reducv^d  (complex)  eigenvalue  problem  and  hence  predict  the  critical, 
i.  e. , flutter,  speed  in  an  iterative  manner. 

B.5.  2 Data  File  Input 

Two  data  files  are  prerequisite  to  performing  a FLAM  anslysis,  in 
addition  to  a minimal  number  of  data  cards: 

1)  "Model  description"  file  created  by  program  HYPALM  previously 
on  TAPES.  This  file  must  be  copied  (or  attached)  to  TAPEl 
before  FLAM  execution. 

2)  "Vibration  mode"  file  created  by  STAGS  in  a previous  analysis 
and  containing  eigenvectors,  eigenvalues  and  generalized  masses 
for  a number  of  modes  evaluated  at  a number  of  pre>stress  states 
(i.e.  , times  along  flight  trajectory).  This  file,  originally  saved 
on  TAPE9  must  be  copied  (or  attached)  to  TAPE2  before  FLAM 
execution.  Note:  Each  pre-stress  state  constitutes  an  individual 
sub-file  on  TAPE9  and  hence  more  than  one  "copy"  command 
may  be  required  in  transferring  the  contents  to  TAPE2. 

B.5.3  Data  Card  Input 

A-1  Case  Title  Card 

Variable  Format  Columns  Description 


TITLE 


12A6 


1-72 


Case  title. 


B-1  Case  Summary  Card 


Variable 

Format 

Column 

Description 

NLD 

15 

1-5 

Number  of  pre-stress  states 
(i.  e.  , points  along  flight 
trajectory)  at  which  a flutter 
analysis  is  desired. 

B-2  Load/ Time 

Step  Specification  Card 

Variable 

Format 

Column 

Description 

PLD(I) 

(1=1,  NLD) 

8E10.  6 

1-80 

Times  (or  load  multipliers) 
denoting  pre-stress  states 

to  be  evaluated.  Vibration 
modes  from  STAGS  analysis 
for  each  must  be  present 
on  TAPE2.  It  is  not  necessary  j j 
to  request  all  of  the  pre-stress  ^ 
states  resident  on  TAPE 2,  but 
PLD(I)  must  be  requested  in  a 
monotonically  increasing  order. 


THE  REMAINING  CARDS  ARE  REPEATED  NLP  TIMES. 

C-1  Reduced  Frequency  Summary  Card 

For  each  pre-stress  state  a modal  flutter  analysis  is  actually  performed 
for  each  of  a number  of  selected  "reduced"  frequencies  and  a tentative  flutter 
speed  is  obtained  graphically  by  means  of  an  interpolative  procedure  described 
elsewhere  in  the  text. 


Variable 

Format 

Columns 

Description 

NREF 

15 

1-5 

Number  of  "reduced"  frequencies 
to  be  considered.  The  maximum 

currently  permitted  is  8.  ^ 

B-1  Case  Summary  Card 


Variable 

Format 

Column 

Description 

NLD 

15 

1-5 

Number  of  pre-stress  states 
(i.e.  , points  along  flight 
trajectory)  at  which  a flutter 
analysis  is  desired. 

B-2  Load/ Time 

Step  Specification  Card 

Variable 

Format 

Column 

Description 

PLD(I) 

(1=1,  NLD) 

8E10.  6 

1-80 

Times  (or  load  multipliers) 
denoting  pre- stress  states 

to  be  evaluated.  Vibration 
modes  from  STAGS  analysis 
for  each  must  be  present 
on  TAPE2.  It  is  not  necessary  j 
to  request  all  of  the  pre-stress 
states  resident  on  TAPE2,  but 
PLD(I)  must  be  requested  in  a 
monotonically  increasing  order. 


THE  REMAINING  CARDS  ARE  REPEATED  NLD  TIMES. 


C-1  Reduced  Frequency  Summary  Card 

For  each  pre-stress  state  a modal  flutter  analysis  is  actually  performed 
for  each  of  a number  of  selected  "reduced"  frequencies  and  a tentative  flutter 
speed  is  obtained  graphically  by  means  of  an  interpolative  procedure  described 
elsewhere  in  the  text. 


Variable 

Format 

Columns 

Description 

NREF 

15 

1-5 

Number  of  "reduced"  frequeni  ies 
to  be  considered.  The  maximum 
currently  permitted  is  8. 
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c 


C-2  Reduced  Frequency  Specification  Card 


Variable 


Format  Column 


Description 


REF(I)  8E10. 6 1-80 

(1=1,  NREF) 


Reduced  frequency  values 
(rad. /sec. ) . 


